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Summary 

The possibility of predicting the ultimate static bearing 
capacity of a pile from its dynamic behavior during driving is 
investigated by the wave equation method of pile driving analysis. 

The method developed by E. A. L. Smith over a period of years is 
followed, but its validity is first shown by proving that the math- 
ematical model used is equivalent to the wave equation. Corputer 
programs developed far purposes of this investigation are presented 
in their entirety together with detailed instructions on their use 
so that they may be utilized by future investigators who may not be 
familiar with computer work. 

The soil engineering aspects of the problem are explored. The 
information available on dynamic soil properties from published test 
data is reviewed to obtain information applicable to pile driving. 
Using values thus deduced and the computer programs, computations of 
ultimate resistance are made from 2k published pile driving records 
and correlated with their load test results. 

Although not conclusive because of the relatively small num- 
ber of correlations attempted, the results are very encouraging. It 
appears that the wave equation method of pile driving analysis may 
become an accurate tool in predicting a pile's static bearing capacity 
from its driving record. This appears to be the case regardless of 
the type or size of pile and driving equipment, and indications are 
that predictions may be possible for piles driven in cohesive as well 
as granular soils. General limitations of the method are summarized. 
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Values deduced for point and side damping and ground quake 
in various soils are presented. Corresponding values of friction 
acting on the pile sides as a percentage of the ultimate resistance 
are given for the piles investigated. More research into this problem 
is indicated, and recommendations are made for such work. 

It is concluded that even with inexact knowledge of dynamic 
soil properties the wave equation method of pile driving analysis 
appears to give good results. 
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Chapter I 
INTRODUCTION 

The possibility of relating the behavior of a pile during 
driving to its static load carrying capacity by the use of the wave 
equation method of pile driving analysis is encouraging. Many 
engineers today realize that the numerous dynamic pile driving 
formulas which have been used, and are still being used, have seri- 
ous limitations and cannot be depended upon to give reliable results, 
but they are still used for lack of an adequate substitute. Some 
feel that the dynamic formula approach is not valid, and have turned 
instead to trying to determine the bearing capacity of a pile from 
purely static soil mechanics and structural considerations for the 
cases where the pile-soil system is required to carry only static 
loads. Such an approach is certainly valid, and it requires a thorough 
investigation of the soil conditions, together with a careful evalua- 
tion of their engineering properties, and the use of a great deal of 
personal judgment in its application. 

It is believed that the wave equation method of analysis of 
the dynamic pile driving problem is an equally valid approach and may 
offer advantages over the static approach if it can be successfully 
applied. This is not to say that one approach should be used to the 
exclusion of another, but that they are working toward a common goal 
and should supplement each other. The wave equation method is not 
just another dynamic pile driving formula, but, rather, it is a 
method of analysis which is well founded mathematically. It requires 
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a knowledge of the soil and its properties, both static and dynamic, 
for its successful application. Current knowledge of soil dynamics 
is incomplete, but an increasing amount of work is being done in this 
field . 

Approaching the problem from the dynamic side, the over-all 
consideration is to relate the dynamic behavior of the driving equip- 
ment-pile-soil system to the static behavior of the structure-pile- 
soil system with due consideration to changes in properties of the 
soil during and after driving of the piles. A part of this broader 
problem is to relate the dynamic resistance to driving of a single 
pile under the last hammer blow to its static bearing capacity, with 
time effects between driving and static loading minimized. It is for 
this more limited problem that the wave equation method of analysis 
is useful. 

The wave equation alone will not provide a solution to the more 
general consideration which includes the group effect of piles, de- 
velopment of negative friction, consolidation of a clay layer beneath 
the pile tip, long-term changes in water table, deterioration of the 
pile due to insect attack or deleterious chemicals, set up or relaxa- 
tion of the soil with time, or development of hydrostatic uplift. 

That it does not provide a solution to such problems may be obvious, 
but the statement is made here in an effort to clear up the type of 
misunderstandings and objections already voiced in opposition to the 
application of the wave equation to pile driving analysis. In short, 
it is not a magical formula, but simply a tool which is well founded 
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mathematically and is available to assist in visualizing and evaluating 
the dynamics of the problem. 

The wave equation can be utilized to advantage for investigation 
of other facets of pile driving such as stresses in the pile and 
selection of appropriate pile driving equipment for a particular set 
of field conditions; but the main points of this investigation will 
be to consolidate the wave equation theory as applicable to pile driv- 
ing, make readily available computer programs to facilitate the work 
of future investigators, explore the unsolved problems of interaction 
between soil and pile, attempt correlation between wave equation solu- 
tions and a few pile driving records and load tests, and to suggest 
possibilities for further research work in this general area. 
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Chapter II 
BACKGROUND 

Although use has been made of piles for at least 2,000 years, 
attempts to determine the bearing capacity of a pile from its driving 
record have been generally unsuccessful in spite of the amount of 
effort spent in this direction. An extensive treatment of the develop- 
ment and interrelatedness of dynamic formulas by Chellis (ref. 5) is 
included, and he shows that they may be grouped as follows: 

The first dynamic formula was proposed by Major Saunders in 
1851 and is perhaps the simplest and most direct. He equated the 
weight of the ram multiplied by the stroke, with the driving resistance 

multiplied by the set, and applied a factor of safety of 8. 

Another group of proposed formulas contained a fixed coefficient 
which was supposed to compensate in some degree for factors affecting 

the results which were not included as terms in the formulas. The 

Engineering News. Wellington, Vulcan, and Bureau of Yards and Docks 
(not now used by Navy) formulas are of this type. 

Another group of formulas attempts to account for the variables 
by using expressions for efficiency of applied energy by including 
relative weights of the pile and the hammer. The Dutch, Ritter, and 
Benabencq formulas are of this type. 

Other similar formulas try to include the effect of variables 
by using both fixed coefficients and expressions for the relative 
weights of the pile and ram. The Eytelwein and Navy-McKay (not now 
used by the Navy) are of this type. 
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The next group contains either all or part of a series of 
terms designed to represent impact losses in the driving cap, soil, 
and pile during driving. The Redtenbacher and Hiley formulas are 
among those which fall into this category. 

Only a few of the formulas proposed have been mentioned. 

Chellis (ref. $) lists 38, while it is understood that the Editors of 
Engineering News Record have on file b$0 such formulas. Some of these 
formulas are extremely sinple, while others are quite conplicated. 

Some try to approximate the dynamic situation, while others are con- 
trived from a purely statistical approach without reference to the 
numerous variables actually involved. For example, Marvin Gates has 
proposed (ref. lH) that for a pile driven with a Vulcan Number 1 Hammer 
the ultimate bearing capacity may be determined by multiplying U8 by 
the log of (10/set in inches). He claims, ". . .this relationship 
gives more consistently accurate results than the most complex dynamic 
formula yet advanced," in spite of the fact that it contains but two 
simple parameters. That this can be true for a formula which seems 
to ignore the basic aspects of the problem is not surprising when one 
considers the shortcomings and omissions of the other dynamic formulas 
proposed. 

Some of the dynamic formulas were based on Newtonian theory 
of impact which has been shown by A. E. Commings (ref. 10) as in- 
applicable to the pile driving problem. Other formulas subtracted the 
same energy losses twice, as also shown by Cummings (ref. 10). The 
Committee on the Bearing Value of Pile Foundations of the American 
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Society of Civil Engineers recognized the many limitations of dynamic 
pile driving formulas (ref. 8). After this point of agreement, en- 
gineers seem to be divided into two general schools of thought. One 
group is optimistic and feels that rational use can be made of such 
formulas. The other group stresses the limitations of the formulas, 
and basically takes the position that, since none of them are reliable, 
you may as well save yourself work and choose a simple formula. The 
latter school of thought may seem rash at first reading; but it does 
have some merit, for it is simply recognizing that all of the dynamic 
formulas are over-simplifications of a very complex problem. 

It is certain that each of the 1;5>0 formulas available to 
choose from vri.ll give reliable results under particular combinations 
of driving equipment and soil conditions for a particular pile, but 
there would seem to be an almost infinite number of combinations of 
the variables involved making the choice of a formula a difficult one. 
That so many formulas have been developed is indicative of the interest 
in the problem and the ingenuity of engineers in attempting to estab- 
lish a relationship between static and dynamic resistance of piles. 
Empirical formulas certainly have their place when their range of 
applicability is clearly established and when their limitations are 
kept clearly in mind; however, this has not always been the case in 
pile driving. 

In spite of the interest exhibited, there has not been a great 
deal of work done theoretically or experimentally in defining or 
attacking the fundamental problems involved. The complexity of the 
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problem, the difficulty in isolating the variables, and the limited 
techniques and equipment available for detailed analysis have generally 
limited advances in this field. Full utilization, however, of the 
facilities available has not been realized. 

The basic problem is to relate the dynamic resistance to 
driving of a single pile to its static bearing capacity with time 
effects between driving and static loading minimized. Available to 
assist in the analysis of this problem is the wave equation which was 
developed over a hundred years ago by De Saint Venant and Boussinesq 
for end impact on rods. It is well founded mathematically. D. V. 
Issaccs, in 1931, was the first to point out that wave action occurred 
during the driving of piles. In 1938, E. N. Fox published a solution 
to the wave equation applied to pile driving, but, as no electronic 
computers were available at that time, he made a number of simplifying 
assumptions. Again in 19U0 and 19U1, A. E. Cummings reported on the 
work of the foregoing investigators. Ten years elapsed before the wave 
equation was again proposed as being applicable to pile driving analysis, 
this time by E. A. L. Smith (ref. 27). Then Tn 1961 Smith published a 
more comprehensive article on application of the wave equation to pile 
driving and described a numerical approach together with a description 
of his efforts in developing a computer program for the solution (ref. 
28). His efforts have been untiring in developing this approach and 
have opened the door for much needed work in this field. 

The wave equation, as will be shown, mimics the behavior of 
a pile and provides a good indication of what will happen to the pile 
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under conditions of impact for specific boundary conditions. For this 
reason, a distinction is made between the wave equation and dynamic 
formulas which appear to be mathematical but are in reality empirical. 
The present state of knowledge up through Smith's work can provide 
a prediction of effects caused by a ram hitting a pile. Still to be 
investigated is the nature of the effects caused by interaction be- 
tween the pile and the soil. This, too, is a dynamic problem, and 
it lies in the new field of soil dynamics. This matter will be ex- 
plored, but, first, the wave equation theory and its numerical solution 
will be developed. 



M* 
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Chapter III 

DEVELOPMENT OF BASIC EQUATIONS AND THEORY 



A. Theory 

A. E. L. Smith, over a period of several years, derived equa- 
tions suitable for numerical calculations and computer adaptation 
for pile driving analysis by the wave equation. His derivations are 
based on a mathematical model using weights and springs. This work 
has generally been accepted as correct and as equivalent to the wave 
equation,, but the full derivation establishing this equivalence has 
not been published. In order to establish this equivalence and the 
validity of the numerical method, a step-by-step derivation from first 
principles is presented. For completeness. Smith's derivations from 
the mathematical model are included. 

The general plan for establishing this link is as follows: 

a. A partial differential equation is developed and is shown 
to be a form of the wave equation. 

b. This equation is converted into a difference equation for 
ease of numerical solution. 

c. To simplify further numerical solution, a mathematical 
model consisting of weights and springs is assumed as being equivalent 
to the elastic rod used as a basis for steps a and b above. Five 
equations are developed from this mathematical model which is based 

on the assumption that all springs are perfectly elastic, and the pile 
is represented typically as shown by figure (3). 




' 

' *• 1 ’ 4 - i • • • ' 1 U C.V M _ , .a 4 

■ r ‘i T. 



. 

* 

■ 



• r r 



• 


• 


.. 


I . •! * * 


- 




1 






' ' 


T • 

... 

« 


1 Iff 














VI vft 








• 





' . 






10 



d. It is then shown that the equations of step c. can be com- 
bined to form the same equation developed in step b. above, and it is 
concluded that the mathematical model and the five equations developed 
from it are equivalent to the elastic rod, and the wave equation is 
applicable to it. 

Beginning with development of the partial differential equation, 
as shown by Timoshenko and Goodier (ref. 32), consider the longitudinal 
waves in a prismatic bar of constant cross section such as shown in 
figure (l). The axis of the bar is taken as the x-axis. It is assumed 
that: 



a. Cross sections of the bar remain plane during deformation. 

b. The unit elongation at any cross section mn, due to a 

longitudinal displacement D, is equal to , and the corresponding 

i> X 
A D 

tensile force in the bar is AE , where A is the cross sectional 
area and E is Young’s modulus. 

c. A simple tension in the x direction and the unit elongation 

4^ is accompanied by lateral contraction of the amount v ^ where v 
dX dX' 

is poisson’s ratio. Inertia forces corresponding to lateral particle 
motion are neglected since the length of the waves is large in comparison 
with the cross sectional dimensions of the bar. 

The foregoing leads to the equation for force at any section 
in the bar of uniform cross section: 

F = EA . 

The net force on a section of bar dX, figure (l), is: 
a P = Fj - P X > (Ej.ll) <ffi) - EjA-jg - Al g (E x g) dX, 
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assuming that E may not be constant for purposes to be shown below. 
Considering the motion of the particle dX, and applying Newton's second 
law that a particle acted upon by an unbalanced force system has an 
acceleration in line with and directly proportional to the resultant 
of the force system. 



"\ 2 

A (E dX = A dX y - — |, simplifying 

0 "fc 



^ t 2 ^>X ^ SX^ 



( 1 ) 



and if E is constant. 



c> 2 P „ _E_ c) 2 P 
c> t 2 f bx 2 

2 E 

Introducing a constant c of a value such that c = , this reduces 

to the familiar form of the wave equation, ref. (lU) : 



£> 2 P = c 2 b 2 P 

bt 2 bx 2 



Thus, it is established that equation (l) is a form of the wave equa- 
tion. It is desired to modify it to account for boundary conditions 
expected in a pile. This can be done in the general case as follows: 

? fv- < E ti> - K < (2) 

with R representing resistance to driving. Equation (2) can be con- 
verted into a difference equation for numerical solution by standard 
methods (ref. 2U)s 
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^t 2 ^ D U,t+/At) " 2D (x,t) + D (x,t-At)-1 



AAX L 



2 i[ K (x+^i)^ D (x + AX } t) " D (x,t)^" t K (x .Ax)K D x,t D (x- Ax,t) 



KL 

inhere E “ -j- and D(x, t) is the displacement of point x during time 
interval t. 

It is convenient at this point to introduce a simplified 
terminology which will facilitate comparison with model equations to 
be developed later, and to correlate with existing literature on this 
subject (ref. 28). 



^(x, t+ ^t) 


— 


D m 


D (x,t) 


= 


d 

m 


D (x, t- At) 


r= 


d* 

m 


K (x + ^£) 


= 


K 

m 


K (X-^) 


s= 


K n 
m-1 


D (x+Ax, t) 


= 


d m+l 


®(x-Ax, t) 


S3 


d m-l 


R 


= 


R m 



Then 

? 



Lft 



o [D - 2d + d~l = — A— >K (d x _ 

2 L m m m J A £0(2 m m+1 



- d ) - K Ad - d ) - R 



£xt 
Simplifying, 
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m-1 m m-1 



m 
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If the acceleration of gravity, g, is expressed in feet per second, 
then 

E m ” ^m + ^m-l ” ^m^^rn-l “ “ c *m+l )^m " (3) 

For purposes of numerical calculation, it is more convenient to use a 
mathematical model consisting of weights and springs such as is shown 
in figure (2). The following notation will be used and is consistent 
with that used in equation.(3) above and in reference (28). 

1, 2, 3, • . • m-1, m, . . . are subscripts designating par- 

ticular weights and associated springs, 

1, 2, 3, . . . n-1, n, . . . are subscripts designating 

specific time intervals. Zero designates the initial instant. 

The following apply to any time interval n : 

D m = displacement of weight m measured from its initial 

position in inches 

C m = compression of spring m in inches 
F m = force exerted by spring m, in pounds 

= net force acting on weight m, in pounds 
= velocity of weight m, in feet per second 
The following apply to the preceding time interval n-1 : 

= displacement of weight m measured from its initial 
position in inches 

c^ *« compression of spring m, in inches 
f m = force exerte.d by spring m, in pounds 
v m = velocity of weight m, in feet per second 
z m = net force acting on weight m, in pounds 
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The following are usually constants : 

W m = magnitude of weight m, in pounds 

** spring constant for spring m, in pounds per inch 

R m = external force or resistance acting on weight m, 
in pounds 

g = acceleration due to gravity (32.2 feet/second) 

E = Young's modulus of elasticity in pounds/square inch 

A = cross sectional area of pile in square inches 

n 

= unit length in inches 

A t = time interval used for numerical calculations, 
in seconds 

T m = critical time interval for spring m, in seconds 

T.= minimum value of T , in seconds 
mxn m-’ 

0 = over-all "factor of safety" = t 

0jnin = individual factor of safety for spring m = T m / A.t 

The numerical calculations will be a step-by-step process in 

which the five variables, D_, C , F , Z , and V , will be calculated 

9 m* m’ m’ m J m’ 

for each weight or spring in each successive time interval. Formulas 

for these five variables are developed by Smith (ref. 28) as follows: 

- dm + v m (12At) (U) 

The coefficient 12 is included since D m and d are expressed in inches 

mm 

and v in feet per second, 
m 

Considering next, and referring to figure (3), let the dashed 
squares represent the initial positions of weights m and m + 1, and let 
the solid squares represent their positions at the end of time interval 
n. The initial length of spring m is il- y and at the end of time interval i 
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it is £. 1 . D and D m+ -| are the displacements of weights m and m + 1 



m+1 



respectively, during time interval n. Then, 

c m - H -l' 



but 



or 



therefore, 



; - Vi ♦ » 



m 



Vi 



C = D - D 
m m m+1 



(5) 



The force F m is directly proportional to C m and to the spring 



constant K^, so that 



F m = C K 
m mm 



(6) 



From figure (2), weight m is acted upon by springs m - 1 and 
m and by the external force R m . The accelerating force is given by 



Z = F , - F - R 
m m-1 m m 



(7) 



From Newton's second law that change in velocity is equal to 
force times time divided by mass, 

( 8 ) 



V - v + Z (-—S) 
m m m W m 7 



It remains to be shown that the model of figure (2) and equa- 
tions (i») through (8) are equivalent to the wave equation. This is 
done by using equations (l) through (8) in the following manner: 

Rewriting equation (8) for the previous time interval gives 

v m " K + } (8a) 

Similarly, equation (7) for the previous time interval is 

z m “ ^m-1 ” ^m " ^m 



(7a) 
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Substituting equation (7a) into equation (8a) gives 

v = v~ + (f , - f - E ) 'Qy — » 
m m v m-1 m m 7 w m 

Substituting equation (9) into equation (U) gives 

^ a ^ + C v m + ( f m-l “ f m " %) ^T 2 ] 12 At 



(9) 






or 



m 



Om + 12 ^ + ^ [ f m-l - fm ' "ml 

m 



(Ua) 



But equation (U) for the previous time interval is 



d*' + v* (12 At) 
mm 



or 






v* s — — 21 
m r; At 



(Ub) 



Substituting (Ub ) into (Ua) yields 

D- - <*» * 12 -t[^] . [f n .! - f m - H.] 



m 



m 



Substituting equations (5) and (6) for the previous time interval 

into the above equation and simplifying gives 

2 

^m = “ “^m + w ~ [(^-1 " ^^^m-l “ ( ( *m ” ^m+l^^i “ 

m 

which is equation (3). Thus equations (U) through (8) are equivalent 
to the wave equation, and the mathematical model of figure (2) has 
been shown to be a valid one. 

Of course, solving the partial differential equation by a 
numerical integration process is an approximate one, but the results 
should be accurate within five per cent (ref. 28) provided that the 
time interval chosen is small enough and is coordinated with the length 
of pile selected to be represented as a block and spring in the mathe- 
matical model so that a stable set of calculations is produced. This 
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matter will be discussed in more detail later, but experience indicates 
that 5- to 10-foot sections of a pile may be so represented with a time 
interval of from 0.00025 to 0.00033 seconds, depending upon the type 
of pile material. As to the numerical method, it may be considered 
sufficiently accurate in light of present knowledge upon which the 
many factors influencing any pile driving operation are evaluated. 

B. Mathematical Model Representation 

That a pile may be treated mathematically as a series of 
weights and springs not only facilitates the computations but also 
assists in visualizing the problem. For that reason, the remainder of 
this paper will deal principally in terms of the mathematical model of 
a pile. It is, therefore, important that this method of pile repre- 
sentation be thoroughly understood at the outset. The following sec- 
tions will adhere closely to reference (28). 

The Hammer -Ram is usually a short, heavy, rigid object which 
may be represented by an individual weight without elasticity. In 
figure (U) the first weight represents the ram. 

The velocity of the ram at the instant of impact is needed to 
start the numerical calculations. Ordinarily the rated foot pounds of 
energy of the hammer is given by the manufacturer. The efficiency is 
sometimes given, and sometimes must be assumed. From these data the 
velocity may be calculated at impact by the following formula: 



velocity at impact /Rated energy (ft-lb) x efficiency x 6U.4 
(ft/sec) 53 V weight (lb) 



( 10 ) 
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The Capblock is a short, springy object of wood, plastic, or 
similar material which is comparatively light, and which may there- 
fore be represented by a spring. In figure (U) this spring is 
identified as K^. The form of the stress-strain diagram (or the 
hysteresis loop) that is produced as the capblock is suddenly com- 
pressed and then allowed to re-expand is assumed to be as shown by 
figure (5). Compression occurs along the line AB whose slope is 
determined by the elastic constant of the capblock Restitution 
occurs first along the line BD, and then, because the capblock cannot 
transmit tension, it is completed along line DA. ABDA is the hyster- 
esis loop. The computations are made so that 

Area BCD / \2 Energy output 

Area ABC "" ^ e l' ” Energy input 

where is the coefficient or restitution of the capblock. This 

relationship may be derived in the following manner: 

Let Vf „ be the weight of the ram 
ram 



W. 



r ca p be the weight of the capblock 



’ram 



be the velocity of the ram before impact 



v cap 138 the ve l°oity of the cap before impact 
be the velocity of the ram after impact 
be the velocity of the cap after impact 



ram 



cap 

Then, from Newton's laws, 



W v + W v 
ram v ram cap cap 



w ram v ram + ^cap v cap 



and by definition, 
e l 



v* -v 1 

cap ram 

v -v 
ram cap 
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If 



V = v 

cap v cap 



= 0 



then 



'ram 



ram 



2 2 

energy output s = (©i v ram) 2 

energy input “ * (W * ’ * 

From figure (5), 



area BCD 
area ABC 



= e 



energy ' output 
energy input 



( 11 ) 



During compression the spring force may be computed from equation 

(6) > , > 

F i - Vi (12> 

Inelasticity of the capblock must be considered during restitution, 
however. From equation (11) and figure (5), 
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. (1/2) BC x CD 


e CD 
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(1/2) BC x AC 


AC 


but, 


AC 
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= e 2 C 





1 1 max 

by similar triangles ADE and BCD, 
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Writing the equation for line DB gives 

1-e? 

7 



F - S C - BC C — 2^) 
A CD ° 



and substituting for CD, 



BC 



1 e 2 C 

1 1 max 



Ct - 



* °1 



BC (IS) 

er 

1 



Since BC = 



^1 max % ’ 



F - — _ c ) 

1 2 _ 1 max l v -=,2 7 



e ? C 
1 1 max 



which reduces to 



K 1 C 1 



•* c, 



1 max 






- 1) 



(13) 



Calculations involving the capblock will, therefore, utilize 
alternately equation (12) during compression and equation (13) during 
restitution. C^ max should not be thought of as a constant, but as 
the maximum compression occurring in spring before restitution be- 
gins. It often happens that recompression occurs and a new value of 
Cf max is reached before restitution begins again. This is illus- 
trated by figure (5A) where the stress-strain diagram begins at A, 
goes to B during compression following equation (12), goes to C 
diring restitution as governed by equation (13), at C recompression 
begins, goes to B following equation (13), continues to D by equa- 
tion (12), and thence to E by equation (13). As previously described 
the return to A is along EA because no tension can be transmitted. 
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That rec expression can occur may not be apparent readily, but a 
feeling for it may be acquired in solving a problem or two manually. 

As will be noted later, equations (12) and (13) are used in 
the case of a cushion block, but with the subscripts changed to 2. 

Few tests have been made to determine the elastic character- 
istics of capblocks under impact conditions. Raymond International, 
Inc., in conducting a small number of tests found that the character- 
istics of a wooden capblock vary during driving, and concluded that 
in order to be on the conservative side in computing pile penetration 
per blow, a hardwood capblock with the grain vertical and of an 
original six-inch height, may be assumed to have the following 
characteristic s : 

Spring constant K = 20,000 A (lb/inch of compression) 

where A is the horizontal cross sectional area in square inches. The 
coefficient of restitution was 0.5. Ihe tests also showed that a 
Micarta capblock (Nema Grade "C" ) of 12 inches 1 height had a spring 
constant given by the following relation: 

K = 1+5,000 A 

and the coefficient of restitution was 0.8. 

Pile cap, or Follower, or Helmet — like the ram of the hammer, 
the pile cap is ordinarily a short, heavy, rigid object that can be 
represented by a single weight without elasticity, such as Wg in 
figure (1+). If the pile cap is long and slender, as is the case when 
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used as a follower to drive piles below water or the ground, it is 
better to represent it by a number of weights and springs as shown 
in figure (6e). In this case elastic constants must be computed by 
the formula 



K 



AE 



(HO 



where A is the cross sectional area in square inches, E is the 
modulus of elasticity in pounds per square inch, and X is the unit 
length in inches represented by one spring. 

Cushion Blocks — (called "head packing in the U.K.) — in figure 
(U) springs to inclusive represent the elasticity of the pile 
itself. However, if a precast concrete pile is being driven, a cushion 
block must be used under the pile cap Wg to protect the concrete from 
shattering. In this case, spring K 2 would represent the cushion in 
combination with the first spring of the pile. Figure (6c) applies 
to the cushion block as well as the capblock, and they are treated 
mathematically in the same manner. Dynamic tests similar to those 
described for capblocks show that U-inch- thick pine boards with the 
grain horizontal as used on top of a precast pile to distribute the 
blow evenly may be assumed to have a spring constant equal to 3,!i80 A 
and a coefficient of restitution of 0.5. Of course, the character- 
istics of the wood vary during driving, but use of the above 
relationship is on the conservative side when used for calculating 
penetration of the pile hammer blow. 

The Pile will be either of wood, concrete, or steel. All of 
these materials have traditionally been treated structurally as behaving 
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elastically within the working load ranges. Because the pile is 
compressible due to its length, it is subject to wave action under the 
blow of a hammer. This wave action may be analyzed by dividing the 
pile into short unit lengths of $ to 10 feet. The weight of each 
length is noted as through W-^q, and the elasticity of each length 
by an individual spring as shown by K 2 through in figure (U). If 
the pile is of uniform section, the weights and springs representing 
it are identical. If the pile is tapered or is a composite pile, or 
is in any way of irregular cross section, then the weights and springs 
are modified to represent the pile under investigation. That com- 
paratively large unit lengths, such as 10 feet, give desired accuracy 
in computing the action of impact waves may not be clear. Although 
water waves are transverse waves and the waves in a pile are longitudinal 
waves, water waves nevertheless can be used as an analogy to illustrate 
the principles involved. 

If a long strip of flexible material were allowed to float on 
the surface of a body of water, it would follow the wave action. If, 
for purposes of mathematical analysis, this flexible strip were rep- 
resented by a number of short floats of rigid construction but con- 
nected by flexible links, the mathematical model would appear as in 
figure (7). This representation would be an approximation, but it 
would involve negligible error because the small floats would ride 
the waves almost exactly like the flexible strip. If, however, the 
rigid floats are made comparatively long and they approach or exceed 
the length of the water wave, then the model would be as shown by 
figure (8). Obviously these long floats cannot follow the wave form 
closely. 
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Applying this analogy to a pile being hit by a ram, the 
importance of dividing the pile into unit lengths considerably shorter 
than the stress wave produced by a hammer blow can be appreciated. 
Fortunately, a pile driving impact usually produces a fairly long 
wave form, and unit lengths of the order of 5 to 10 feet will produce 
acceptable accuracy. In the special case where the exact form of 
the impact wave is being investigated, a smaller unit length of 1 or 
2 feet may be advisable. Under these conditions, it may also be ad- 
visable to divide the ram into a series of weights and springs as is 
done with the pile. For normal pile driving analysis, such refinements 
are not considered justified* 

There is a very important relationship between the unit length 
of pile selected and the time interval chosen for the calculations, 
the time interval to be used in the calculations is a function of the 
weight of the unit length selected and of the elasticity and density 
of the pile material. This relationship is derived and presented 
later, and governs the upper bound of the time interval to be used. 

On the other hand, choice of too small a time interval results in 
excessive computation for the accuracy required. 

The spring constants are determined by the formula given for 
K under "Capblock" above. 

The Ground — as indicated by figure (ii) provision, is made for 
ground resistance to be applied on each section of the pile. In 
order to make a calculation using the wave equation, an assumption 
must be made as to the amount of resistance which will be offered by 
the ground at each of the pile sections. Furthermore, a stress-strain 







. . - i~ •• . «■%** # «■» i- v4" 4 ■■ lw| ! UmM 










25 



relationship of the soil and viscous damping must be considered. For 
the purpose of this section, it will be assumed that these factors 
can be evaluated so that the method of bheir inclusion in the analysis 
may be developed. 

Resistance at the pile point— according to Chellis (ref. 5) 
the ground compresses elastically for a certain distance and then 
fails plastica l ly with constant ultimate resistance R u . This concept 
is illustrated by figure (9). The elastic deformation will be called 
n quake," and is denoted by the letter Q.' The plastic deformation will 
be referred to as "permanent set," and is identified as 's'; in figure 
(9) it is the distance AB or OC, There is a factor of importance 
which has not yet been included, and this is the rate of stress 
application. It is clear that the ground will offer more instanta- 
neous resistance to a rapid advancement of the pile tip than to a slow 
advancement. Viscous damping, therefore, will be included in the 
analysis to provide for the effect of rate of penetration. 

As has been shown, the wave equation calculation will give the 
instantaneous velocity of the point of the pile in any time interval. 

If J is designated as a damping constant and v^ the velocity of 
the pile point, then their product, JVp, can be used to increase or 
decrease ground resistance at the point of the pile so as to produce 
damping. The instantaneous damping resistance at point x, for example, 

in figure (9) is JvX . This value would be used for R^o ^ or a P ar " 

a ^ 

ticular time interval applied as shown by figure (U). This damping 
resistance is a resistance to driving the pile and does not, of course. 
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contribute to the bearing capacity of the pile in supporting static 
loads . 

Resistance Along the Sides of the Pile — side resistance is 
calculated in the same manner as at the pile point, except that a 
different damping constant J* is used, and it applies to R^ through 
R^ in figure (U). As the pile is driven, the soil under the pile 
point is displaced and caused to flow aside very rapidly. On the 
other hand, the soil at the sides of the pile, as it is being dirven, 
is not displaced nearly as much. It may be assumed, then, that J 1 
should be less than J. 

Distribution of Ground Resistance Between the Pile Sides and 
Point — the wave equation computation will utilize any distribution 
chosen. The problem of what distribution to choose will be discussed 
at length later. 

Figure (10) shows the model representation of the ground 
resistance alorg the sides of the pile. It must be remembered that 
ground quake is an elastic,, or recoverable, deformation and that the 
set associated with it is of a temporary nature. The deflection of 
the spring must, therefore, be limited to a maximum value less than 
or equal to Q, the ground quake. As the maximum allowed deflection 
of spring is reached, provision is made for plastic deformation 
of the ground, D^, in the normal downward direction, or in an upward 
direction in the case of pile rebound. Figure (ll) illustrates the 
way in which is handled in the calculations. It may be noted that 
does not take on a value until D m exceeds Q. From that point on, 

D^ cannot be less than D m - Q or more than D m + Q. 
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The value of the ground resistance on a particular pile sec- 
tion is computed by the equation 

\ - < D m ’ D mK + («) 

The point of the pile is dealt with in a similar, but slightly 
different, manner. Only downward plastic ground movement Dp is con- 
sidered. The maximum value of D^ is the permanent set of the pile 
which is designated as s . Figure (12) illustrates the way in 
which the calculations are handled. The value of the ground re- 
sistance at the point is computed by the equation 

Rp - (Dp - Dp (1 + Jv p ) (16) 

C. Illustrative Problem 

In order to bring the theory more sharply into focus, a 
hypothetical problem will be set up and solved manually using the 
general format described by E. A. L. Smith (ref. 27). 

Let figure (13) represent a steel pile 30 feet long. The 
weight of the ram is 5,000 pounds, and its velocity at impact as 
determined by equation (10) is found to be 10 feet per second. The 
capblock is hardwood, has a coefficient of restitution of 0.5, and 
its spring constant as determined by equation (5) is 2,000,000 pounds 
per inch. The pile cap weighs 500 pounds, aid is not able to transmit 
tension. No cushion block or head packing will be used, so K 2 will 
be the elasticity of the first pile section whose coefficient of 
restitution will be taken as 1.00. The pile will be divided into 
10-foot unit lengths, each of which weighs 500 pounds. The spring 
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constants K^, and Kj^ for the pile as computed by equation (6) 
are each taken to be U, 000, 000 pounds per inch. Ground quake will be 
0.10 inch. The ultimate ground resistance has been determined to be 
220,000 pounds. It is distributed as follows: 10,000 pounds on each 

of the two top sections of the pile and 200,000 pounds on the bottom 
10-foot section of pile. Dividing the ultimate ground resistance by 
ground quake will give the soil spring constants as follows: 

® = 100,000 pounds per inch 

and Kji = 2,000,000 pounds per inch 

The time interval for the calculations will be 1/^000 seconds. The 
dairping constants J and J' will both be taken as equal to 0.1 for 
purposes of this example, although, as noted before, J* may be smaller 
than J. 

The calculations may then be performed by slide rule and 
recorded in tabular form as shown on the following pages. Although 
space does not permit here, it is recommended that when making such a 
table a single large sheet be used. In referring to Table A, on each 
line values are first computed for D^, D^, D y D^, D^j then C^, C^, 

C 3, C U> C^j etc., until all computations have been completed. The 
next line is approached in the same order, and so on. 

It may be of'interest to note from Table A that reaches a 

first maximum at time interval 8; restitution takes place from time 

interval 9 through 22. Recompression occurs from time interval 23 

through 2? at which time a second C n is reached. Restitution 

x max 

continues until the end of the calculations. 
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TABLE A 



C;L = ( 2 x 10 6 )( C 1 ), z 



Time 

Interval 


d -^+ 0.003 vx 


D 1 - d 2 


except during 
restitution 


~ F 1 


1 

v l + 621,000 


0 


0 


0 


0 


0 


10.00 


1 . 


0.03 


0.03 


60,000 


- 60,000 


9.9035 


2 . 


0.06 


0.057 


114,000 


- 114,000 


9.7202 


3 . 


0.0892 


0.0784 


156,800 


- 156,800 


9.468 


4 . 


0.1176 


0.0934 


186, 800 


- 186,800 


9.167 


5 . 


0 . 1U51 


0.1025 


205,000 


- 205,000 


8.837 


6 . 


0.1716 


0.1073 


214,600 


- 214,600 


8.491 


7 . 


0.1971 


0.1092 


218,400 


- 218,400 


8.139 


8 . 


0.2215 


0.1096 


219,200 


- 219,200 


7.786 


9 . 


0.2449 


0.1091 


215,200 


- 215,200 


7.440 


10 . 


0.2672 


0.1081 


207,200 


- 207,200 


7.106 


11 . 


0.2885 


0.1069 


197,600 


- 197,600 


6.788 


12 . 


0.3088 


0.1058 


188,800 


- 188,800 


6.484 


13 . 


0.3283 


0.1053 


184,800 


- 184,800 


6.186 


14 . 


0.3469 


0.1051 


183,200 


- 183,200 


5.891 


15 . 


0.3646 


0.1049 


181,600 


- 181,600 


5.599 


16 . 


0 . 381U 


0.1043 


176,800 


- 176, 800 


5.314 


17 . 


0.3973 


0.10334 


169,120 


- 169,120 


5.042 


18 . 


0.4124 


0.10224 


160,320 


- 160,320 


4.784 


19 . 


0.4267 


0.10164 


155,520 


- 155,520 


4.534 


20 . 


0.4403 


0.10238 


161,440 


- 161,440 


4.274 


21 . 


0.45312 


0.10479 


180, 720 


- 180.720 


3.984 


22 . 


0.46507 


0.10855 


210, 800 


- 210,800 


3.644 


23 . 


0.47600 


0.11266 


225,320 


- 225,320 


3.281 


2 k . 


0.48584 


0.11687 


233,740 


- 233,740 


2.905 


25 . 


0.49456 


0.12095 


241,900 


- 241,900 


2.516 


26 . 


0.50211 


0.1243 


248,600 


- 248,600 


2.116 


27 . 


0.50846 


0.12612 


252,240 


- 252,240 


1.710 


28 . 


0.51376 


0.12580 


249,680 


- 249,680 


1.309 


29 . 


0.51769 


0.12294 


226,800 


- 226, 800 


0.944 


30 . 


0.52052 


0.11925 


197,280 


- 197,280 


0.626 


31 . 


0.52240 


0.11684 


178,000 


- 178,000 


0.340 


32 . 


0.52342 


0.11697 


179,040 


- 179,040 


0.052 


33 . 


0.52357 


0.11913 


196,320 


- 196,320 


- 0.264 


3 k . 


0.52278 


0.12142 


214,640 


- 214,640 


- 0.610 


35 . 


0.52095 


0.12185 


218,080 


- 218,080 


- 0.962 


36 . 


0.51806 


0.11971 


200,960 


- 200,960 


- 1.284 


37 . 


0.51421 


0.11586 


170,160 


- 170,160 


- 1.558 


38 . 


0.50954 


0.11247 


143,040 


- 143,040 


- 1.789 


39 . 


0.50417 


0.11129 


133,600 


- 133,600 


- 2.004 
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TABLE A (cont'd. ) 




.003 v 2 




F = 

*2 

l*xl0 6 C„ 




V, 



v + 



2 62,100 



0 

0.003 
0.01081 
0.021*23 
.01*257 
. 06U3U 
.08788 
.11190 
.13575 
.15905 
.18157 
.20297 
0.22302 
. 21*183 
0.25972 
0.27706 
0.29396 
0.31016 
0.32506 
0.33792 
0.31*833 
0.35652 
0.3633U 
0.36897 
0.37361 
0.37781 
0.38231* 
0.38796 
0.391*75 
0 . 1*0127 
0.1*0556 
0.1*061*5 
0 . 1 * 01 * 1 * 1 * 
0.1*0136 
0.39910 
0.39835 
0.39835 
0.39707 
0.39288 



0 

0.003 
0.0102 
0.02123 
0.03357 
0.01*1*5 
0.05208 
0.0557 
0.0 5655 
0.05595 
0.05517 
0.051*17 
0.05262 
0.05053 

0.01*825 

0.01*65 

0.01*59 

0.01*68 

0.01*91*8 

0.05303 

0.05667 

0.05980 

0 . 0621 * 1 * 

0.06359 

0.06275 

0.0601*6 

0.0571*1 

0.05635 

0.05808 

0.06092 

0.06208 

0.05977 

0.051*66 

0.01*91*1 

0.01*671 

0 . 01*636 

0.01*917 

0.0508 

0.01*769 



0 

12,000 
1*0,800 
81*, 920 
13 k , 280 

178,000 
208,320 
222,800 

226,200 
223,800 

220,680 

216,680 
210,1*80 
202,120 

193,000 

186,000 
183,600 

187,200 
197,920 
212,120 

226,680 

239.200 
21*9,760 
251*, 360 

251, 000 
21*1, 81*0 
229,61*0 
225,1*00 
232,320 
21*3,680 
21*8,320 
239,080 
218,61*0 
197,61*0 
186, 81*0 
185,1*1*0 

196,680 

203.200 
190, 760 



60', 000 

102,000 

116,000 
101,880 
70,720 
36,600 
10,080 
-3,600 
-11,000 
-16,600 
-23,080 
-27,880 
-25,680 
-18,920 
-11,1*00 
- 9,200 
-ll*,l*80 
-26, 880 
-1*2,1*00 
-50,680 
-1*5,960 
-28,1*00 
-21*, 1*1*0 
-20,620 
- 9,100 
+ 6,760 
+22, 600 
+21*, 280 
- 5,520 
-1*6,1*00 
-70,320 
-60, 01*0 
-22,320 
+17, 000 
+31,21*0 
+15,520 
-26,520 
-6 0,160 
-57,160 



0.965 
2.605 
1*.1*75 
6.115 
7.255 
7.81*5 
8.007 
7.91*9 
7.772 
7.505 
7.131* 
6.681* 
6.270 
5.965 
5.781 
5.633 
5.1*00 
l*. 967 
U.285 
3.1*69 
2.729 
2.272 
1.878 
1.51*6 
1 . 1*00 
1.509 
1.873 
2.261* 
2.175 
1.1*29 
0.298 
-0.669 
-1.028 
-0.751* 
-0.251 
- 0.001 
-0.1*28 
-1.398 
-2.320 
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38 
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TABLE A (cont'd. ) 



+0.003 




V 

lpcl0 6 C_ 




0 

0 

0.0006 

0.0030 

0 . 0090 : 

0.0198 

0.0358 

0.0562 

0.0792 

0.1031 

0 . 1261 ; 

0.1U88 

0 . 1701 ; 

0.1913 

0.2111;7 

0.23056 

0.2U800 

0.26332 

0.27558 

0.28U89 

0.29166 

0.29672 

0.30090 

0.30538 

0.31086 

0.31735 

0.32U93 

0.33161 

0.33667 

0.3U035 

0.3U3U8 

0.3U668 

0.3U978 

0.35195 

0.35239 

0.35200 

0.3U918 

0.3U627 

0.3U519 



0 

0 

0.0006 

0.0029 

0.00823 

0.01693 

0.02809 

0.0396 

0.0U883 

0. 051+2I; 

o.o551i7 

0.05U02 

0.05175 

0.05003 

O.OU965 

0.05095 

0.0539 

0.05791 

0.06122 

0.06289 

0.062U9 

0.06151 

0.05803 

0.05557 

0.05U55 

0.05181 

0.05893 

0 . 06166 

0,06228 

0.06096 

0.05898 

0.05751 

0.05670 

0.05571 

0.05363 

0.05139 
o.ohlhh 
0.0UU19 
0 . 0 hb66 



0 

0 

2,1;00 

11,600 
32,920 
67,720 
112,360 
158, UOO 
195,320 
216,960 
221, 880 
216,080 
207,000 
200,120 

198,600 
203,800 

215,600 
231,61*0 
21;!;, 880 
251,560 
2U9,960 
2U6,0U0 
232,120 
222,280 
218,200 
207, 2U0 
235,720 
2l;6,61;0 
2U9,120 
2U3,8Uo 
235,920 

230, 0U0 
226, 800 
222, 8i;0 
21i;,520 
205,560 

189.760 

156. 760 

178, 61;0 



0.0031 
0.0261; 
0. OI4.88 
0.070 H 
0.0913 
0.1111*7 
0.13056 
0.11*800 
0.16332 
0.17558 
0.181*89 
0.19166 
0.19672 
0.20090 
0.20538 
0.21086 
0.21735 
0.221*93 
0.23161 
0.23667 
0.21*035 
0.21*31*8 
0.21*668 
0.21*978 
0.25195 
0.25299 
0.25200 
0.21*918 
0.21*627 
0.21*519 
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TABLE A (cont'd.) 



Time 

Inter- 

val 




D 3 ‘ D 3 


K'(D 3 -Dp 


1 + J'v^ 


nr 3 3 

(1 + J'v 3 ) 


Z 3 " 

F 2 - F 3 - R 3 


3 

+ Z 3 


v 3 62.100 


1 










0 


0 


0 


2 


0 


0 


0 


0 


0 


12,000 


0.1931 


3 


0 


0.0006 


60 


1.0193 


60 


38,31*0 


0.8091 


l* 


0 


.0030 


300 


1.0809 


321* 


72,996 


1.985 


5 


0 


0.0090 


900 


1.1985 


1,079 


100, 281 


3.601 


6 


0 


0.01980 


1,980 


1.3601 


2,690 


107,590 


5.336 


7 


0 


0.03580 


3,580 


1.5336 


5,1*90 


90,1*70 


6.795 


8 


0 


0.0562 


5,620 


1.6795 


9,1*30 


51*, 970 


7.680 


9 


0 


0.0792 


7,920 


1.7680 


11*, 000 


16,880 


7.952 


10 


0.0031 


0.10 


10,000 


1.7952 


17,952 


-11,112 


7.773 


11 


0.0261* 


0.10 


10,000 


1.773 


17,730 


-18,930 


7.1*58 


12 


0.01*88 


0.10 


10,000 


1.71*58 


17,1*58 


-16,858 


7.186 


13 


O'. 0701* 


0.10 


10, 000 


1.7186 


17,186 


-13,706 


6.965 


11* 


0.0913 


0.10 


10,000 


1.6965 


16,965 


-11*, 965 


6.721* 


15 


0.1111*7 


0.10 


10,000 


1.6721 


16,721* 


-22,321* 


6.361* 


16 


0.130^6 


0.10 


10,000 


1.6361* 


16,361* 


-31*, 161* 


5.811* 


17 


0.11*800 


0.10 


10, 000 


1.581U 


15,811* 


-1*7,811* 


5.01*1* 


18 


0.16332 


0.10 


10,000 


1.50UU 


15,01*1* 


-59, 1*81* 


U.087 


19 


0.17558 


0.10 


10,000 


1.1*087 


11*, 087 


-61,01*7 


3.103 


20 


0.181*89 


0.10 


10,000 


1.3103 


13,103 


-52,51*3 


2.257 


21 


0.19166 


0.10 


10,000 


1.2257 


12,257 


-35,537 


1.686 


22 


0.19672 


0.10 


10,000 


1.1686 


11,686 


-18,526 


1.391* 


23 


0.20090 


0.10 


10,000 


1.1391* 


11,391* 


+ 6,21*6 


1.1*91* 


2i* 


0.20538 


0.10 


10,000 


1.1U9U 


11,1*91* 


+20,586 


1.826 


25 


0.21086 


0.10 


10,000 


1.1826 


11,826 


+20,971* 


2.163 


26 


0.21735 


0.10 


10,000 


1.2163 


12,163 


+22,1*37 


2.525 


27 


0.221*93 


0.10 


10,000 


1.2525 


12,525 


-18,605 


2.225 


28 


0.23161 


0.10 


10,000 


1.2225 


12,225 


-33,1*65 


1.686 


29 


0.23667 


0.10 


10,000 


1.1686 


11,686 


-28,1*86 


1.227 


30 


0.21*035 


0.10 


10,000 


1.1227 


11,227 


-11,387 


1.01*1* 


31 


0.21*31*8 


0.10 


10, 000 


l.iol*!* 


11,01*1* 


+ 1,356 


1.066 


32 


0.21*668 


0.10 


10,000 


1.1066 


11,066 


- 2,026 


1.033 


33 


0.21*978 


0.10 


10,000 


1.1033 


11,033 


-19,193 


0.721* 


3 b 


0.25195 


0.10 


10, 000 


1.0721* 


10,721* 


-35,921* 


0.11*6 


35 


0.25239 


0.10 


10,000 


1.011*6 


10,11*6 


-37,826 


-0.1*61* 


36 


0.25239 


0.09961 


9,961 


0.9536 


9,1*90 


-29,610 


-0.91*1 


37 


0.25239 


0.09679 


9,679 


0.9059 


8,760 


- 1,81*0 


-0.971 


38 


0.25239 


0.09388 


9,388 


0.9029 


8, 1*1*0 


+38,000 


-0.359 


39 


0.25239 


0.09280 


9,280 


0.961*1 


8,91*2 


+ 3,178 


-0.308 
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TABLE A (cont ‘d. ) 



D i* 



dj^ +0.003 




F a * 

ipcio 6 ^ 



- Q 



0 0 0 
0.00012 0.00012 1*80 

0.00077 0.00075 3,000 

0.00287 0.00269 10, 760 



0.00771 


0. 


.01660 


0. 


.03037 


0. 


0.0U886 


0. 


o. 07093 


0. 


0.09U79 


0. 


0.11865 


0. 


0.1U127 


0. 


0.16182 


0. 


0.17961 


0. 


o.i9iao 


0. 


0.205U1 


0. 


0.211;36 


0. 


.22200 


0. 


0.22917 


0. 


0.23521 


0. 


0.2l;287 


0. 


0.21*981 


0. 


0.25631 


0. 


0.2619U 


0. 


0.26600 


0. 


0.26995 


0. 


0.27U39 


0. 


0.27939 


0. 


0.281*50 


0. 


0.28919 


0. 


0.29308 


0. 


0.29621; 


0. 


0.2987 6 


0. 


0.30061 


0. 


0.3017U 


0. 


0.30208 


0. 


0.30053 


0. 



0069 27, 

01382 55 , 

02329 93, 

03376 135, 

01*313 172, 

01*978 199, 

05365 212*, 

05636 225, 

05975 239, 

061*02 256, 

06663 266, 

06668 266, 

061*77 259, 

06219 21*8, 

06008 21*0, 

05518 220, 

05865 23U, 

051*75 219, 

05603 221;, 

05696 227, 

05670 226, 

05627 225, 

05651* 226, 

05 71*5 229, 

0583U 233, 

05851* 23U, 

05777 231, 

05626 225, 

051*1*0 217, 

052U3 20 9, 

01*898 195, 

01*61*1 185 , 

OU358 17U, 



600 


- 


280 


- 


160 


— 


01*0 


am 


520 


mm 


120 


mm 


600 


0.01865 


hho 


O.OU127 


000 


0.06182 


080 


0.07961 


520 


o.o9iao 


720 


0.1051a 


080 


0. 111*36 


760 


0.12200 


320 


0.12917 


720 


0.13521 


600 


0.11*287 


000 


0.11*981 


120 


0.15631 


81*0 


0.1619U 


800 


0.16600 


080 


0.16995 


160 


0.171*39 


800 


0.17939 


360 


0.181*50 


160 


0.18917 


080 


0.19308 


ouo 


0.19621* 


600 


0.19876 


720 


0.20061 


920 


0.20171* 


61*0 


0.20208 


320 


0.20053 
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TABLE A (cont'd.) 



Time 

Inter- 

val 










H), - 

K u^h- D u>- 

1 + J'v^ (1 + j'v^) 



1* 

F 3" F U" R t 



z i* 



1 

2 



3 












2,1*00 


0.03865 


1* 


0 


0.00012 


12 


1.003865 


12.05 


11,108 


0.21765 


5 


0 


0.00077 


77 


1.0218 


78.75 


29,81*1 


0.699 


6 


0 


0.00287 


287 


1.0699 


307 


56,653 


1.613 


7 


0 


0.00771 


771 


1.1613 


895 


83,865 


2.963 


8 


0 


0.01660 


1,660 


1.2963 


2,11*5 


100,975 


1*.590 


9 


0 


0.03037 


3,037 


1.1*59 


1*,1*20 


97,71*0 


6.163 


10 


0 


0.01*886 


I*, 886 


1.6163 


7,900 


71*, 020 


7.356 


11 


0 


0.07093 


7,093 


1.7356 


12,300 


37,060 


7.953 


12 


0 


0.091*79 


9,1*79 


1.7953 


17,000 


-1*0 


7.952 


13 


0.01865 


0.10 


10,000 


1.7952 


17,952 


-25,550 


7.51*1 


1U 


0.01*127 


0.10 


10,000 


1.751*1 


17,51*1 


-1*2,861 


6.851 


15 


0.06182 


0.10 


10,000 


1.6851 


16,851 


-57,251 


5.929 


16 


0.07961 


0.10 


10,000 


1.5929 


15,929 


-68,209 


U.830 


17 


0.091*10 


0.10 


10,000 


1.1*830 


11*, 830 


-65,750 


3.770 


18 


0.105m 


0.10 


10,000 


1.3770 


13,770 


-1*8,850 


2.981* 


19 


0.111*36 


0.10 


10, 000 


1.2981* 


12,981* 


-27,181* 


2.51*7 


20 


0.12200 


0.10 


10,000 


1.251*7 


12,51*7 


- 9,71*7 


2.390 


21 


0.12917 


0.10 


10, 000 


1.2390 


12,390 


- 2,750 


2.31*6 


22 


0.13521 


0.10 


10,000 


1.231*6 


12,31*6 


+12,971* 


2.555 


23 


0.11*287 


0.10 


10,000 


1.2555 


12,555 


-15,035 


2.313 


21* 


0.11*981 


0.10 


10, 000 


1.2313 


12,313 


- 9,033 


2.1678 


25 


0.15631 


0.10 


10,000 


1.21678 


12,168 


-18,088 


1.8768 


26 


O.1619I* 


0.10 


10,000 


1.18768 


11, 877 


-32,1*77 


1.351*8 


27 


0.16600 


0.10 


10, 000 


1.1351*8 


11,355 


- 2,1*35 


1.3156 


28 


0.16995 


0.10 


10,000 


1.13156 


11,316 


+10,21*1* 


1.1*806 


29 


0.171*39 


0.10 


10,000 


1.11*806 


11,1*81 


+11,1*79 


1.6656 


30 


0.17939 


0.10 


10,000 


1.16656 


11,666 


+ 2,371* 


1.7038 


31 


0.181*50 


0.10 


10,000 


1.17038 


11,701* 


- 9,11*1* 


1.5568 


32 


0.18917 


0.10 


10,000 


1.15568 


11,557 


-15,677 


1.301*8 


33 


0.19308 


0.10 


10,000 


1.1301*8 


11,305 


-15,585 


1.0538 


3 k 


0.19621* 


0.10 


10, 000 


1.10538 


11,051* 


-13,251* 


0.81*03 


35 


0.19876 


0.10 


10,000 


1.081*03 


10,81*0 


-13,920 


0.6163 


36 


0.20061 


0.10 


10, 000 


1.0616 


10,616 


-li*, 776 


0.378 


37 


0.20171* 


0.10 


10, 000 


1.0378 


10,378 


-16,538 


0.112 


38 


0.20208 


0.10 


10,000 


1.0112 


10,112 


-38,992 


- .515 


39 


0.20208 


0.0981*5 


9,81*5 


0.91*85 


9,31*0 


- 5,020 


-0.596 



29 f 



d 5 = 

'ime + 

" S r ' °- 003 v 5 b 5- <3 



TABLE A (cont'd. ) 



D' 

P 



Rj - 

(K5) 1*J v ; (KjKDj-Dp)- 

%-Dp (%-Dp) (1 +Jtj ) 




VS,. SzTIT 



1 

2 

3 



b 


0 


0 


5 


0.00002 


0 


6 


0.00019 


0 


7 


.00086 


0 


8 


.00278 


0 


9 


.00708 


0 


10 


.01510 


0 


11 


.02780 


0 


12 


.0b50i 


0 


13 


.06500 


0 


lb 


.08b 91 


0 


15 


0.10207 


0.0207 0.0207 


16 


0.11559 


0.01559 0.01559 


17 


0.127b7 


0.027b7 0.027b7 


18 


0.13873 


0.03873 0.03873 


19 


.lb959 


0.0b9 59 O.Ob959 


20 


.15981 


0.05981 0.05981 


21 


0.16909 


0.06909 0.06909 


22 


0.17733 


0.07733 0.07733 


23 


.I8b22 


0.08b22 0.08b22 


2b 


0.19506 


0.09506 0.09506 


25 


0.20028 


0.10028 0.10028 


26 


0.20b98 


0.10b98 0.10b98 


27 


0.20931 


0.10931 0.10931 


28 


0.21368 


0.11368 0.11368 


29 


0.21785 


0.11785 0.11785 


30 


0.2219b 


0.1219b 0.1219b 


31 


0.22616 


0.12616 0.12616 


32 


0.23063 


0.13063 0.13063 


33 


0.23531 


0.13531 0.13531 


3b 


0.23998 


0.13998 0.13998 


35 


0.2bb36 


0.1bb36 0.1bl36 


36 


0.2b8l8 


0.1b3l8 0.1b 8l 8 


37 


0.25276 


0.15276 0.15276 


38 


0.25567 


0.15567 0.15567 


39 


0.25695 


0.15695 0.15695 



0 



0.00002 


bo 


1.0008 


bo. 03 


0.00019 


380 


1.005bb 


382 


0.00086 


1,720 


1.022b 


1,693 


0.00278 


5,560 


1.0639 


5,920 


0.00708 


lb, 160 


l.lb3b 


16,190 


0.01510 


30,200 


1.267b 


38,250 


0.02780 


55,600 


l.b23b 


79,1°° 


0.0b5oi 


90, 020 


1.5737 


lbl, 700 


0.06500 


130,000 


1.6662 


216,200 


0.08b9l 


169,820 


1.6636 


282,300 


0.10 


200, 000 


1.5721 


31b, b20 


0.10 


200,000 


l.b508 


290,160 


0.10 


200, 000 


1.3959 


279, 18° 


0.10 


200,000 


1.3755 


275,100 


0.10 


200,000 


1.3620 


2 72, boo 


0.10 


200, 000 


1.3b06 


268,120 


0.10 


200,000 


1.309b 


261, 880 


0.10 


200, 000 


1.27b7 


25b, 9b0 


0.10 


200,000 


1.2296 


2b5,920 


0.10 


200, 000 


1.211b 


2b2, 280 


0.10 


200,000 


1.1739 


23b, 780 


0.10 


200,000 


1.1567 


231, 3b0 


0.10 


200,000 


1.1511 


230,220 


0.10 


200, 000 


l.lb56 


229,120 


0.10 


200,000 


1.1391 


227,820 


0.10 


200,000 


1.136b 


227,280 


0.10 


200,000 


l.lb05 


228,100 


0.10 


200,000 


l.lb90 


229,800 


0.10 


200,000 


1.1560 


231,200 


0.10 


200,000 


1.1558 


231,160 


0.10 


200,000 


l.lb6o 


229,200 


0.10 


200, 000 


1.1273 


225, b60 


0.10 


200, 000 


1.1527 


230, 5b0 


0.10 


200,000 


1.0969 


219,380 


0.10 


200,000 


1.0b25 


208,500 



U80 0.007' 
2,960 0.055i 
10,378 .222: 
25,907 .639! 

U9,360 l.U3ii 
76,970 2.67b 
96,790 b.23b 
93,b20 5.737 
57,U20 6.662 

- 1,600 6.636 
-56,860 5.721 
-75,U20 U.508 
-3b, 080 3.959 
*12,660 3.755 

- 8,380 3.62C 

-13,320 3.U06 
-19,360 3.09b 
-21,560 2.7b7 
-3b, 220 2.296 
-11,320 2 . llit 
-23,280 1.739 
-10,660 1.567 

- 3,500 1.511 

- 3,U20 1.U56 

- b,0b0 1.391 

- 1,660 1.36b 
+ 2,520 1.U05 
+ 5,260 1.U90 

+ U, 360 1.560 
-120 1.558 

- 6,120 I.U60 
-11,600 1.273 
-l5,7bO 1.527 
-3b, 620 0.969 
-33,7bO 0.b25 
-3b, 180 -0.125 
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It may also be noted that permanent ground displacement was 
not recorded until it exceeded the value assigned ground quake which 
was 0,10 in this problem. The time intervals at which these occurred 
were 

for , time interval 10 
, time interval 13 
D(t , time interval 15 

The calculations were terminated when all the velocities became nega- 
tive j effects after this are considered to be secondary and of little 
importance for the purposes of this study. 

A graphical presentation of the results is shown in figures 
(lit) through (16) inclusive. 

The reason that a time interval of 0.00025 seconds was chosen 
may not be clear. Generally, it may be said that the best time inter- 
val to use is one that produces a completely stable calculation. This 
means that when the results are plotted there should be no sharp 
wiggles or peaks. If such irregularities are present, either an error 
has been made in the calculations or else the time interval chosen is 
too large. Once the unit lengths of the pile have been chosen, the time 
interval must be selected to correspond. Choosing too small a time 
interval would result in an excessive amount of calculation with only 
a very limited increase in accuracy. 

Each spring in the figure (13) has a critical time interval 
which is the time that is required far a stress wave to traverse the 
spring and its associated weight. Stress waves travel in both directions 
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at a speed of j 3 where P is the mass per unit of volume. The 



following equations may be derived for the critical time interval for 
spring Kj^ which is called T m . 

For wave motion downward: 



The smallest value of T m obtained by using these formulas is the 
critical one. About one half of the critical value should be selected 
as the time interval to be used in any particular problem so as to 
prevent instability from arising due to other factors not included in 
the above formulas such as quake, damping, and coefficients of restitu- 
tion. When dividing the pile into ten-foot sections for purposes of 
mathematical representation, the recommendations of Smith (ref. 28) 
have been found to be generally satisfactory: 



During some of the test cases min, instability resulted in using this 
time interval for wood piles, and recomputations were made using smaller 
time intervals. Instability can generally be detected by plotting 
the results as noted previously or by comparing the velocity of the 
pile cap and the pile tip with the velocity of the ram at impact. If 
either of these velocities exceeds twice the ram velocity at impact, in- 
stability is probably present. The velocity check should not be relied 
upon alone, however, to establish the stability of the calculations. 




For wave motion upward: 




For steel piles — 0.00025 seconds 

For wood piles — 0.00025 seconds 



For concrete piles — 0.00033 seconds 
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IV. COMPUTER ADAPTATION 



A. General 

As shown by the illustrative example, pile driving may be 
analyzed by the wave equation using a numerical solution method and a 
slide rule. 'This approach is a valid one and will produce satis- 
factory results. The main problem, of course, is the amount of time 
and work required to make the calculations. The example problem was 
selected as being about as simple as possible, with convenient figures 
to deal with numerically, and still illustrate all aspects of the cal- 
culations. Nevertheless, it required about twelve hours of manual 
computation to get the final set. It should also be noted that an 
error made at the beginning of a solution will carry throughout the 
calculations. The usual problem encountered in practice would probably 
take several man-days to complete. 

Fortunately, there is a quicker and better way of making these 
calculations. Use may be made of high-speed digital computers which 
can solve up to 800 such practical problems in about 5 minutes, de- 
pending upon the computer used. As a consequence, use of the wave 
equation for pile driving analysis has become a practical matter. 

Although it is not the purpose here to go into the details of 
programming, it may be well to mention some of the broader aspects for 
those totally unfamiliar with computers. For a computer to be effective 
in solving engineering- type problems, the problem-solving procedures 
must be presented to the computer in a language which it can "understand." 
The computer's basic language consists of elementary instructions such 
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as add or subtract. The problem-solving procedure must be translated, 
then, into instructions which the computer can obey. This translation 
is called programming and can be carried out entirely by a person, or 
the computer can assist in the process by use of a compiler. 

A compiler is a large set of computer instructions which can 
accept a problem-solving procedure, written in a form resembling the 
procedure, and produce from it the proper elementary machine instruc- 
tions that will solve the problem. The algebraic compiler system used 
for the programs here is called Fortran, which stands for formula 
translation. A procedure to be followed in solving a problem with 
Fortran is specified by a series of statements, of which there are 
several types. One type specifies arithmetic operations: another calls 
for reading a data card or printing results; a third specifies the 
sequence in which the statements are to be executed; and a fourth type 
is made up of statements which provide information about .the procedure 
without themselves causing any action. All of these statements taken 
together form a source program. The programs included in the appendices 
of this writing are source programs. When punched on cards and put 
through the machine, it is translated by the Fortran system into simple 
sets of machine instructions which is called an object program. For 
those who wish to pursue the subject further, reference (21) provides 
in the space of 83 pages all the information necessary for a person to 

learn to develop programs of his own. 

B. Description of Programs Developed 

1. Basic Programs 

A basic program in Fortran language for the IBM 7090 computer 
was prepared, and it follows the method of calculation illustrated in 



3U 



the example problem. This program requires that input data be supplied 
to indicate the number and size of the weights in the mathematical 
model, the time increment to be used, the spring constants of the pile, 
the spring coefficients of the soil, the coefficients of restitution 
of the capblock and cushion block (if used), the velocity of the ram at 
the instant of impact, the value of ground quake, values of point and 
side damping, whether or not the pile cap can transmit tension, and 
whether or not side resistance is present. Output of the basic program 
includes the input data for purposes of identification and checking, 
and then for each time increment the force and resistance applied to 
each block of the model together with displacement and velocity of each 
block at that particular instant. The permanent set of the tip of the 
pile is also recorded. The program is arranged to stop the computer 
after 300 time- intervals or before if 

(a) the velocity of the pile cap exceeds twice the velocity the 
ram had at the instant of impact; 

(b) the velocity of the pile tip exceeds twice the velocity the 
ram had at the instant of impact; 

(c) the velocities of all the blocks are simultaneously negative 
or equal to zero. Velocities toward the top of the pile are taken to 
be negative. 

The basic program can be used to study stresses, velocities, 
accelerations, and permanent set of the pile as they vary with time. 

The program may be run with one or more sets of data at a time. This 
program in its entirety, together with instructions on how to use it, 
is included as appendix A. 
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Having developed this program, the question was how best to 
use it in trying to detemine if correlations could be achieved with 
available field data on driving records, load tests, and borings. Of 
the variables that must be provided as input data, the greatest un- 
certainty was with the values to assign ground quake, point and side 
damping, and distribution of resistance along the sides of the pile. 

2. Vary RU Program . 

Correlation was approached in two different ways. The first 
method was particularly suited to the very few cases in which careful 
and complete field tests were reported and analyzed comprehensively. 
In these cases, reasonably accurate estimates of the variables except 
ground quake and damping factors could be made. Another program was 
developed for this approach which was basically the original program, 
but automatically for the one set of input data varies ground quake, 
ultimate ground resistance, and damping factors over a wide range of 
values. The output includes far each variation of basic information 
indicated the ultimate ground resistance, the blows per inch at final 
set, maximum tension and compression at the head, mid-length, and tip 
of the pile, the time interval the calculations stopped and why, and 
the original input data. -This program in its entirety is included 
with instructions on how to use it as appendix B. 

3. Researcher Program . 

The second method of approaching correlation was more of a 
statistical approach, and it was designed to try and glean as much 
information as possible from published field data which were not as 
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comprehensively reported as data used in the first method of approach. 
References (15) and (16) were used extensively in this attempt. Certain 
pile sections, lengths, and driving hammers were chosen because more 
reasonably good data were available for them than for other types. 
Another computer program was developed around the original program and 
designed to systematically produce in quantity data from which curves 
of ultimate ground resistance vs. number of blows per inch could be 
plotted for various values of damping constants and ground quake and 
side resistance. This program, together with instructions for its use, 
is included as appendix C, It is designed to be run with one or more 
sets of data. Its output includes for each of the above automatic 
variations of input, the ultimate ground resistance, values of soil 
spring coefficients, the permanent set, the blows per inch at the final 
set, maximum tension and compression at the head, mid-length, and tip 
of the pile, the time interval the calculations stopped and why, and 
the original input data. Thus about eighty complete problem solutions 
are obtained for each set of data run. Ordinarily three sets of data 
can be run at a time in this manner and take no more than five minutes' 
IBM 7090 computer time. Manual solution would require at least one man- 
year for just one set of data. Running time on a 709 computer would 
be 20 to 2$ minutes, but only about a minute if the 709U is used. 

U. Validation of Computer Programs 

Before using any computer program, great care must be taken to 
prove that it is correct and that all of its features behave as they 
are intended to. Making a manual solution is the only way to be cer- 
tain that the program is correctly written. To prove the original, or 
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basic, program, the data of the illustrative example presented 
herein was used as input, and the solution compared with the manual 
solution. A tabulation of the results of the two solutions for 
displacements of the ram in each time interval are found as Table B. 
It can be readily seen that the results compare as closely as one 
would expect considering that a slide rule was used in the manual 
computation. As a further check of the program, the illustrative 
problem presented by E. A. L. Smith (ref. 28) was computed using the 
basic program, and the results plotted. They compared in every 
detail with Smith's solution. 

As additional programs were developed around this basic 
program, they were tested with both of these illustrative examples 
and verified to insure that subsequent work would be based on sound 
programs . 

5 . Computer Program for the Hiley Dynamic Formula 

The Hiley type formula is thought by many to be one of the 
better dynamic formulas in general use today. Although not accepted 
by all, it is becoming more popular in the United States and has been 
included in some building codes. In order to facilitate comparison 
of results obtained by the wave equation and those obtained by using 
the Hiley formula, a computer program was developed for it also. 

This program requires that input data be provided to give the hammer 
data; pile material; the pattern, either rectangular or triangular, 
of side resistance distribution; coefficient of restitution; weight 
of the pile; area of the pile at the top, mid-length, and bottom; 
the ground quake; the modulus of elasticity of the pile material; 
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Table B 

COMPARISON OF COMPUTER SOLUTION WITH MANUAL SOLUTION 
FOR ILLUSTRATIVE EXAMPLE 



Time Interval 

1 

2 

3 

1 * 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Hi 

15 

16 

17 

18 

19 

20 
21 
22 
23 
21 * 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3U 

35 

36 

37 
33 
39 



Displacement 
Computer Results 



0.030 
0.060 
0.089 
0.117 
0.11i5 
0.171 
0.197 
0.221 
0.21*5 
0.267 
0.288 
0.309 
0.328 
0.31*7 
0.361* 
0.381 
0.397 
0.1|12 
0.1*26 
OJjl*0 
0.1*53 
0 .1*65 
0.1*76 
0.1*85 
o.l*9l* 
0.502 
0.508 
0.513 
0.517 
0.520 
0.522 
0.523 
0.523 
0.522 
0.522 
0.517 
0.513 
0.508 
0.502 



of Block 1 
Manual Solution 



0.030 
0.060 
0.089 
0.118 
0.11*5 
0.172 
0.197 
0.221 
0.21*5 
0.267 
0.288 
0.309 
0.328 
0.31*7 
0.365 
0.381 
0.397 
0.1*12 
0.1*26 
0. 1*1*0 
0.1*53 
0.1*65 
0.1*76 
0.1*85 
0.1*95 
0.502 
0.508 
o.5H* 
0.518 
0.520 
0.522 
0.523 
0.521* 
0.523 
0.521 
0.518 
o.5H* 
0.509 
0.50I* 



END OF CALCULATIONS 
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and the pile length in feet. ,_The program will automatically range 
over ultimate ground resistances from 20 to 380 tons in increments 
of 20 tons, writing out the set and blows per inch corresponding to 
each ground resistance as well as repeating the input data for each 
case. It computes for each value of the ground resistance the set 
associated with all the load carried by the point, one half of the 
load carried by the point, and none of the load carried by the point. 
Programming the Hiley formula may not be justified under usual cir- 
cumstances, but it was done in this case as a convenience. The 
program and instructions for its use are included as Appendix D. 

It should be noted that the Hiley formula is discriminatory to long, 
heavy piles and is subject to the same general limitations as any 
dynamic formula. 
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V. SOIL ENGINEERING CONSIDERATIONS 

As has been mentioned in the development of this method of 
pile driving analysis, values must be assigned to the dynamic soil 
properties of ground quake and damping. The distribution of forces 
acting on the pile tip and sides must also be evaluated to complete 
the information required for a solution. Furthermore, in order to 
approach correlation between dynamic driving resistance and static 
resistance, the manner in which frictional resistance along the 
sides of the pile acts in both the static and the dynamic cases must 
be explored. The subject of load testing must also be investigated 
insofar as the way in which it is done and the interpretation of 
results affect correlations. 

A. Side Friction Under Static Loading 

1. Granular Materials . 

Some recent work has been done on load transfer in the case 
of end bearing piles in which two instrumented H piles were tested 
(ref. 12). ^he piles were litf3P89 and 1UBP117 and were driven ap- 
proximately UO feet through granular soil varying from sandy silt 
to sand and gravel. The piles were end bearing on weathered silty 
shale. In spite of the fact that the piles were short and rela- 
tively stiff, they transferred as much as one third of their load 
to the surrounding soil by friction. Figure 17 shows the soil con- 
ditions and driving records. 

Upward shear along the pile develops as the pile shortens 
elastically under the applied static load at the top of the pile. 
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This movement of the pile causes distortion of the surrounding soil 
resulting in development of a frictional force resisting the dis- 
placement of the pile. Generally, it would be assumed that the 
greatest side resistance would be developed near the top of the pile 
and would diminish with depth. This action may be visualized by 
thinking of the pile as a spring embedded in the ground with a load 
applied at the top of the spring. One might also expect a nearly 
linear relationship between load applied and load transferred as long 
as no relative movement occurs between the soil and the pile. How- 
ever, a loading is reached at which the soil begins to slip along the 
pile. This zone of shear failure begins at the top of the pile and 
progresses downward. If the pile tip never moved, the entire applied 
load would be carried by friction along the pile sides, but even in 
the case of progressive downward failure, the failure zone would 
never extend down to the tip. In the actual case the pile tip does 
move downward as the applied load is increased so that it accepts a 
portion of it. 

The test results reported in reference 1 were presented from 
a viewpoint of over-all transfer of load to the soil by friction. 
Taking the data presented, and dividing each of the test piles into 
sections about $ feet long. Table C has been prepared in an effort 
to show the variation in friction force on each section of the pile 
as the test load was increased. In each of the test piles, section 
1 is the top section. A study of this table will show the validity 
of the description of the interaction between pile and soil, i.e., 



the frictional forces have a tendency to develop more at the top of 
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Table C 



FRICTION FORCE ON SECTION 
(IN TONS) 



Increment of Total 


0-75 Tons 


75-150 Tons 


150-225 Tons 


225-300 Tons 


Applied Load 


Change Total 


Change Total 


Change 


Total 


Change Total 






(at 75) 




(at 150) 




(at22^) 




(at 300 ) 


Section of IUBP 89 


















Pile and Soil 
Description 


















1— Sand & Gravel 


8 


8 


2 


10 


5 


15 


a 


19 


2— Sand k Gravel 


7 


7 


2 


9 


h 


13 


2 


15 


3 — Fine to Med. Sand 


6 


6 


1 


7 


-2 


5 


3 


8 


ii—Fine to Med. Sand 


0 


0 


3 


3 


a 


7 


1 


8 


5 — Sand k Gravel 


5 


5 


2 


7 


-3 


a 


a 


8 


6 — Sand k Gravel 


h 


h 


0 


h 


1 


5 


8 


13 


7— Sand k Gravel 
8~Sandy Silt, tip 


3 


3 


-1 


2 


5 


7 


5 


12 


in weathered 
shale 


1 


1 


5 


6 


0 


6 


a 


10 


TOTAL FRICTION FORCE 
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U8 




62 




93 


TOTAL FORCE ON HP 
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102 




163 




207 


Section of 1UBP117 


















Pile and Soil 
Description 


















1 — Sand k Gravel 


a 


h 


-1 


3 


-3 


0 


0 


0 


2 — Sand k Gravel 


a 


h 


1 


5 


-3 


2 


0 


2 


3— Fine to Med. Sand 


0 


0 


h 


a 


-a 


0 


2 


2 


U— Fine to Med. Sand 


2 


2 


3 


5 


-5 


0 


1 


1 


5— Sand k Gravel 


3 


3 


1 


h 


13 


17 


11.5 


28.5 


6 — Sand k Gravel 


3 


3 


1 


h 


a 


8 


- 1.5 


6.5 


7— Sand k Gravel 

8— Sandy Silt, tip 


6.5 


6.5 


3.5 


10 


-1 


9 


5 


ia 


in weathered 
shale 


0.5 


0.5 


1.5 


2 


a 


6 


-a 


2 


TOTAL FRICTION FORCE 




23 




37 




a2 




5 a 


TOTAL FORCE ON TIP 




52 




113 




183 




2a6 
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the pile. It also shows that local ikilure extends over a wide range 
of total applied load. For example, as the applied load is increased 
from 75 to 150 tons, the friction forces on section 7 of the Ik BP 89 is 
decreased by one ton; from 150 to 225 tons' applied load, local failure 
occurs at sections 3 and 5. After this loading the granular structure 
of the soil seems to have been readjusted by the shearing forces and 
displacements so as to have an increased load-carrying capacity. 

In the case of the Ik BP 117 pile, as the applied load is in- 
creased from 75 to 150 tons, local failure has taken place on section 
1 only. As the applied load is increased to 225 tons, failure continues 
on section 1 and is expanded to include sections 2, 3, U, and 7. Upon 
further increase of the applied load to 300 tons, it is evident that 
failure has generally progressed down to and including section U even 
though there is a very slight increase in friction force on the top 
half of the pile. Sections 6 and 8 are experiencing failure while 
sections 5 and 7 are attempting to assume the load. Even at this 
depth it is interesting to note that section 5, the uppermost of these 
two sections, has assumed an appreciably higher part of the load. 

Both of these piles were designed for 150- ton static loads. 

When test-loaded to 150 tons, it is noted that the general shape of 
side frictional resistance was triangular with the base of the triangle 
at the top of the pile and the apex toward the lower sections. In the 
case of the lU BP 117 at a test load of 150 tons, however, the general 
shape of the frictional resistance acting along the length of the pile 
was generally of a rectangular shape down almost to the pile tip. Thee 
difference may be explained by considering the loads required for a 
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unit axial shortening (spring constant) of each of the pile shapes. 

AE 

From the expression K - , and the fact that the length and modulus 

of elasticity of each pile is the same, their spring constants' ratio 
will vary as their area ratio or 3U.Ult/26.19. It can be seen that the 
1U BP 11? is 1,32 times as "stiff" as the lU BP 89 . Since the over-all 
dimensions of the piles, the applied load, and the surrounding soil 
are essentially the same, this relative "stiffness" must account for 
the difference found in the pattern of frictional resistances developed. 

As a consequence of this, it is interesting to note that the lU BP 117 
was nearer total failure at a test load of 300 tons than was the 
lighter lU BP 89 pile. 

In summary, the following factors may be important in con- 
sidering the static behavior of piles in granular material: 

a. The time which has elapsed since driving. In some cases a 
quick condition may have been caused during driving which will dis- 
appear as pore pressures dissipate. Relaxation may also occur. 

b. The general manner in which load is transferred to the 

soil may vary considerably and is generally dependent upon the "stiffness" 
of the pile in resisting axial deformations as related to the strength 
characteristics of the soil. The less "stiff" piles tend to give up 
ttieir loads to the soil at the upper strata while the very "stiff" 
piles tend to carry their load by end bearing. 

c. Failure of piles in granular materials is often not too well 
defined in the case of friction piles. Small local failures along the 
pile may result in a decrease in load carrying capacity, but then the 
soil particles in undergoing shear may assume different packings and 



become denser with the result that their strength has been restored 
or even increased over its original value. The pile load may be 
increased while these readjustments are taking place until a load is 
finally applied which is sufficient to push the pile through the soil 
continuously. With each of these readjustments the soil takes on 
different properties causing the pile to undergo variations in its 
bearing capacities. 

d. Unusual conditions which may be pertinent to a particular 
problem such as a quick condition or hydrostatic uplift. 



As in the case of piles in granular material, some research 
work has also been done with instrumented piles in cohesive materials. 
H. B. Seed and L. C. Reese report in reference (2) that several 6-inch- 
diameter pipe piles from 20 to 22 feet long were driven approximately 
15 feet into a stratum of soft, saturated clay. Electric strain gages 
were installed on one of these piles. Test loadings to failure were 
made at various intervals after driving as follows: 



2. Cohesive Materials. 



Interval After Driving 



Failure Load 
(pounds) 



3 hours 
1 day 
3 days 
7 days 
111 days 
23 days 
33 days 



1070 
2800 
3200 
5U00 
5800 
6100 
6 200 



This data show that the ultimate bearing capacity increased very 
rapidly with time for the first week, but that thereafter the increase 
was very gradual. 
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Most of the load was removed from the pile by side friction 
with only about 10-15$ of the applied load reaching the pile tip. 

The increase in supporting capacity was reflected by an increase in 
the load transfer to soil at all depths as shown by Table D which was 
derived from published data (ref. 26). The pile was divided into 
three-foot sections with section 1 being at the top of the pile. From 
this table it can be seen that the resistance is generally increasing 
with depth for the first 9 to 12 feet of embedment, and then it tends 
to become constant with depth for the final 3 to 6 feet of embedment. 

It is interesting to contrast this behavior with that already described 
for granular materials. 

In this series of tests relationships were established for 
undisturbed soil, and fully remolded soil. Tests of soil samples taken 
next to the pipe wall showed that the strength loss caused by dis- 
turbance of driving the pile was about 6C$to 7C$ of that which the soil 
would have lost due to conplete remolding. The clay next to the pile 
wall subsequently gained strength, and at the end of the test period 
it had a strength of 60$ higher than the strength of the undisturbed 
soil. This increase in soil strength appeared to parallel the decrease 
in excess hydrostatic pressure indicating a relationship between the 
two. For all vaLues of time at all depths the failure of the friction 
pile did not occur in the soil, but at the interface between pile and 
soil. 

In another series of tests made by Vey (ref. 35), with H piles, 
he found that in some cases the friction force reversed itself when the 
applied load was removed and held the stress in the pile. He also 



Table D 



FRICTION FORCE ON SECTION 
(IN POUNDS) 



Test 7 — 33 Days After Driving 



Increment of Total 


0-2000 


lbs. 


2000-1*000 lbs. 


1*000-6000 lbs. 


Applied Load 


Change 


Total 


Change 


Total 


Change 


Total 






(for 




(for 




(for 


Section of 6" Pipe 




2000) 




1*000) 




6000 ) 


Pile and Soil 














Description 














1— Soft sat. clay 


110 


no 


0 


no 


50 


160 


2--Soft sat. clay 


290 


290 


270 


560 


150 


710 


3 — Soft sat. clay 


U80 


L80 


U00 


880 


1*20 


1300 


Ij.— Soft sat. clay 


U 80 


U80 


520 


1000 


500 


1500 


5 — Soft sat. clay 


h 75 


U75 


525 


1000 


560 


1560 


TOTAL FRICTION FORCE 




1835 




3550 




5230 


TOTAL FORCE ON TIP 




165 




U5o 




770 



First Loading to 0-1070 lbs. 



Failure (3 hrs. 


Change 


Total 


after Driving) 


(for 

1070) 


1 — Same soil 


0 


0 


2— 


200 


200 


3— 


0 


0 


h- 


330 


330 


5- 


280 


280 


TOTAL FRICTION FORCE 




810 


TOTAL FORCE ON TIP 




260 



found that adherence would increase between the soil and the pile with 
time, and reconsolidation of the disturbed clay would increase its 
shear strength. 

The findings of these two series of tests regarding the effect 
of pile driving on clay cannot, of course, be assumed to be valid under 
all conditions. Tscheb otari off has shown in reference (3I4) that the 
remolding effect caused by driving pipe piles closed end in 100 feet 
of varved clay can be great. Remolded clay layers can lead to the 
development of negative friction with the result that end bearing 
piles may become overloaded. It is generally understood that displace- 
ment piles should not be driven in highly sensitive clays; therefore 
this part of the problem will not be discussed further. 

In summary, it may be said that the static behavior of piles 
driven in cohesive material may depend upon: 

a. The time which has elapsed since driving. In sane cases 
the increase in bearing capacity for even short periods of time may 
be great; in other cases a loss of bearing capacity may result. 

b. The friction force acting on the pile tends to increase 
with depth. It is believed that the load transfer would decrease with 
depth if the shear strength of the clay were sufficient to accept the 
loads at the upper strata and if there were a motivating force avail- 
able to cause these upper strata to come into intimate contact with 
the pile. Since neither of these conditions usually exists, the load 
transfer is found to generally increase with depth. 

c. Failure of piles in clay may be more dramatic than those in 
granular materials. In the case of a pile not resting on an extremely 
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hard strata such as rock, once the shear strength or adhesion has been 
exceeded, the pile may begin to punch its way through the supporting 
strata and continue downward until the load is reduced or the pile 
encounters a firmer strata. 

B. Dynamic Factors Involved in the Pile-Soil Interaction 

The main obstacle in attempting to relate the penetration 
record obtained during driving of a pile to the static load carrying 
capacity is that of obtaining a reasonably accurate estimate of the 
dynamic resistance to penetration between the pile and soil. The prob- 
lem is the same whether using the wave equation or a dynamic pile 
driving formula. 

The dynamic resistance to penetration is related to one or 
more of the following factors : 

a. The static resistance. 

b. The increase in strength observed to occur in the soil 
due to rapid loading. 

c. Frictional resistance between the pile and soil. 

d. Viscous resistance developed as the pile penetrates the 
soil. 

The basic s trength of the soil, whether due to cohesion or 
internal friction, constitutes the static load capacity of the pile 
and soil. At least, this represents the limiting value since the pile 
may not develop the full shear strength of the soil when acting as a 
friction pile. The value of the static resistance can be thought of 
as the base value of the dynamic resistance in the idealized case 
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where the strength properties of the soil are not changed by the 
process of pile driving. This is rarely the case, and usually the 
driving of the pile acts to decrease or increase the natural strength 
of the soil. 

Several investigators have studied the effect of rapid loading 
times on the compressive strength of the soil (ref. U, 11, 25, 37), 
and some discussion of the results of such investigations appears 
below. It seems logical to assume that in the early time intervals 
of pile penetration most of the dynamic resistance is a result of 
this type of phenomenon, the magnitude of which would be directly 
related to the static strength of the soil. 

Frictional resistance in nearly all cases develops between the 
pile sides and soil as a result of the penetration of the pile. The 
nature and magnitude of this frictional resistance depends upon the 
type of soil, moisture content, type of pile and pile material. If 
the adhesion or friction acting between pile and soil is greater than 
the shear strength of the soil, the frictional resistance may be de- 
veloped within the soil itself. The dynamic frictional resistance 
would be of about the same value as the static resistance since the 
dynamic coefficient of friction between the pile and soil materials 
is abait equal to the static coefficient of friction (ref. 2U). Fric- 
tion as a dynamic factor will be neglected since there is some doubt 
as to the extent that actual relative movement takes place between the 
pile and soil at the pile-soil interface. 

No published work has been found that deals with viscous 
forces as developed by penetration into soils. One investigator 
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attributes part of the influence of high rates of strain on the increase 
in compressive strength of soil to viscous action (ref. 37). Since 
viscosity is defined as a property of internal friction (ref. 2), it 
will be considered as included as one of the factors contributing to 
the increase of compressive strength experienced by a soil subjected 
to high rates of strain. 

During the driving of a pile, it is assumed that the resistance 
developed at the tip and along the sides can be expressed in terms of 
the resistance developed under static loading. This is to say that 
the dynamic resistance to driving is composed of the static resistance 
plus an increment of resistance that develops under dynamic loading 
and is expressed as a percentage of the static value. The problem is 
to determine what form the relationship between static and dynamic 
resistance takes and how to quantitize it. 

Another problem, which will be taken up later, is how much of 
the static resistance is present during driving. It caanot be generally 
assumed that all of the static resistance is present during driving, 
and Chellis points out ways in which this may vary (ref. 6). On the 
other hand, dynamic resistance may develop in soil strata that offer 
only very small or negligible static resistance. 

In the computer solution, the total resistance to driving at 
the pile point is represented by the expression: 

R = (D m - Djjj) K' (1 + Jv) (15) 

where (D m - D^) is the "elastic” strain of the soil with 

a limiting value of Q, 

K 1 is the spring constant of the soil. 
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v is the pile tip velocity at any instant, and 
J is the damping constant. 

This relationship is shown by Smith (ref. 28). 

Two questions arise with respect to this mathematical representa- 
tion of the dynamic resistance to driving: 

(a) Does this equation approximate the actual form of the pile 
resistance; i.e., does the resistance-time relationship approximate 
that of the actual case? 

(b) What values of J should be assigned for different soil and 
driving conditions? 

Smith proposed the use of a value of 0.15 for J; however, since the 
purpose of his published work did not include the correlation of re- 
sults with load tests, he gave no basis for the use of this value. Some 
preliminary attempts to correlate wave equation results with the results 
of pile load tests and driving records using values of 0.15 and 0.05 
for J and J' indicated that this value for J was not appropriate for 
clay soils. 

In an attempt to answer the questions posed above, a search of 
the literature was made for any experimental work that might lead to 
a confirmation of the nature of the dynamic resistance as well as to a 
determination of more appropriate values for the coefficients J and J'. 

Several research programs over the past 15 years have investi- 
gated the influence of time of loading, or strain rate, on the compres- 
sive strength of soils. All such projects have used rapid loading in 
triaxial and unconfined compression tests, or dynamically loaded footing 
models as the method of investigation. The application of results fran 
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such tests to pile driving must be approached with caution and should 
be considered only as indicative of the general nature of the dynamic 
resistance to pile driving ■ and of the general order of magnitude for 
values of the coefficients of damping. 

Before describing the tests made by individual investigators, 
several terms will be defined: 

(1) Time of loading: The time from the beginning of a test 
to that at which the maximum force is achieved. Figure (21 ) shows 
this quantity graphically. 

(2) Strain rate: The rate, expressed in per cent per second, 

of strain imposed on the sample at the time the compressive strength 
is reached. 

(3) Transient load: A dynamically applied load in an uncon- 

fined compression or triaxial test that reaches a maximum in times of 
less than 30 seconds and as short as 0.005 second. 

(U) Static test: For comparative purposes, a test that is made 
with a time of loading of ten minutes or more, unless otherwise indicated. 

(5) Strength ratio: The ratio of the compressive strength for 

a given time of loading to the strength for a static test. 

(6) Strain-rate effect: Same as strength ratio. 

2. Work at Harvard 

In a test program sponsored by the Panama Canal in 19US, 

Casagrande and Shannon made a series of static and transient unconfined 
and triaxial tests on soft to medium clays and a soft rock (ref. I:). 
Triaxial tests of the vacuum type were also made on a dry sand. Times 
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of loading for the transient tests ranged from 0.01 seconds to 30 
seconds . 

In all the tests on clays, the transient compressive strength 
was greater than the compressive strength as determined by static tests. 
The transient strengths, taken at a time of loading of 0.02 second, 
were from l.U to 2.6 times the static compressive strength. The tests 
also showed that the increase of the fast transient compressive 
strength over the static strength is greatest for specimens at the 
highest water content and least for the specimens at lowest water 
content. (Figure (17-A) is a summary of the average results from the 
clay tests. 

Vacuum-type triaxial tests were performed on a clean, medium 
sand with a void ratio of 0.62 and at a minor principal stress of 
0.3 kg/cm^. The test results indicated that the compressive strength 
of a dry sand in fast transient tests is about greater than the 
static compressive strength. 

3. Work of Whitman 

Whitman measured the increase in compressive strengths of 
soils using a triaxial machine, and testing at high rates of strain 
(ref. 37). Table E shows his results for several types of cohesive 

soils. The last column of this table (S-R) is the ratio of the 

failure load at a strain rate of 1000 per cent per second to the 

failure load at a strain rate of 0.03 per cent per second. It is 

observed that the strain-rate effect varies from 1.3 to 2.0 for these 
soils. In his tests, Whitman achieved times of loading as short as 
.005 second. Confining pressures up to 100 p.s.i. were used. 



Table E 



INCREASE OF COMPRESSIVE STRENGTH 
IN COHESIVE SOILS DUE TO HIGH STRAIN RATES 

(from ref. 37) 



1 . 



2 . 

3. 

U. 



5 . 



6 . 

7. 

8 . 



9 . 



10 . 



DESCRIPTION 

Normally consolidated, sensitive 
ocean sediment, undisturbed 

Very plastic clay, remolded 

Plastic clay, remolded 

Medium soft, slightly sensitive 
clay, undisturbed 

Slightly plastic clay loam, 
remolded 

Plastic clay, remolded 

Moderately sensitive silty clay, 
undisturbed 

Impervious compacted fill 
Tough compacted fill 
Stiff, dry clay, undisturbed 



P.I. 


P.L. 


w 


-2H. 


S-R 


63 


h9 


92 


0.3 


2.0 


27 


38 


h8 


7 


1.6 


17 


11 


16 


10 


1.7 


2U 


26 


27 


10 


2.0 


23 


22 


21 


13 


l.U 


27 


38 


kh 


15 


1.7 


21 


22 


35 


22 


1.6 


17 


11 


12 


25 


1.8 


Ui 


21 


26 


35 


2.0 


23 


30 


20 


250 


1.3 



In comparing results of tests on clay with and without confining 
pressure. Whitman observed two different time effects. The first effect 
can be likened to that encountered in an extremely viscous fluid, and 
is associated with a continuous, plastic deformation of the soil. The 
second effect is associated with the formation of discontinuities such 
as shear planes and splits, aid occurs when the soil sample is tested to 
failure without (or small values of) confining pressure. From these 
results, he concluded that when cohesive samples were subjected to a 
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lateral pressure, the strain-rate effect was the result of viscous 
phenomena, and was the same regardless of the magnitude of strain. 

The same effect occurred in very plastic soils regardless of confining 
pressure. 

In dry or moist sands an increase in compressive strength of 
10-15 per cent was observed. Only limited tests were conducted with 
saturated sands, and figure ( 18 ) shows results for saturated Ottawa 
sand and a well graded sand. Whitman observed that the strain-rate 
effect in saturated sand must result from differences between the pore 
water pressures in slow and rapid tests. When sand is tested at 
constant volume, pore water migrates to the central part of the sample. 
As speeds of deformation increase, more energy must be expended to 
overcome resistance to the flow of water. 
k. Work of Seed and Lundgren 

Seed and Lundgren investigated the effect of transient loadings 
on the strength and deformation characteristics of saturated sands 
(ref, 25). Their tests were also made with the triaxial apparatus, and 
three different categories of tests were conducted: 

a. Static tests in which the loading times were about 
10 minutes; 

b. Slow transient tests with loading times of about 
U seconds; 

c. Rapid transient tests with loading times of about 
.02 second. 

In the rapid transient tests a constant rate of deformation of U0 inches 
per second was used. All tests were conducted wi th a confining pressure 
of 2 kg/sq.cm., and the samples were tested at various void ratios. 
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Results of the tests on dense, fine sand indicated that the 
compressive strength of a sample subjected to rapid transient loading 
was about UO per cent greater than that obtained by static drained 
tests. Due to the high rate of loading the rapid tests were actually 
undrained tests, since drainage could not occur within the very short 
time of the test. About half of the increase in strength is attributed 
to a dilatancy effect and half to the high rate of loading alone. 

It is a well known fact that when a dense sand is sheared it 
increases in volume. If shear occurs under conditions that prevent 
the movement of pore water, such as in the undrained triaxial test or 
in very rapid loading, the increase in volume will result in a decrease 
in pore pressure and an increase in strength, as shown by the equation 

s = ('5 - u)tan0 

in which the value of u (pore pressure) is decreased. These results 
are summarized in figure (19) and apply to fine and coarse sands far 
very short loading times. 

If the results of the tests made by Seed and Lundgren are 
applied in equation (15) the following results are obtained: 

Dynamic Resistance (R^) a 1.U0 x Static Resistance (Rg) 
or, Rp = Rg(l + JV) 

where J is the damping factor and V the velocity of deformation, 

1.U0 Rg - R s (l + JV) 

T _ 1.U0-1.00 

j = 

From figure (20), the load curve from one test, the velocity of 
deformation is estimated to be about 35 inches per second, and 



58 



substituting this value in the last equation gives 

J - .133 

for the case of rapid loading in saturated sand. This value agrees 
closely with the value of 0.15 for J suggested by Smith (ref. 28). 

5. Work of Cunny and Sloan 

Some preliminary tests were conducted by the Waterways Experi- 
ment Station, Corps of Engineers, Vicksburg, Mississippi in which small 
footings were dynamically loaded in a test tank filled with clay or 
dry sand (ref. 11). In these tests a loading machine capable of 
applying a load in a time of from 3 to 150 milliseconds was programmed 
to apply a preset load diring a predetermined loading time. As a 
control, static tests were also performed on the same soil under 
identical conditions except for loading time. 

Dynamic and static tests were conducted on a sand compacted to 
a relative density of 96 per cent. Figure (21) shows the relationship 
between load and deformation for the static and one dynamic test on a 
9-inch-square footing. The load time curve for the dynamic test is 
shown in figure (22). From this curve it is seen that the loading 
time was 11 milliseconds and the maximum load reached was 3500 pounds. 
Since the maximum value of the load was limited by the setting of the 
machine, it is not necessarily the full strength of the soil under 
dynamic loading. 

By relating the values of static and dynamic loading at equal 
deformations, it is possible to arrive at an indication of the value 
of the damping factor under these conditions of loading. At a deformation 
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of 0.10 inch the static load is found to be 1600 pounds and the dynamic 
load 2650 pounds from the load-deformation curve of figure (21). The 
velocity of deformation at this point is estimated to be about UO inches 
per second, or 3.33 feet per second. The corresponding value of J 
using equation (15) as before is found to be 0.197. 

6. Comparison of Dynamic Tests with Computer Solutions 

In order to determine whether the dynamic tests described above 
are qualitatively comparable with the equations used for computation 
of the total resistance, an analysis has been made of several sets of 
computer results in the same terms as the data obtained f rcm the 
laboratory results. This analysis could only be carried to the point 
of maximum deformation, since the program was established to stop at 
this point. 

Figures (23) and (2h) are plots of the computer results for 
pile tip resistance-time, deformation- time, and load-deformation 
relationships for a 53-lb. H pile UO feet long and a concrete pile 36 
inches in diameter and 1?6 feet long. The H pile is driven with a 
Vulcan #1 hammer and the concrete pile is driven with a Raymond U/0 
hammer. The tip resistance-time curve for the H pile, figure (2ii), 
compares closely in shape with same curve for the dynamic laboratory 
test of the 9-inch-square footing shown in figure (22). This would 
tend to indicate that the load-time relationship resulting from 
equation (15) is a fair representation of the actual dynamic load 
mechanism, at least with respect to a dynamic test on a small footing. 

The pile tip velocity-time curves in figures (23) and (2h) 
show a velocity of about 2 feet per second at the time of maximum 
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resistance for the UO-foot H pile and a velocity of about one foot 
per second for the 176-foot concrete pile. If it is accepted that 
most of the increase of dynamic resistance over the static resistance 
is due to the increase in compressive strength of the soil due to 
rapid loading, and if this increase amounts to a factor of 1.0 for 
days and O.U for sands, then the damping coefficient must have a 
value such that when multiplied by the velocity of deformation at 
the time of maximum load results in values of 1.0 for clays and 
O.U for sand. Figures (23) and (2U) indicate that tip velocities at 
the time of maximum resistance are between 1.0 and 2.0 feet per second, 
and based on these values, the damping factor (J) would have a value 
of 0.5 to 1.0 for clay and 0.1 to O.U for sand. 

If load vs. deformation is plotted from the information 
available in the curves of load and deformation vs. time for the 
12 BP 53 pile shown in figure (23), a curve shown in figure (25) is 
obtained. This curve is very similar in shape to the curve of dynamic 
load vs. deformation from rapid load tests performed on sarid in the 
laboratory as shown in figure (21). The fact that these curves are 
similar in shape tends to indicate that the variation of load with 
deformation in the wave equation solution is of the same fcrm as 
that occurring in laboratory load tests. From this it is concluded, 
tentatively, that the computation of total driving resistance at the 
tip by the equation 

R ** (Dp - D^) K« (1 + Jv) 

is a valid representation of the actual physical process of rapid 
penetration of the soil by the pile. 
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7. Ground quake, -which has been defined as the limit of elastic 
deformation of the soil, is normally thought of as a static concept, 
but it is an essential element in wave equation computations. The 
elastic ground compression occurring during static loading is not 
necessarily the same that occurs during pile driving, and no attempt 
has been made to relate the two. Therefore, quake as used in this 
investigation is a dynamic property. 

A method of measuring quake during pile driving is described 
by Chellis (ref. 6). Essentially, this method measures the total 
elastic compression of the pile and soil. Pile compression alone is 
computed from the formula ** R u L/AE, where is the ultimate 
bearing load on the pile, L is the length of the pile, A is the cross- 
sectional area, and E ± s the modulus of elasticity. This value is 
subtracted from the measured compression of pile and soil to obtain 
the value of quake. This method of measuring quake assumes that the 
pile is in compression simultaneously along its entire length, which 
is somewhat erroneous in view of the longitudinal stress wave concept. 

In extensive studies made of pile driving, Hiley made measure- 
ments of quake and reported that quake varies considerably for different 
conditions of ground and may average 0.1 inch in firm gravel, up to 
0.2 inch in firm clay and may have values as high as 0.5 inch in 
peaty material (ref. 17). With respect to driving conditions, Hiley 
summarized his observations as follows: 

Driving Conditions Quake Stress in Pile 
Easy driving 0.0 - 0.10 500 psi 

Medium driving 0.10 - 0.20 1000 psi 

Hard driving 0.10 - 0.30 1500 psi 

Very hard driving 0.0$ - 0,20 2000 psi 
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It was also stated that quake may be as much as double these values 
if soft ground is immediately below the pile tip. These values 
for quake are used as a limiting range for investigations of load 
tests. 
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Chapter VI 

CORRELATION OF WAVE EQUATION RESULTS WITH PILE LOAD TESTS 

A. Methods of Approach and Effects of the Various Parameters 

Although there are many load tests published, there are com- 
paratively few which have not left out some vital part of the in- 
formation required for a thorough evaluation of the total pile driving 
and loading process. Partly due to this limitation, two general 
methods of approaching the problem of correlating wave equation 
computations with load test data were tried. 

The first method was particularly applicable to those tests 
for which very complete information was published. In these cases 
it was thought that most of the variables could be fairly accurately 
evaluated so that the investigation could be limited to a trial and 
error approach in which the principal unknowns could be varied over a 
wide range of values and tried in various combinations . Generally 
this means that ground quake and damping constants and side friction 
were varied with all other factors kept constant for the case in 
question. 

The second method of approach was adopted to try and derive 
as much information as possible from published load test data which 
were not reported as completely as those investigated by the first 
method of approach. It also considered desirable to run a series of 
data to study more fully the effects of the various parameters on the 
resistance-set relationship for particular piles and hammers. It 
was thought that both of these objectives could be met by selecting 
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a particular pile type for which a number of reasonably good test 
data were available. Furthermore, it was desired that the pile meet 
the following requirements: 

(a) Be constructed of a material of which definite information 
was available as to its properties. 

(b) Be constructed of a material which of itself would least 
affect the soil around it and vice versa. 

(c) Have uniform properties throughout its length. 

(d) The pile as driven would be required to sustain the static 
loads to be imposed upon it. 

The reason for these requirements was to eliminate as many of 
the unknowns from the problem as possible. The steel H pile and steel 
pipe pile each satisfied these requirements better than other types. 
Because there were more usable load tests readily available for steel 
H piles, particularly the 12-inch H bearing pile weighing 53 pounds 
per foot, this shape was selected. It had the added advantage that 
most of the records available indicated that one type of hammer had 
been used, the Vulcan Number One; however, tests were available with 
other size hammers so that a further check of the method could be made. 

The main effort has been directed toward establishing correla- 
tion for each load test case within a limiting range of values of each 
of the variables involved. Correlation of wave equation computations 
with the actual pile load test involves the evaluation of five different 



variable factors. These are: 
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(a) Ground quake; 

(b) Point damping factor; 

(c) Side damping factor; 

(d) Distribution of resistance between point and sides of 
the pile; 

(e) Per cent of static resistance acting during driving. 

These factors have been discussed in above sections with the exception 
of the last named variable, the per cent of static resistance acting 
during driving, which will be discussed below. 

In correlating theoretical calculations with load tests and 
the corresponding driving records, all variables except one must be 
held constant while varying one factor at a time. From practical 
considerations it was not possible to treat all five factors alike as 
variables. Since little is known of the nature of the side damping 
factor and since the effect of side damping in most cases is relatively 
small compared to point damping, the side damping factor was eliminated 
as a separate influence by assuming it equal to one third of the point 
damping factor for most cases. Smith indicated his belief that side 
damping would be less than point damping for a given soil, and he 
recommended a value which was one third of the point damping coeffi- 
cient (ref. 28). 

Although up to five different distributions of resistance were 
investigated for each correlation case, this factor was assumed known 
when it came to the actual correlation using the first method. For the 
other method, distribution was treated within a range of values. The 
various distributions of resistance used for investigation are: 
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(a) End bearing 

(b) 75 per cent end bearing and 25 per cent side friction 

(c) 50 per cent end bearing and 50 per cent side friction 

(d) 25 per cent end bearing and 75 per cent side friction 

(e) All side friction 

Since the pile is generally divided into sections 5 to 10 feet long 
for purposes of calculation, the side friction acting on the lower 
section of the pile is indistinguishable from the point bearing resist- 
ance, except for damping, and adds to it for all practical purposes. 

Three methods of distributing side friction along the embedded 
length of the pile have been used. The nethod used most, and the 
simplest, assigns a uniform value of resistance for a unit length of 
pile regardless of depth. This method is termed rectangular distribu- 
tion. In another method the side friction is distributed in a manner 
so that the unit resistance is directly proportional to the depth, and 
is called triangular distribution. Both of these methods are included 
in the computer program, and selection of the applicable method is 
made by use of a control statement in the progran. Except for very 
long piles, the difference in effect on the resistance-set relationship 
of the two methods is small, so the first method was generally used. 

The third method of distributing frictional resistance is by estimating 
the actual resistance developed in each of the soil strata that the 
pile penetrates. 

The per cent of static resistance acting during driving is a 
factor that is influenced by time effects. In a sensitive clay, pile 
driving remolds the soil resulting in a reduction in strength during 
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and immediately after driving of the pile. If the load test is con- 
ducted after the soil has had time to regain its strength, the value 
of the static resistance should be reduced when making the pile-driving 
analysis. Another type of time effect may result when a pile is driven 
into a sand layer through overlying layers of soft or medium clay. 

When the pile is test loaded, resistance initially developed in the 
clay may be transferred by plastic flow to the tip and lost from the 
standpoint of static load-carrying capacity. During driving, however, 
resistance in the upper clay layers will be present and will appear 
to result in heavy damping. Another effect is the liquefaction of 
saturated silts due to a quick condition being created by the driving 
action and resulting in greatly reduced resistance to penetration of 
the pile. The consideration of the effects of such factors in the 
calculations is very approximate and without complete soil data cannot 
be accounted for. 

In correlating wave equation computations with load tests in 
which sufficient data was available on which to base a reasonable 
estimate of the distribution of load between point bearing and side 
friction, the approach was to determine the various values of quake 
and damping that produced a resistance-set curve that included the 
point representing the ultimate load and final penetration of the test 
pile. Hypothetically, there would be an infinite number of such 
curves, but a limiting range of values was used for these variables 
that was determined on the basis of Smith's work and the experimental 
work on the dynamic properties of soil discussed above. The range of 
values for quake and damping investigated were: 
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(a) Quake 

(1) sand 0.05 to 0.20 

(2) clay 0.05 to 0.30 

(b) Damping 

(1) sand 0.10 to 0.20 

(2) clay O.UO to 1.00 

These ranges of values were not rigidly adhered to, and in 
some instances outside values were investigated. The load test in- 
formation, input data used for the computer program for each load 
test, and results of the correlation attempts are included in tables 
F to BB. The accompanying figures are the ultimate resistance versus 
set (in blows per inch) curves drawn from the computer results and 
used as the basis for correlation. 

Similar curves used as the basis of correlation by the second 
method are included in appendix E. The length of piles used for the 
calculations was chosen to generally fit the load test data available, 
but since the effect of length was not clearly established the choice 
of U0-, 80- aid 120-foot piles was somewhat arbitrary. Calculations 
were made for 12 BP 53 piles in these lengths using a Vulcan Number One 
hammer with a hardwood capblock, a pile cap weighing 700 pounds, for 
all permutations of the ranges of ground quake, point and side damping 
factors and side friction distributions of a rectangular shape. For 
each combination of the variables, set in blows per inch was computed 
for each value of ground resistance ranging from liO to and including 
280 tons in 20- ton increments. The results were plotted so that on 
any graph there are five curves, one for each percentage of side 
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resistance assumed. For any one curve of ultimate resistance vs. set, 
all other variables are constant. 

These results show the following generalized effects of 
changing the variables indicated below: 

(a) Ground quake. Increasing ground quake increases the set 
in blows per inch, i.e. harder driving, for a given ultimate ground 
resistance for all values of side resistance including the point 
bearing case. 

(b) Increasing both the point and side damping increases the 
set in blows per inch for a given ultimate ground resistance for al . l 
values of side friction, including zero. 

(c) Side distribution-- increases in per cent of side friction 
cause dramatic decreases in set in blows per inch for a given ultimate 
ground resistance, 

(d) Increased pile length — here the results are mixed. In the 
zero friction or end bearing case, the ultimate resistance tends to 
become the same for 1*0-, 80-, and 120-foot cases after about 2h blows 
per inch. Before that point is reached, the 80- and 120-foot piles 
have essentially the same curve of ultimate resistance vs. set, and 
the UO-foot pile for a given resistance within this range tends to have 
an appreciably decreased set in blows per inch. For about $0$ friction, 
the results may be practically identical for all three lengths. For 
higher values of friction the longer piles may have decreased set for 
given resistance. 

(e) Increased pile weight may result in decreased set in blows 
per inch for a given ultimate resistance; however, in the lower ranges 
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of resistance the difference may be slight. For exarrple, for a side 
resistance of 25 per cent and all factors the same far two U5-foot 
H piles driven with a Vulcan 0 hammer, the difference in set for a 
lU BP 89 and a lU BP 117 is less than 2 blows per inch. The values 
are 8 and 6.6 respectively. 

(f) Increased hammer weight, but with the same stroke, results 
in decreased set in blows per inch for a given resistance and all 
other factors constant. 

(g) A softer capblock results generally in increased blows 

per inch for a given resistance, but the effect is really quite variable. 
Comparisons were made for an ultimate ground resistance of 100 tons for 
the following piles driven with a Vulcan Number One hammer: 

1) Concrete Piles — 18 inches square with lengths 
of l»0, 80, and 120 feet. It was assumed that head packing 
would be used, but this same soft material was assumed for 
the capblock also and the results compared with capblocks 
of hardwood and micarta. There was no appreciable differ- 
ence in sets for any of these conditions for any of the con- 
crete piles. 

2) Wood Piles — for the 80- and 120-foot lengths, the 
difference was slight for the three capblocks. For the UO- 
foot pile the difference between a micarta and a hardwood 
capblock were slight, but the very soft capblock made about 
10 blows per inch difference in the results. 

3) Steel piles— the capblock had essentially no effect 
on the 80- and 120-foot length piles, but the very soft cap- 
block made about 2 blows per inch difference with the UO-foot 
pile. 
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The effects produced by varying some of these parameters are 
not surprising, but in some instances they may be, especially if one's 
judgment is based on results produced by popular dynamic formulas. 

For example, the results produced by the Hiley formula may indicate 
very slight differences between the cases when side friction is 
present and when it is totally absent. On the other hand, it is 
interesting to note that when side friction is a large percentage of 
the ultimate resistance to driving, the set vs. ultimate resistance 
curves may tend toward the curve produced by the Engineering News 
formula. The foregoing results produced by the wave equation method 
of analysis would seem to bear out its general suitability under widely 
varying conditions of driving equipment, pile type, and ground conditions. 

As can be seen by examining the following case studies, some 
of the data upon which to base correlation is rather complete, and 
some is not. The cases studied include various size piles, assorted 
pile driving equipment, in soils varying from gravel to clay, piles 
constructed of concrete and steel, different ratios of embedded length 
to driven length, and for varying combinations of soil strata. It is 
apparent that for many of the cases correlation can be obtained in 
several ways by various combinations of values for ground quake, point 
and side damping, and side friction. Because these factors are not now 
known with a high degree of certainty, the problem becomes one of 
searching for a pattern within the likely ranges of the variables for 
the various types of soils encountered. 

In some of the following cases it was necessary to estimate 
the values of Q and J which produce correlation, since the range of 
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values used in the computations did not extend far enough to include 
the load test point within computed load-set curves. The values for 
Q and J so estimated are based on the curves nearest to the load test, 
and the trend of variation between curves. It is believed that the 
estimated values are reasonably accurate considering the nature and 
accuracy of the basic data. 

Before examining the case studies in detail, the manner of 
identifying the set versus resistance curves is noted: 

Curve A — 0$ friction, i.e., an end bearing pile with no 
resistance to driving offered by the soil along 
the pile sides. 

Curve B — 25 % friction, i.e., 75 % of resistance to driving 
offered by the soil at the pile point and 2 5 % 
along the pile sides. 

Curve C — 50$ friction, i.e., 5C$ of resistance to driving 
offered by the soil at the pile point and £0$ 
along the pile sides. 

Curve D — 75 % friction, i.e., 25 % of resistance to driving 
offered by the soil at the pile point and 75 % 
along the pile sides. 

Curve E — 100$ friction. All resistance to driving 
distributed along the pile sides. 
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TABLE F 



Case Number 2001 
Reference: (18) 

Pile: Prestressed concrete, 36 in. diameter, 176 ft. long 

Hammer: Raymond 3/0, 12500 lb. ram, 2h in. stroke 

Embedded length: UO ft. 

Soil: Soft to medium clay with point in stiff clay. 

Final penetration: 20 blows per inch 



Load test: Tested to failure at 190 tons including the w eight of 

pile, tested one week after driving 



Input data: 



W(1) 


12,500 


Ram 


W(2) 


8,800 


(includes follower and driving cap) 


W(3)-W(20) 


U, 955 


(based on 9.775 ft. segments) 


V 


10.15 


(ram velocity at impact for 2-ft. 
drop) 


K(l) 


10,600,000 


(capblock spring constant) 


K(2) 


2,180,000 


(cushion and pile) 


K(3) - 
K(19) 


20,750,000 


(Pile) 


K«(m) 


Soil spring constant depends on value of 
ultimate resistance which is varied. 



Distribution of resistance: 50$ end bearing and 50$ friction, based 

on computation of point capacity by the 
formula: P = 9 AC . Friction considered 
to have a rectangular distribution. 



Results: See figures (27) and (28) 

Correlation is obtained for the following values 
of Q and J : 

Q J J 1 

0.1 0.9 0.3 



0.2 



0.3 



0.1 
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TABLE G 



Case Number 2002 
Reference: (l8) 

Pile: Prestressed concrete, 36 in. diameter, 176 ft. long. 

Embedded length: UO ft. 

Hammer: Raymond b/0, 15,000 lb. ram, 3b in. stroke. 

Soil: Soft to medium clay with point in stiff clay. 

Final penetration: Refusal, last 3A in. at 99 blows. 

Load test: Tested to failure at 230 tons including the weight of pile. 

Input data: W(l) 15,000 

W(2) 8,800 

V 12.08 

(Rest of input data same as for case 2001) 

Distribution of resistances Estimated between 75 % and 50$ point bearing. 

Results: See figures (29) and (30). 

Correlation is obtained for the following values 
of Q and J; 



Q 






0.2 



J 

0.8 



0.27 




I 
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TABLE H 

Case Number 2003 
Reference: (18) 

Pile: Prestressed concrete, 36 in. diameter, 176 ft. long. 

Embedded length: UO ft. 

Hammer: Raymond 3/0, 12,500 lb. ram, 39 in. stroke 

Soil: Soft to medium clay with point in stiff clay 

Final penetration: 10 blows per inch 

Load test: Tested to failure at 165 tons including weight of pile. 

Input data: W(i) 12,500 

V 12.9 

(Rest of input data same as for case 2001) 

Distribution of resistance: Estimate 5Q& end bearing and 50# friction 

Results: See figures (31) through (3U). 

Correlation obtained for the following values 
of Q and J; 



Q 


J 


j« 


0.1 


1.0 


0.33 


0.2 


0.5 


0.17 




I 



I 







* 
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TABLE I 



Case Number 200U 
Reference: (18) 

Pile : Prestressed concrete , 

Embedded length: lj.3.5 ft. 



36 in. diameter, 80 ft. long. 



Hammer: Raymond 3/0, 12,500 lb. ram, 39 in. stroke. 

Soil: Medium and stiff clay with point in very stiff clay. 

Final penetration: 20 blows per inch. 



Load test: Tested to failure at 365 tons. 



Input data: W(l) 

W(2) 

W(3) - 

V 

K(l) 
K(2) 
K(3) - 
e l 

e 2 



12,500 

8,800 

W(10) 5,070 

12.9 

10,600,000 
U, 020, 000 
K(9) 20,300,000 

0.8 

0.5 



Distribution of resistance: From computations based on load test and 

soil data, it is estimated that 87$ of the load was 
carried by the pile point at failure. Correlation should 
be made for distributions between 100 per cent and 75 per 
cent of the load carried by the point. 

Results: See figures (35) through (37) 



Correlation is obtained for the following values 
of Q and J : 



Q J J* 



0.1 0.1;5 0.15 

.05 0.6 0.2 

Remarks: In this case, a hole 2k ft. in depth and lj.2 in. diameter 

and below that for a depth of 27 ft., a hole 35 in. 
diameter was pre-excavated for the pile. The effect of 
the 35 in. pre-excavation is uncertain, but appears to 
have considerably reduced the dynamic resistance to driving. 
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TABLE J 



Case Number 2005 
Reference: (l 8 )^Test 7-1 

Pile: Prestressed concrete, 36 in. diameter, 68 ft. long. 

Hammer: Raymond 3 / 0 , 12,500 lb. ram, 39 in. stroke 

Embedded length: 31 ft. 

Soil: Stiff clay. 

Final penetration: UO blows per inch. 

Load test: Tested to failure at 260 tons. 

Input data: Same as case 200U, except: 

W(3) - W(9) h,930 

K(3) - K( 8 ) 20,900,000 

Distribution of resistance: 

From calculations based on load test and soil data it is 
estimated that b3% of the load was carried by the point 
of the pile failure. For correlation, 50$> of the load is 
considered to be carried by the point. Friction is con- 
sidered to be distributed over the lower 31 ft. of the pile. 

See figure ( 38 ). 

Correlation is obtained for the following values 



of Q and J : 






Q 


J 


j» 


0.10 


1.0 


0.33 


0.20 


0.5 


0.17 



Results : 




I 

I 
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TABLE K 



Case Number 6006 
Reference: (l) 

Pile: Prestressed concrete, 5U in. diameter, 96 ft. long. 

Hammer: Raymond 3/0, 12,500 lb. ram, 39 in. stroke, wood capblock. 

Embedded length: 60 ft. 

Soil: Soft, medium and stiff clay and silty clay with some sand 

with pile of point in dense, fine, gray sand. See figure 
(39) for soil profile. 

Final penetration: l;5 blows last 3/U in. (60 blows per inch). 



Load Test: 



Input data: 



Tested to 1*20 tons, failure had not occurred. Bayliss cal- 
culated a minimum ultimate load capacity of 1*98 tons based 
on load test data and soil characteristics using static 



formulas (ref. l). 




W(l) 


12,500 


W(2) 


6,000 


W(3) - W(12) 


6,200 


V 


12.9 


K(l) 


2,715,000 


K(2) 


2,025,000 


K(3) - K(ll) 


27,200,000 


e l 


0.50 


e 2 


0.50 


T 


0.00033 



Distribution of resistance: 



Based on Bayliss' computations (ref. 1) it is estimated 
that 50 per cent of the load was carried as erri bearing 
with friction distributed over the lower 60 ft. of the 
pile. 



Results : 



See figures (39) through (U2). 

Correlation is obtained for the following values 
of Q and J: 



Q 


j 




0.1 


0.15 


0.05 


0.15 


0.10 


0.05 


o.q5 


.2 


. 05 extrapolated 
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I 

TABUS L 



Case Number 6007 



Reference; (19) 

Pile; Pipe pile, 18 in. diameter, 3/8 in. wall thickness, 80ft. 

long. 

Hammer; Vulcan No. 1, 5,000 lb. ram, 36 in. stroke. 

Embedded length; 

Soil; Predominantly plastic clays to a depth of 80 ft. underlain 

by uniformly graded fine sand having a relative density of 
approximately 80 per cent. Pile tip stopped 5 ft. above sand. 

Final penetration; 2 blows per inch. 

Load test; Tested to failure at 81 tons. 



Input data; 


W(l) 


5,000 




W(2) 


1,000 




W(3) - W(12 ) 


562 




V 


12. i; 




K(l) 


2,715,000 




K(2) - K(ll) 


6,520,000 




e l 


0.50 




e 2 


1.0 




T 


0.00025 



Distribution of resistance; Load carried as friction. 
Results; See figure (U3). 

Correlation is obtained for the following values 



of Q and J ; 






Q 


J 


J' 


0.1 


0.15 


o.o5 




I 

I 
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TABLE M 



Case Number 6008 



References (19) 

Pile: Steel pipe pile, 18 in. diameter, 3/8 in wall thickness, 

96 ft. long. 

Hammer: Vulcan No. 1, 5,000 lb. ram, 36 in. stroke. 

Embedded length: 75 ft. 

Soil: Plastic clays to a depth of 80 ft. underlain by a uniformly 

graded fine sand having a relative density of approximately 
80 per cent. The pile tip was driven about 5 ft. into the 
sand layer. 

Final penetration: 226 blows for the last 6 inches (37.7 blows/in.) 

Load test: Tested to failure at 2I4.U tons. 



Input data: W(l) 

W(2) 

W(3) 

V 

K(l) 

K(2) 

e l 

e 2 

T 



W(12) 


5.000 

1.000 
6 76 


K(ll) 


12. U 

2,715,000 

5,U20,000 




0.50 

1.0 

0.00025 



Distribution of resistance: 



Based on two different load tests, one with the pile tip 
in sand and the other (case 6007) with the tip 5 feet 
above the sand layer, it was determined that 159 tons 
or 65 per cent of the load was carried by the point (the 
lower 5 ft.). It is assumed the frictional load is 
distributed over the embedded length. 

Results: See figure (UU). 

Correlation is obtained with the following values 
of Q and J: 

Q J J' 

0.05 0.10 0.03 

0.10 0.05 0.05 (estimated) 
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TABLE N 



Case Number 6011 

Reference: (20) 

Pile: lU-inch H pile weighing 73 lb. /ft. and 8l ft. long. 

Hammer: Vulcan OR, 9,300 lb. ram, 39 in. stroke. 

Embedded length: 80 feet. 

Soil: Alternating strata of silts, sandy silts, and silty sands 

with some interspersed clay strata for a total thickness of 
about 50 to 60 feet. Clean sands from UO to 60 feet thick 
lie beneath the silts. The pile penetrated 32 feet into 
the sand layer. 

Final penetration: 28 blows per foot (2.3 blows/inch) 



Load test: The pile buckled under a test load of 3hh tons; at this 

point the pile had reached a gross settlement of 1.2 inches. 
From the load- settlement diagram it appears that the pile 
capacity had been reached, based on the usual criteria. 



Input data: 



W(1) 


9,300 


W(2) 


700 


W(3) - W120 


6hQ 


V 


12.9 


K(l) 


2,715,000 


rH 

H 

1 

CM 
' — ' 


7,270,000 


e l 


0.50 


e 2 


1.0 


T 


0.00025 



Distribution of resistance: 

This pile was equipped with strain rods in order that the 
distribution of load on the pile could be determined during 
the load test. The high value of set (low value of blows 
per inch) for this pile tends to indicate the possibility 
of liquefaction in the silt layers. Assuming this to be 
the case, correlation is made on the basis of the load 
carried to the sand layer which was 170 tons. 



Results: See Figures (U5) and (U6). 

Correlation is obtained for the following values 
of Q and J : 

Q J J l 

.1 .15 .05 
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TABLE 0 



Case : 
Source : 
Pile : 

Embedment : 
Hammer : 



31 

Page 31, reference (16). 
65' 12BP53 
55. U feet. 

Vulcan #1 



Soil: 



10' silty sand 
lit 1 loam 

26' fine sand, silty sand, and clay 
10* fine to medium sand 
10* fine to medium sand and gravel 
then fine sand 



Final penetration: 1.63 blows per inch. 

Load test results: k$ tons 17 days after driving. 

Results of calculations by the- wave equation method: 

a. If Q = 0.10 

J = 0.10 
J'= 0.03 

correlation is obtained with a side friction of 2$. 

b. If Q = 0.10 

J = 0.20 
J«= 0.067 

correlation is obtained with a side friction of 17$. 

c. If Q = 0.20 

J = 0.10 
J'= 0.03 

correlation is obtained with a side friction of 2%. 

d. If Q = 0.20 

J = 0.20 
J'= 0.067 

correlation is obtained with a side friction of 32$. 



Comments: 



Hiley Formula gives ultimate resistance of 52 tons. 
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TABLE DD 



Case: 
Source : 
Pile : 

Embedment : 

Haraner: 

Soil: 



3U 

Page 3U, reference (16). 

12BP53 

U7 feet. 

Vulcan #1 

35* silty clay 
then sand 



Final penetration: U.57 blows per inch. 

Load test results: 100 tons, but ultimate resistance not reached. 

Results of calculations by the wave equation method: 

a. If Q - 0.10 

J = 0.10 

J« = 0.03 

correlation cannot be obtained for these values for a 
failure load of 100 tons; however, these are the results: 

122 tons for no side friction 
1U6 tons for 25$ side friction 
185 tons for 50$ side friction 

b. If Q » 0.10 

J = 0.20 

J 1 = 0.067 

correlation is obtained with no side friction. 

c. If Q = 0.20 

J = 0.10 

J' = 0.03 

correlation is obtained with no side friction. 

d. If Q = 0.20 

J = 0.20 

J» » 0.067 

correlation is obtained with a side friction of 25$. 

Comments: It is not known just what the ultimate resistance would have 

been if the load test had been carried all the way to failure. 



The driven length of the pile was not given; 57' was assumed. 
Hiley Formula gives ultimate resistance of 96 tons. 
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TABLE P 



Case : 35 

Source; Page 35, reference ( 16 ). 

Pile ; 12BP53 

Embedment; $h.h feet. 



Pile 



Hamr.ier ; Vulcan #1 



Soil: 



15' sandy, silty, clay and loam 

10* sand 

10' sandy loam 

2k 1 sand 

then sandy loam 



Final penetration: U .57 blows per inch 



Load test results: 100 tons, but ultimate resistance not reached. 

Results of calculations by the wave equation method: 



a. If Q = 0.10 

J = 0.10 

J' = 0.03 

correlation cannot be obtained for these values for a failure 
load of 100 tons; however, these are the results: 



b. If Q = 0.10 

J = 0.20 

J« - 0.067 

correlation is obtained with a side friction of 8# for 
a 100- ton ultimate resistance. Other results are: 



c. If Q = 0.20 

J = 0.10 

J* = 0.03 

correlation is obtained with a side friction of 2% for 
a 100- ton ultimate resistance. Other results are: 



113 tons for no side friction 
ll*l tons for 25# side friction 
182 tons for 50 # side friction 



llU tons for 25# side friction 
151 tons for 50# side friction 



122 tons for a 25# side friction 
163 tons for a 50 # side friction 
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Case 35 (continued): 

d. If Q = 0.20 

J = 0.20 

J« = 0.067 

correlation is obtained with a side friction of 26 % for a 
100-ton ultimate resistance. Other results are: 

103 tons at 25 % side friction 
129 tons at 50$ side friction 

Comments: It is not known just what the ultimate resistance would 

have been if the load test had been carried all the way 
to failure. 

The driven length of the pile was not given; 65' was 
assumed. 

Hiley Formula gives ultimate resistance of 92 tons. 
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TABLE Q 



Case: 36 

Source: Page 36 , reference (16) 

Pile : 12BP53 

Embedment: U9.8 feet. 

Hammer : Vulcan #1 



Soil: 



15 1 sandy, silty, clay and loam 
then sand 



Final penetration: Iul5 blows per inch. 

Load test results: 100 tons, but ultimate resistance not reached. 

Results of calculations by the wave equation method: 



a. If Q = 0.10 

J = 0.10 

J 1 = 0.03 

correlation cannot be obtained for these values for a 
failure load of 100 tons. Other results are: 

108 tons for no side friction 
135 tons for 2 $% side friction 
178 tons for 50 $ side friction 

b. If Q « 0.10 

J = 0.20 

J' = 0.067 

correlation is obtained with a side friction of 7 % for 
a 100-ton ultimate resistance. Other results are: 

115 tons for 2 % side friction 
1U3 tons for side friction 

c. If Q = 0.20 

J = 0.10 

J* = 0.03 

correlation is obtained with a side friction of 1 % for 
a 100-ton ultimate resistance. Other values are: 



115 tons for 2 % side friction 
152 tons for side friction 
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Case 36 (continued): 

d. If Q = 0.20 

J = 0.20 

J* = 0.06? 

correlation is obtained with a side friction of 32 % fer 
a 100- ton ultimate resistance. Other values are: 

91 tons at 2 $% side friction 
121 tons at 50 % side friction 

Comments: It is not known just what the ultimate resistance would 

have been if the load test had been carried all the way to 
failure . 

The driven length of the pile was not given; 60' was 
assumed. 

Hiley Formula gives ultimate resistance of 92 tons. 
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TABLE R 


Case: 


67 


Source : 


Page 6?, reference (16) 


Pile: 


11U» 10HU2 


Embedment: 


95 feet 


Hammer: 


Steam, double acting, 950 pounds with stroke of 16 inches. 


Soil: 


Silt and fine sand, pile hit hard strata 



Final penetration: 32 blows per inch 

Load test results: 60 tons 

Results of Calculations by the wave equation: 



a. 


If Q - 0.10 
J = 0.15 
J* = 0.05 

correlation cannot be obtained; however, 100 tons was the 
calculated value for no side friction. 


b. 


If Q = 0.10 
J = 0.20 

J‘ » 0.067 

correlation cannot be obtained; however, 88 tons was the 
calculated value for no side friction. 


c. 


If Q = 0.20 
J* = 0.05 
J = 0.15 

correlation is obtained with a side friction of 1$%. 


d. 


If Q = 0.20 

J = 0.20 

J« = 0.067 

correlation is obtained with a side friction of U0$>. 


Comments : 


Time between driving and load testing is unknown. 
Hiley Formula gives an ultimate resistance of 25 tons. 
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TABLE S 



Case : 

Source: 

Pile: 

Embedment: 
Hammer : 
Soil: 



71 

Page 71, reference (16) 

20* 10HU2 
19 feet 

Gravity, 3,150-pound ram and 5-foot drop 
Sand and Gravel 



Final penetration: 2 blows per inch 

Load test results : 6i+ tons 

Results of Calculations by the Wave Equation: 

a. If Q » 0.10 

J - 0.15 

J‘ = 0.05 

correlation is obtained at a side friction of both 75$ 
and 100$. 

b. If Q = 0.10 

J = 0.20 

J' = 0.067 

correlation is obtained for a side friction of 100$. 

c. If Q = 0.20 

J = 0.15 

J* = 0.05 

correlation is obtained for a side friction of 100$. 

d. If Q = 0.20 

J = 0.20 

J* = 0.067 

correlation cannot be obtained; however, the ultimate 
resistance is calculated to be 60 tons for 100$ friction, 
or U0 tons for the end bearing case. 

Comments: Time between driving and load testing is unknown. 

Hiley equation gives an ultimate resistance of 80 tons. 
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TABLE T 



Case : 
Source : 
Piles 

Embedment : 

Hammer: 

Soil: 



75 

Page 75, reference (16) 

30 « 10Hl*2 
23 feet 

Drop — 3000-pound ram, 5-foot drop 
Fine Sand 



Final Penetration: 5 blows per 2 inches 
Load Test Results: 90 tons 

Results of Calculations by the Wave Equation: 

a. If Q = 0.10 

J = 0.15 

J« = 0.05 

correlation is obtained for 10C$ side friction. 

b. If Q ; - 0.10 

J = 0.20 

J« = 0.067 

correlation cannot be obtained for a failure load of 
90 tons; however, these are the results: 

55 tons for no side friction 
6k tons for 50$ side friction 
70 tons for 10G$ side friction 

c. If Q = 0.20 

J = 0.15 

J* = 0.05 

correlation cannot be obtained for a failure load of 
90 tons; however, these are the results: 

5U tons for no side friction 
60 tons for 50$ side friction 
60 tons for 1005& side friction 

d. If Q = 0.20 

J = 0.20 

Jt = 0.067 

correlation cannot be obtained for a failure load of 
90 tons; however, these are the results: 






ti 
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TABLE T (cont*d.) 

U8 tons for no side friction 
55 tons for 5Q& side friction 
60 tons for 10Q& side friction 

Comments: Time between driving and load testing unknown. 

Hiley formula gives Ru of 80 tons. 
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TABLE U 



Case : 
Source : 
Pile : 

Embedment : 

Hammer: 

Soil: 



76 

Page 76, reference (16) 
kO' 10Hi|2 
3k feet 

Drop — 3000-pound ram, 5-foot drop 
Fine Sand 



Final Penetration: 2 blows per inch 

Load Test Results: 75 tons on second loading, 50 tons first loading 

Results of Calculations by the Wave Equation: 



a. If Q = 0.10 

J = 0.15 

J* = 0.05 

correlation is obtained with side friction of 25$. 

b. If Q = 0.10 

J = 0.20 

J« = 0.067 

correlation is obtained for a side friction of 35$ (the 
difference in this case between 25 $ and 50$ friction 
is only 1 ton). 

c. If Q = 0.20 

J » 0.15 

J* = 0.05 

correlation is obtained with a side friction of 75$. 



d. If Q - 0.20 
J = 0.20 

J 1 = 0.067 

correlation is obtained with a side friction of 90$. 
Comments: Time between driving and load testing unknown, 

Hiley formula gives an ultimate resistance of 60 tons. 
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TABLE V 



Case: 
Source : 
Pile: 



89 

Reference (12 ) 

U 5 * U4BP89 



Embedment: I4J4. * 

Hammer: Vulcan 0 

Soil: Same as case # (l) 

Final Penetration: 30 blows per inch 

Results of Load Test: 300 tons 

Results of Calculations by the Wave Equation Method: 

Correlation is obtained with 25 % friction and 



Q = 0.20 
J = 0.20 
J* = 0.067 



Comment: This pile was tested 5 days after driving. The 300-ton 

load appears to be close to the ultimate resistance of 
the ground. Load increments were 75 tons, l£0 tons, 

225 tons, 225 tons (repeated), and 300 tons. 

Hiley formula gives ultimate resistance of 230 tons. 
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TABLE ¥ 



Case : 


113 


Source : 


Page 113, reference (16) 


Pile: 


UO* precast c oncrete weighing 2.7 tons 


Embedment: 


31.5' 


Hammer: 


Vulcan #1 



Driving Cap: 0.525 tons 



Soil: 


18' coarse sand 
26 * fine, yellow sand 



Final Penetration: 10 blows per 2.05 inches 

Load Test Results: 100 tons 

Results of Calculations by the Wave Equation: 



a. 


If Q = 0.10 
j = 0.15 
J' = o.o5 

correlation is obtained for 5C$ friction. 


b. 


If Q - 0.10 
J - 0.20 

J' = 0.067 

correlation is obtained with a side friction of 85$. 


c. 


If Q - 0.20 
J = 0.15 
J» = 0.05 

correlation is obtained with a side friction of 75$. 


d. 


If Q = 0.20 
J = 0.20 

J' = 0.067 

correlation is obtained with a side friction of 95$. 


Comments : 


Time between driving and load testing is unknown. 



Results using the Hiley equation are 122 tons. 
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TABLE X 



Case : 


117 


Source : 


Reference (12) 


Pile: 


b$ x HiBP117 


Embedment : 


feet 


Hammer: 


Vulcan 0 



Soil: 9* sand & gravel 

12* fine to medium sand 
15’ sand and gravel 
6* sandy silt 
2* weathered, silty shale 
then hard, silty shale 

Final Penetration: 22 blows per inch 

Results of Load Test: 300 tons 

Results of Calculations by the Wave Equation Method: 

Correlation is obtained with 25# friction and 

Q = 0.20 

J = 0.20 

J* = 0.06? 

t 

Comments: This pile was test loaded 2 days after driving. The 

300-ton load appears to be close to the ultimate 
resistance of the ground. Load increments were 75 tons, 
150 tons, 225 tons, 225 tons (repeated), aid 300 tons. 

I 

Hiley formula gives ultimate resistance of luO tons. 

i 
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TABLE Y 



Case : lUU 



Source s 
Pile : 

Embedment : 
Hammer: 



Page lUii, reference (16) 

120 » 12BP53 

81 . 5 * 

Vulcan #1 



Soil: 



5* clay loam 
11* sandy, silty, clay 
12' medium sand 
5* coarse sand 
39 * medium clay 
10' stiff clay 
10' very stiff clay 
10* medium clay 
11 1 hard clay with pebbles 
then medium clay 



Final Penetration: 8.0 blows per inch 



Load Test Results: 117 tons ultimate resistance 



Results of Calculations by the Wave Equation Method: 



a. If Q =■ 0.10 

J = 0.10 

J* = 0.03 

correlation cannot be obtained for a failure load of 117 
tons. At full point bearing (0$ side friction) the 
ultimate resistance is calculated to be 132 tons. 

b. If Q = 0.10 

J = 0.20 

J' = 0.067 

correlation is obtained with a side friction of $%. 

c. If Q = 0.20 

J = 0.10 

J» = 0.03 

correlation is obtained with a side friction of 0$. 

d. If Q = 0.20 

J « 0.20 

J‘ = 0.067 

correlation is obtained with a side friction of 20$. 



Comments: Hiley Formula gives ultimate resistance of 90 tons. 
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TABLE Z 



Case : 11*5 

Source: Page 1 1*5, reference (16) 

Pile: 106 « 12BP53 

Embedment : 105 * 

Hammer : Vulcan #1 

Soil: Same as case # 1)|)| 

Final Penetration: 10,3 blows per inch 

Load Test Results: 182 tons ultimate resistance 



Results of Calculations by the Wave Equation Method: 

a. If Q - 0.10 

J = 0.10 

J» = 0.03 

correlation is obtained with a side friction of 19$. 

b. If Q - 0.10 

J - 0.20 

J» - 0.067 

correlation is obtained with a side friction of 35$. 

c. If Q = 0.20 

J = 0.10 

J* = 0.03 

correlation is obtained with a side friction of 36$. 

d. If Q - 0.20 

J - 0.20 

J« « 0.067 

correlation is obtained with a side friction of bl%. 



Comment: Hiley Formula gives ultimate resistance of 86 tons. 
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TABLE AA 



Case: 


11*6 


Source : 


Page ll*6, reference (16) 


Pile: 


100 « 12BP53 


Embedment : 


97.1* 


Hammer: 


Vulcan #1 


Soil: 


Same as Case #!liU 



Final Penetration: U.75 blows per inch 

Load Test Results: 91 tons 

Results of Calculations by the Wave Equation Method: 

a. If Q » 0.10 
J - 0.10 





J» - 0.03 

correlation cannot be obtained for a failure load of 
91 tons. At full point bearing (0$ side friction) the 
ultimate resistance is calculated to be 110 tons. 


b. 


If Q - 0.10 
J = 0.20 

J* = 0.067 

correlation is obtained with a side friction of 2$. 


C 0 


If Q 8 0.20 

J - 0,10 

J* - 0.03 

correlation is obtained at 2 % side friction. 



d. If Q « 0.20 
J = 0.20 

J> = 0.067 

correlation is obtained with a side friction of 17$. 
Comment: Hiley Formula gives ultimate resistance of 72 tons. 
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TABLE BB 



Case: 


161 


Source : 


Page 161, reference (16). 


Pile: 


SO* 12 BPS 3 


Embedment : 

I 


32 feet 


I 

Hammer: 


Vulcan #1 


Soil: 


3* fill 

2 S* coarse sand 
18 1 medium sand 
lS* fine sand 
9 X medium clay 
20 * very stiff clay 
22 * medium clay 
then very stiff clay 



Final Penetration: It . 83 blows per inch. 

Load Test Results: 167 tons 

Results of Calculations by the Wave Equation Method: 



a. 


If Q = 0.10 
J = 0.10 

J* = 0.03 

correlation is obtained with a side friction of 3S$. 


b. 

| 


If Q * 0.10 
J = 0.20 
J* = 0.067 

correlation is obtained with a side friction of SS$. 


c. 


If Q = 0.20 
J = 0.10 

J* = 0.03 

correlation is obtained with a side friction of S3$. 


d. 


If Q - 0.20 
J = 0.20 

J* » 0.067 

correlation is obtained with a side friction of 70$. 


Comment : 


Hiley Formula gives ultimate resistance of 96 tons. 
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C. Discussion of Results . 

Each of the cases correlated by the first method are dis- 
cussed below: 

1. Case number 2001, table (F). This and the next two cases 
are 36-inch diameter prestressed concrete piles driven at Lake Maricaibo, 
Venezuela. The load carrying capacity of the pile tip was calculated 
using the equation, P = 9AC, and a value of cohesion (C) of 2560 p.s.f., 
giving a capacity of 163,000 pounds. This value is 1+2 per cent of the 
total load carried by the pile at failure, and the assumption that the 
load is carried one half in friction and one half by point bearing 
would be reasonable. If the load carried as friction is considered to 
be evenly distributed on the embedded length of the pile, each ten-foot 
section would develop a resistance equal to 1+7,500 pounds. If the 
distribution of the side friction resistance were assumed to vary 
directly as the depth, the distribution would be as shown in figure 
(26), which also compared both distributions with that estimated from 
the shear strength of the soil. In the calculated values the full 
value of the shear strength of the soft clay layers was assumed to be 
effective and a reduced value of 800 p.s.f. was used for the stiff 
clay layer between elevations -153 and -li+3 ft. This value is con- 
sistent with Nordlund’s estimate (ref. 23) and the findings of others 
that have investigated the friction values of piles in stiff clay 
(ref. 33, 38). Comparing the load distribution in figure (26), it 
appears that the rectangular distribution is as good a representation 
of the estimated actual load distribution as the triangular method 
would be. From figure (27) it is seen that values of Q of 0.1 and 
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J of 0 . 9 ) and Q of 0.2 and J of 0.3 produce nearly identical curves 
that give a resistance value slightly higher (200 tons vs. 190 tons) 
for the set of 20 blows per inch. 

Case 2002 . This is the same pile as case 2001 redriven to 
refusal with a heavier hammer. The load capacity was not increased 
a great deal by the redriving, going only frem 190 tons to 230 tons. 
Correlation was obtained only with a quake of 0.20 using a value of 
J of 0.80 as shown in figure (29). This is quite a departure from the 
value of J obtained in the previous case and would tend to indicate 
that for harder driving, quake increased. Correlation could not be 
obtained for a quake of 0.1 with values of damping as high as 1.2. 

Case 200U . The tip of the pile in this case penetrated into 
clay having a shear strength of 10,000 p.s.f. as determined from 
unconfined compression tests. Results of correlation taid to indicate 
a smaller value of quake than experiencedin the previous cases. 

Figure (38) shows the results from wave equation computations. 

Case 2003 . In this case correlation was obtained with values 
of Q and J of the same approximate values as in case 2001 as indicated 
in figures (3U) and (35). Figures (31) to (33) are included as a matter 
of interest, and show that correlation could be obtained with different 
distributions of resistance. 

Case 2005 . This pile is of the same type as the previous one 
but was driven into a somewhat softer clay having a shear strength of 
about 3000 p.s.f. Figure (39) is the basis of correlation, and indicates 
values of J in the range more-or-less expected. 



: ' 

:-‘ v 

• » 



■ 







.. • 



. i 

. • 









102 



Case 6006 . Tfte pile in this case was not loaded to failure, 
but, from analysis of the load test and driving records, the engineer 
supervising this test estimated the value of side friction and point 
bearing resistance developed under static loading. From this the 
ultimate load carrying capacity of the pile was computed to be a mini- 
mum of U98 tons. This value was used in correlating the load test 
with wave equation computations. The correlation for a value of Q of 
0.20 was obtained by extrapolation, but it is believed sufficiently 
accurate considering the general order of accuracy of the other vari- 
ables and data. Figure (1|0) is the load settlement curve obtained 
from the load test, and figures (1+1 ) and (1+2 ) show the basis for 
correlation. 

Cases 6007 and 6008 . These are two pipe piles driven through 
80 feet of plastic clay underlain by a thick sand stratum. One pile 
was driven about 5 feet into the sand and the other was stopped 5 feet 
above the sand layer. The object of the load tests was to determine 
the load carried by the point in sand. Correlation for the pile that 
extends into the sand (case 6008) was obtained in the lower end of the 
expected range of values for the point damping factor and quake, with 
values of 0.05 and 0.10 respectively. No attempt was made to account 
for the fact that the pile itself is largely embedded in clay. For 
the pile stopped above the sand layer, correlation was obtained for 
values of Q of 0.10, J of 0,15 and J‘ of 0.05. Since the load is 
carried entirely by friction, all of the damping is due to J‘, the 
side damping factor. These values are not consistent with the expected 
range of values for clay, and indicate that the value of quake is smaller 
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than expected, which would be consistent with the easy driving condi- 
tions, that the relationship between J and J' is in error, or that 
side damping does not develop when a steel pipe pile is driven through 
plastic clay. With respect to case 6008, pile point in sand, it is 
possible that the clay was somewhat sensitive and regained strength 
during the period before the load test was made. If the sensitivity 
of the clay were such that one half of its strength were lost due to 
remolding and then regained after a period of time, the ultimate re- 
sistance during driving would be reduced by about IjO tons. This would 
produce correlation at a value of Q of 0.10 and J of 0.15. 

Case 6011 . The high load capacity developed by this pile and 
the low resistance to penetration indicate that a quick condition -in 
the overlying silt layers occurred. If this were the case, the resistance 
to driving would be developed on the portion of the pile embedded in 
sand with no resistance occurring in the silt layers. Correlation is 
bbtained on this basis as shown in figures (U$) and (U6). 

The results of correlation attempts for these cases are sum- 
marized in figure (6?). The piles driven with the tip in sand are 
closely grouped and fall within a range of values for J of 0.05 to 0.2, 
and for values of Q of 0.0$ to 0.10. This result tends to confirm 
Smith’s recommended value for J of 0.1$, for piles driven in sand. 

Results of correlation for piles driven in clay are spread 
over a range of values for J of 0.15 to about l.U, and for values of 
Q of 0.1 and 0.2. Case 6007, which has a J value for correlation of 
0.1$, is a 10<$ friction pile, and the point damping factor (J) has no 
influence. This possibly indicates that the one-third relationship 
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Summary of Correlation Results — First Method 
Figure (6?) 



between J and J 1 is too gross an assumption in cases where all or most 
of the load is carried by friction. In general, the results in clay 
are consistent with values of J expected, based on dynamic tests made 
on soils in the laboratory. 



D. Summary of Correlation Results for Second Method of Approach 

For the cases studied, the following values of ground quake, 
point damping, and side damping seem to be appropriate for the soil 
types indicated. The percentages of side friction given are for steel 
H piles unless otherwise indicated. Furthermore, the principal case 
studies used in arriving at the values are shown: 

1. Coarse Sand: (Case Numbers 113 and l6l) 

Q = 0.10 
J = 0.15 
J« = 0.05 

Side Friction is 35$. 

2. Sand and Gravel Mixed: (Case Number 71) 

Q = 0.10 
J = 0.15 
J» = 0.05 

Side Firction is between 75$ and 10C$. 

3. Fine Sand: (Case numbers 75, 76, and 113) 

Q = 0.15 
J = 0.15 
J« = 0.05 

Side Friction is 100$ (use 50$ for straight-sided 

concrete pile). 

U. Sand layers combined with layers or strata of clay or 

loam or both, but sand layers predominate: (Case numbers 

31, 3h, 35 , 36 , liiU, U£, H*6, and l6l apply) 

Q = 0.20 

J = 0.20 

J* = 0.067 
Side Friction is 25$ 
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5* Silt and Fine Sand underlain by Hard Strata: (Case 
number 67) 

Q - 0.20 

J = 0.20 

J* = 0.067 

Side Friction is 1 jQ£. 

6. Sand and Gravel Underlain by Firm Strata such as Shale: 
(Case numbers 89 and 117) 

Q = 0.10 to 0.20 

J = 0.15 to 0.20 

J* = 0.05 to 0.067 

Side friction of 25$. 

Using the average values of Q, J, and J* and side friction 
obtained for the various types of soil conditions, the variations of 
calculated ultimate resistances from load test results for the specific 
cases will be investigated. 

Hiley equation results are also shown with this comparison in 
Table CC. 

It should be noted that the values of ultimate resistance shown 
to be calculated by the Hiley formula were obtained by using values of 
C^, C 2 , and C^ recommended by Chellis (ref. 5) for use in the absence 
of specific information. No attempt was made to determine what changes 
in these coefficients would be necessary to achieve closer agreement 
with load test results. Comparison of Hiley formula ultimate re- 
sistance calculations has therefore been restricted to those case 
studies where complete data were not provided. Whether someone having 
experience with the Hiley formula could have justified use of different 
values for these coefficients and thereby achieved closer agreement 
with load tests is not known. 
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TABLE CC 

COMPARISON OF WAVE EQUATION SOLUTION 
WITH HILEY FORMULA AND LOAD TESTS 



Ultimate Resistance 




Per Cent 


Error 


Computed by 
Wave Equation 


Computed by 
Hiley Formula 


Load 

Test 


Wave 

Equation 


Hiley 

Formula 


(tons) 


( tons ) 


(tons) 




T~ 


UO 


52 


U5 


-11.1 


-15.5 


100 


96 


100 


0.0 


- U.o 


103 


92 


100 


3.0 


- 8.0 


91 


92 


100 


- 9.0 


- -8.0 


60 


25 


60 


0.0 


-58. k 


6k 


80 


6h 


0.0 


25.0 


90 


80 


90 


0.0 


-11.1 


50 


60 


50 


0.0 


20.0 


300 


230 


300 


0.0 


*23.3 


300 


lUO 


300 


0.0 


-53.3 


100 


122 


100 


0.0 


22.0 


122 


90 


117 


U.3 


-23.0 


133 


86 


182 


-26.9 


-52.7 


97 


72 


91 


6.6 


-20.8 


103 


96 


167 


-38.3 


-U2.5 
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It must be remembered that the number of cases investigated 
here is small, and the values of point and side damping, ground quake, 
and side friction reported as giving good correlation may be replaced 
by more accurate values as more driving records and load tests are 
analyzed and as more basic information about the dynamic properties 
c£ soils becomes known. Furthermore, side resistance was taken as 
rectangular, or uniform distribution, along the side of the piles 
with about the top 10 feet of embedment discounted in each case. 
Different side distribution patterns would, of course, alter the values 
found above somewhat. 

On the other hand, the values reported are within the ranges 
that one might anticipate, and they do give good correlation as well 
as providing results generally more accurate than predictions from 
the Hiley formula. 
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Chapter VII 

RECOMMENDATIONS FOR FURTHER RESEARCH 

No attempt will^ be made here to make an exhaustive and detailed 
list of further research possibilities, but several will be discussed 
briefly. Probably as many have already occurred to the reader, and 
certainly many more will present themselves to anyone beginning work 
in this general area. The main factors reqi iring further evaluation 
are ground quake, point and side damping, and distribution of side 
frictional resistance. Specifically, the recommendations are as follows 

A. Construct an apparatus in which model piles may be loaded 
both statically and dynamically. If economically feasible, the model 
pile may be instrumented and recording pressure cells embedded in the 
soil. Also it would be desirable to be able to enclose the entire 
soil-pile system so that a variable confining pressure could be 
applied. With this set-up the model pile could be encased to avoid 
side friction effects in the initial experiments. The model pile 
could then be struck by a small drop hammer and the set determined. 

This procedure could be repeated for different sized hammers and 
strokes, and the results compared with computer results for the test 
set-ups. The ultimate resistance could be determined frctn test load- 
ing statically an identical pile in the same controlled soil. In this 
way values of quake and damping might be bracketed more closely for 
various soil-pile systems. 

B. Obtaining a picture of the interaction between pile and 
soil may prove useful. It is thought that this may be done for both 
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the static and dynamic cases. A photosensitive material could be 
used for a model pile which would be embedded in another photosensitive 
material representing the ground. The stress waves propagated by 
striking the model pile at its head by a small drop hammer could be 
captured in successive time intervals by using a high-speed camera. 

Even though such a study may give only qualitative results, it may 
assist in understanding how much of the ground is acting with a pile 
and the order of magnitude of stresses within this zone. Such informa- 
tion should prove helpful in evaluating both ground quake and the 
distribution of side frictional resistance. 

C. Study the problem from both an analytical and statistical 
viewpoint. Seek the cooperation of major construction corporations, 
government bodies, and others in obtaining as many good load tests 

as possible. For example, it is understood that the University of 
Michigan or the Bureau of Public Roads is on the verge of publishing 
a number of very carefully controlled tests. Run a wide range of 
curves for various pile and hammer combinations as was done herein for 
the ^3-pound bearing pile, and see if dynamic soil properties cainot 
be deduced from them which will give good correlation. If sufficient 
driving records and load tests become available, it may be advisable 
to write a computer program which will assist in matching up wave 
equation results with field data. 

D. A wide range of wave equation solutions may be run off 
for the common soil sampling rigs with a view toward correlating the 
information with the soil data and with pile driving and testing 
records. In this connection it should be noted that the pipe section 
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of the sampler should be broken down into one— foot increments and 
the time intervals of the calculations proportionately reduced in 
order to obtain stable calculations. It should also be kept in mind 
that since lateral support of the long pipe section is lacking in 
this case, some energy will be lost in transverse vibrations. How 
seriously this will affect the results is not known. 

E. The computer program included herein may be modified 
so that it can also handle the case of a long hammer, which should 
be represented as several weights in the mathematical model, and the 
case where the ram strikes the pile other than at its top. A, E. L. 

Smith illustrates how the mathematical model for these cases may be 
represented (ref. 28). 

F. More detailed study should be made of pile tip velocities 
at the time of maximum tip force, as determined from the computer re- 
sults, for a wide range of pile types, lengths, and types of hammers 
in order to establish a more valid basis for relating results of rapid 
loading tests made in the laboratory with pile driving action. 

G. With respect to the computer program, investigation should 
be made into the factors that cause instability in the computations. 

In a few computer results for this work, instability was encountered 
unexpectedly. These instances involved light piles, driven with 
relatively heavy hammers, or low values of ultimate resistance relative 
to the nonnal capacity of the pile. Investigations should be made to 
determine additional criteria for use in making up the input data so 

as to avoid instability. 

As previously indicated, these are but a few of the possibilities 
for further workj however, it is hoped that these listed will be suffi- 
cient to stimulate the interest of potential investigators. 
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CONCLUSION 

The prediction of a pile*s ultimate static bearing capacity 
from its dynamic behavior during driving has been an elusive goal 
of civil engineers for many years. Hundreds of ingenious empirical 
dynamic formulas have been devised which may give good results under 
particular conditions of driving equipment, pile, and soil combina- 
tions; but the difficulty has been in selecting the proper formula 
to obtain reliable results. Furthermore, these dynamic formulas 
have limited ranges of applicability since each ignores important 
aspects of the problem. Because of these limitations, their validity 
as a general method of approaching the problem at hand is open to 
question. Civil engineers have recognized these shortcomings and 
have tended to treat the dynamic formulas as crude guides to the 
ultimate bearing capacity of a pile in granular soils, recognizing 
their complete inadequacy for piles driven in cohesive soils. 

The application of the wave equation to the pile driving 
problem, as envisioned here, represents an attempt to obtain a 
general method of analysis which will assist in predicting the 
ultimate bearing capacity of a pile from its observed behavior under 
the last hammer blow during driving subject to the limitations 
enumerated in paragraph E below. The wave equation method of analysis 
is also useful for studying stresses occurring during driving and 
for selecting appropriate pile driving equipment to meet economically 
specific field conditions. The major contribution by A. E. L. Smith, 
a pioneer in this work, was in devising a mathematical representation 
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of the pile driving problem suitable for digital computer solution. 

He then developed such a computer program, and his published work 
which followed emphasized the stress determination problem. He did, 
however, suggest that the method could be used in relating the static 
tearing capacity of a pile to its dynamic behavior, and he offered 
values of damping and ground quake which he believed to be accurate 
enough to make use of the method practical. 

To examine the possibility of correlating pile driving 
records with load tests by this method of analysis has been the 
principal objective of this thesis. Before making detailed case 
studies, it was necessary to examine the validity of the method 
itself, to review the soil engineering aspects of the problem, and 
to develop the necessary computer programs. These findings, together 
with the results from 2U case studies, are as follows: 

A. Validity of the mathematical model used . The model de- 
vised by A. E. L. Smith has been accepted as correct by those com- 
menting on his work, but it had not been proven to be equivalent 

to the wave equation. The formal proof of this equivalency has been 
provided in Chapter III. 

B. Resistance to driving . It is shown in Chapter V that 

the expression used for the resistance of a pile to driving corresponds 
approximately to published experimental findings of others on the re- 
sistance developed in a soil due to high rates of loading. 

C. Computer programs . Although A E. L. Smith developed 
computer programs for the wave equation method of analysis, they were 
not available as he had not published them. Computer programs are 
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published herewith which were developed independently. Their validity 
is shown by comparing with both a manual solution for a sample problem 
which utilizes all parts of the program and with results published by 
A. E. L. Smith. Complete instructions for usirg these programs are 
provided in Appendices A, B, and C so that someone unfamiliar with 
computer work may utilize them. 

The computer program developed for the Hiley formula is also 
subject to the limitations inherent in an empirical dynamic pile 
driving formula. 

D. Ground damping, quake, and side friction . From published 
experimental work by others, approximate values of damping to be used 
with the wave equation method of analysis for some soils have been 
derived in Chapter V. Starting with a range of quake reported by 
others, a range of values of damping were found which caused the 
solution of the wave equation method to correlate with load tests. 

This range of values of damping was compatible with the values derived 
from the published experimental work. Using the values of quake and 
damping thus established, values of uniformly distributed side fric- 
tional resistance were determined for various soil and pile combinations. 

E. Usefulness of the wave equation method for predicting 
static bearing capacity . Attempts to correlate 2h pile driving records 
with load tests have shown surprising results. It appears that the 
wave equation method of pile driving analysis may bee one an accurate 
method which may be used together with other factors in predicting the 
ultimate static bearing capacity of a pile, as indicated by the results 



summarized in figure (67) and table (CC). Among these other factors 
which must be evaluated are the relaxation or set-up of the soil with 
time after driving, the group effect of piles, development of negative 
friction, consolidation of a clay layer beneath a pile tip, long-term 
changes in water table, corrosion or deterioration of the pile from 
other causes, and any other special problems which may be associated 
with a particular construction project. 

Furthermore, the problem has been treated as one-dimensional; 
i.e., lateral vibrations have not been included in the analysis. This 
should not, however, be a serious limitation on the method for most 
pile driving situations since a large proportion of the pile will be 
embedded as final penetration is reached. Hysteresis losses in the 
pile have not been included in the analysis, but their neglect is 
considered justified at this stage of development of the wave equation 
method. Since good correlation was possible using values derived 
from the published dynamic soil tests, it appears that the wave equation 
method is well suited for analysis of any size or type of pile being 
driven by any hammer in not only granular soils, but also cohesive soils. 

F. The dynamic interaction between pile and soil is incom- 
pletely understood and offers a fruitful area for further research. Even 
a qualitative picture of stress transfer should prove helpful. Addi- 
tional experimental work to isolate the parameters affecting ground 
quake, point damping, side damping, and distribution of point and side 
resistance to driving is needed as well as the evaluation of these 



factors. 
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G. Correlation with many driving records and load tests 
should be made before taking the values reported herein as correct. 
Future investigators may use freely the computer programs and plotted 
results of calculations already made, and presented here, to facilitate 



such work. 
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APPENDIX A 

DETAILS OF THE "BASIC COMPUTER PROGRAM" 

WHICH USES THE WAVE EQUATION 

General : This program is written in Fortran language for the 

IBM 7090 computer, but it can be used equally well with the IBM 709 or 
IBM 709I4. computer. It has as its objective the detailed analysis of 
stresses, velocities, accelerations, displacements, and permanent set 
of the pile as these factors vary with time. For details of the theory 
and the mathematical model, refer to chapter III. 

Termination : Analysis is terminated automatically when the 

velocities of all the blocks of the mathematical model are simul- 
taneously negative or equal to zero. Negative velocities are taken 
as being in the direction from the bottom of the pile toward its top. 
Termination is also scheduled when the permanent set of the pile be- 
comes constant. Effects beyond either of these times are considered 
to be of secondary importance for purposes of this study. Generally 
these conditions will be met before calculations have been made for 
300 time intervals, and in most cases much sooner. The time interval 
is usually taken as 0.00025 or 0.00033 seconds because of the high 
velocity of stress wave travel in materials used for piles; therefore, 
this program investigates what happens during the first 0.1 second 
or less after the hammer hits the pile. In the usual case the pile 
will achieve its permanent set well within this time. 

The program will also stop automatically if the velocity of 
either the pile cap or pile tip exceeds twice the velocity the ram 
had at the moment of impact. This feature is designed to help detect 
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instability (see Chapter III, C) in the calculations, but should not 
be relied upon entirely. The surest way of detecting instability is 
by plotting the displacements and velocities of each block in the 
model against time and then examining the curves so obtained for very 
sharp peaks or discontinuities. Should such irregularities be present, 
the time interval may be reduced and the program rerun and rechecked 
for stability. The importance of this check cannot be over-emphasized. 

Should the program run the full 300 time increments, no 
indication is given in the output as to why the program stopped. In 
such a case the results should be examined especially carefully, for 
in all probability either instability was present or the program 
stopped before the final permanent set was reached. In such a case 
the variable DP may be redimensioned to a larger value (second state- 
ment above statement number 207 in the beginning of the program), the 
statement immediately after statement number $00 should be changed to 
read DO 501 I = 1,X where X will be the new dimension of variable 
DP, and the statement immediately after statement number 302 should 
be similarly changed to read DO 101 N - 3,X. In all other cases of 
program termination the reason for stopping will be printed out auto- 
matically in plain language; for example, •'VELOCITY OF PILE CAP EXCEEDED 
TWICE THE RAM VELOCITY." 

Preparation of Input Data : All data is punched on standard 

data cards in the exact order indicated below and then placed immedi- 
ately after the * DATA card at the end of the source program. 
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CARD ORDER 


DESCRIPTION 


1 


Case number— any number up to 10 places may be 
assigned, and it should be punched within the 
first 10 places on the card. Blanks will be read 
as zeros; use no decimal points. 


2 


If there is another case with more data following 
this set of data, punch +1 in spaces 9 and 10 on 
this card; if this is the last case for this run, 
punch -1 in places 9 and 10. 


3 


In the first U places of this card punch the num- 
ber of weights in the mathematical model. Within 
spaces 5 through 8 punch the number of weights 
in the mathematical model minus 1. Use no decimals 


h 


Punch within the first 10 spaces the decimal form 
of the time interval to be used in seconds. 


The next M cards 


M is the number of weights in the mathematical 
model. One card is to be used for each weight. 
Punch the value in pounds of each weight in the 
first 10 spaces of each card. The first card is 
for weight 1 (the ram) and then consecutively 
down to the pile tip. Use a decimal point and 
the number can be located anywhere within the 
first 10 spaces. 


The next (M-l) cards 


Punch within the first 10 spaces of each of the 
cards the spring constants beginning with the cap- 
block and extending to the top of the pile tip 
section of the model. The values are to be in 
pounds per inch. Use a decimal point and the num- 
ber can be located anywhere within the first 10 
spaces. 


The next M cards 


Punch within the first 10 spaces of each card the 
soil spring coefficients beginning with block 1, 
At least the first 2 or 3 will be zero normally. 
Use a decimal point and the number in pounds per 
inch can be placed anywhere within the first 10 
spaces. 


The next card 


Punch with a decimal point anywhere within the 
first 10 spaces on this card the value of the 
capblock coefficient of restitution. 


The next card 


Punch with a decimal point anywhere within the 
first 10 spaces on this card the value of the 
coefficient of restitution of the head packing 
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The next card 

The next card 
The next card 

The next card 

The next card 

The next card 



in the case of a concrete pile. For all other 
cases without head packing punch 1.00. 

Punch with a decimal point anywhere within the 
first 10 spaces the velocity of the ram at the 
instant of impact in feet per second. 

Punch with a decimal point anywhere within the 
first 10 spaces the value of ground quake in inches. 

Punch with a decimal point anywhere within the 
first 10 spaces the value of point damping to be 
used. 

Punch with a decimal point anywhere w ithin the 
first 10 spaces the value of side damping to be 
used. 

If the pile cap is allowed to transmit tension, 
enter a +1 in spaces 9 and 10; if not, enter a 
-1 there. Use no decimals. 

If side friction is present, enter a +1 in spaces 
9 and 10;. if not, enter a -1 there. Use no 
decimals . 



Running the Program ; The data deck of cards just described 
is placed immediately after the source program deck and is ready to be 



turned over to a machine operator for processing. 



Results : The computer will clearly write out in the output 

the input data; the number of each time interval; and identify the dis- 



placements, velocities, permanent sets, forces, and ground resistances 
acting on each block. The results are printed in decimal form with a 
multiplier notation. For example, 0.2207 E06 means 0.2207 x 10^ or 
220,700. 



The Program : The complete program is presented so that any 

who wish to use it may do so freely. Each line of the program is punched 
on a separate card and must be assembled in the exact order shown. Any- 
one preparing data cards from this program should verify painstakingly 
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that no mistakes have crept in during punching, for the slightest dis- 
crepancy can lead to totally erroneous results. For that reason those 
preparing cards from the program given here are cautioned to try it 
out with test cases before relying on the results. 

The notation used in the program is the same as was used in the 
derivation of the basic equations except as dictated by Fortran require- 
ments. The more important variations in terminology are: 

S instead of K for spring constants 

SP instead of K* for soil spring coefficients 

RES1 instead of e-^ for coefficient of restitution for the 

capblock 

RES2 instead of e ^ for coefficient of restitution for 

head packing 

Z instead of J for point damping factor 

ZP instead of J’ for side damping factor 

Addition all y, the distinction between the present time interval, the 
last time interval, and the second previous time interval is made by 
use of the letters L, M, and S respectively. For example, in the case 
of velocities, 



VL is used instead of V 

VM is used instead of v 

VS is used instead of v* 

Other variables introduced in the program are defined as they occur 
within the program or have been explained under preparation of data 



above. 
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*ID REESE C2165 CIVIL ENGR G WAVE EQ 5 $Q 

% XEQ 

* LABEL 

C WAVE EQUATION 

DIMENSION DS(IOO), DM(IOO), DL(IOO), VS(lOO), VM(IOO), VL(IOO) 
DIMENSION DPSM(IOO), DPSL(lOO) 

DIMENSION S(IOO), W(lOO), DP(300), R(lOO), SP(lOO) ,CS(100) 
DIMENSION CM(IOO), CL(lOO),FS(lOO), FM(IOO), FL(lOO),DPSS(lOO) 
207 DO 500 I = 1,100 

DS(I) = 0.0 
DM(I) = 0.0 
DL(I) =0.0 
VS(I) = 0.0 
VM(l) = 0.0 
VL(I) = 0.0 
DPSM(I) =0.0 
DPSL(I) =0.0 
S(I) = 0.0 
W(I) = 0.0 
R(I) = 0.0 
SP(I) = 0.0 
CS(I) = 0.0 
CM(I) =0.0 
CL(I) = 0.0 
FS(l) =0.0 
FM(I) =0.0 
FL(I) =0.0 

500 DPSS(I) = 0.0 

DO 501 I = 1,300 

501 DP (I ) =0.0 

READ INPUT TAPE 5,27,NCASE,ITURN 

27 FORMAT (110/110) 

0 READ INPUT TAPE 5, 70,M,TM1,T, (W(I),I=1,M), (S(l),I=l,IMl) 

1 ,(SP(I),I-1,M), RES1, RES2, V,Q,Z,ZP,CAP,RUB 
70 F0RI4AT (2IU/(F10.U)) 

iffilTE OUTPUT TAPE 6, 30, NCASE 
30 FORMAT (12H CASE NUMBER 15) 

WRITE OUTPUT TAPE 6,28, (I,W(I),S(I ),SP(I),I=1,M) 

28 FORMAT (5H DATA//58H M W S 

1 SP/(IU,3F20.2)) 

• WRITE OUTPUT TAPE 6,29,RES1,RES2,ZP,Z,Q,V 

29 F0RI4AT (36H COEF. OF REST. OF CAPBLOCK (HES1) = FU.2/ 

1 36 H CCEF. OF REST. OF PILECAP (RES2) = FU.2/ 



2 36H SIDE DAMPING FACTOR (ZP) = FU.2/ 

3 36H POINT DAMPING FACTOR (Z) = FU.2/ 

U 36H GROUND QUAKE (Q) = FU.2/ 

5 36H INITIAL RAM VELOCITY (V) = F6.2) 

WRITE OUTPUT TAPE 6, 31, T, CAP, RUB 
31 FORMAT ( 36 H TIME INTERVAL = F8.6/ 
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1 36H CAP = F6.2/ 

2 36H RUB = F6.2) 

WRITE OUTPUT TAPE 6,1 

1 FORMAT (8ijH M DL VL DP 

1 FL R//) 

DS(l) = V*12.*T 

VS(1) = V + ( -DS ( 1 )-*S ( 1 ) )-::-Tx-32 . 2/W( 1 ) 

VS(2 ) = (DS ( 1 )*S ( 1 ) )-*T-::-32 . 2/W( 2 ) 

DM(1) = DS(l) + VS(l)*12.-*T 
DM(2 ) = VS(2)-;:-12.-x-T 
CM(2) = DM(2 ) 

VM(l) = VS(l) -((DM(l)-DM(2))*S(l))*T*32.2/W(l) 

VM(2) = VS(2)+((DM(l)-DM(2))-*S(l)-(DM(2)-DM(3))*S(2))*T*32.2/W(2) 
VM(3) = VS(3) +((DM(2)-DM(3))*SJ(2)-R3)^32.2M3) 

CM(1) = DM(1) - DM(2) 

CHECK » 2.-X- V 

WRITE OUTPUT TAPE 6,301,DS(1), VS(1) 

301 FORMAT (8H N « 1,6H 1 2E15.U//) 

WRITE OUTPUT TAPE 6,302,DM(l),VM(l),DM(2), VM(2) 

302 FORMAT (OH N » 2,6H 1 2El5 . U/ll|H 2 2E15.1*//) 

DO 101 N = 3, 300 

WRITE OUTPUT TAPE 6,303,N 

303 FORMAT (5H = 13) 

MLESS1 = M-l 

DO 130 I = 1,M 

130 DL(I) = DM(I) + VM(I)*12.#T 
DO 131 I = 1,MLESS1 

131 CL(I) = DL(I) - DL(I+1) 

13 IF(DL(M)) 12,12,11* 

12 DE = 0.0 
GO TO 2i* 

1U IF(DL(M) - Q) 16,16,18 
16 DE = 0.0 
GO TO 21* 

18 VALUE » DL(M) - Q 

IF (VALUE-DE) 22,22,20 
20 DE = VALUE 
GO TO 2k 
22 DE = DE 

2k DP(N) = DE 

R(M) - (DL(M)-DP(N) )-x-SP(M)-x-(l. +Z*VM(M) ) 

IF (DP(N)) 25,25,26 
26 IF (DP(N) - DP(N-l) ) 197,197,25 
25 CONTINUE 
3 VAL = S(l)*CL(l) 

IF(CL(l)-CM(l)) 5,5 , h 
1* FL(1) = VAL 

GO TO 33 

5 CL(l) = CM(l) 

FL( 1 ) -VAL/(RES1**2 ) - (l . /(HESS**2 ) -1 . )*3(l.)*CL( 1 ) 

33 VAL2 - S(2)*CL(2) 
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IF(CL(2)-CM(2)) 35,35,3U 
3k PL(2) = VAL2 
GO TO 37 

35 CL(2) = CM(2 ) 

FL(2) = VAL2/(RES2-::-x-2 )-(l./(RES2-«-«2 )-l. )-x-S(2 )*CL(2 ) 

37 IF (CAP) 38,36,36 

38 IF(FL(2) ) 39,36,36 

39 FL(2) = 0.0 

36 CONTINUE 

7 IF(FL(l)) 8,133,133 

8 FL(1) = 0.0 

133 DO 132 I = 3.MLESS1 
132 FL(I) = CL(I)*S(I) 
kZ IF(RUB)U9,U9,Ui 
10; DO U8 I = 3,MLESS1 
DPSL(I) = DPSM(I) 

CHANGE = DL(I) - Q 

SUM + DL(I ) + Q 

IF ( DPSL( I ) -CHANGE )H5 , U6, 46 

45 DPSL(I) = CHANGE 

46 IF (DPSL( I )-SUM)U8, 48, 47 

47 DPSL(I) = SUM 

U8 R(I) = (DL(l) - DPSL(l))*SP(lMl. + ZP*VM(l)) 

49 DO 50 I = 1,M 

VL(I) = VM(I) + (FL(l-l) - FL(I ) - R(I))-*T*32.2/W(I) 

50 CONTINUE 

DO 100 I = 1,M 

703 WRITE OUTPUT TAPE 6,304, 1,DL(l),VL(l),DPSL(I),FL(l),R(l) 
304 FORMAT (11H I3,5El5.4) 

100 CONTINUE 

DO 75 K = 1,M 
X = 1. 

IF (VL(K)) 75,75,800 

75 CONTINUE 

GO TO 199 
800 CONTINUE 

IF(VL(2 ) - CHECK )99, 99,190 
99 IF ( VL (M ) -CHECK )102 , 102 , 19-U 

102 DO 98 K = 1,M 
DS(K) = DM(K) 

DM(K) = DL(K) 

DL(K) = 0.0 
VS(K) = VM(K) 

VM(K) = VL(K) 

VL(K) = 0.0 
CS(K) = CM(K) 

CM(K) = CL(K) 

CL(K) = 0.0 
FS(K) = FM(K) 

FM(K) = FL(K) 

FL(K) = 0.0 
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DPSS(K) = DPSM(K) 

DPSM(K) = DPSL(K) 

98 DPSL(K) =0.0 

WRITE OUTPUT TAPE 6,3(£,DP(N) 

305 FORMAT (litH E15.U) 

101 CONTINUE 
GO TO 196 

190 WRITE OUTPUT TAPE 6,191,N 

191 0 FORMAT (6I4H VELOCITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS 13) 

GO TO 196 

19U WRITE OUTPUT TAPE 6,195,N 

195 0 FORMAT ( 6 I 4 H VELOCITY OF PILE TIP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS 13) 

GO TO 196 

197 WRITE OUTPUT TAPE 6,198,N 

198 FORMAT (2oH DP BECAME CONSTANT AT N = 13) 

GO TO 196 

199 WRITE OUTPUT TAPE 6,200,N 

200 FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = IU ) 

196 IF(ITURN)205,20$,207 
205 CALL EXIT 

END 

* DATA 
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APPENDIX B 

DETAILS OF THE "VARY Ru PROGRAM" WHICH USES 1HE WAVE EQUATION 

General . This program is identical to the BASIC PROGRAM except 
that it provides the capability of varying certain input data in order 
to obtain wave equation solutions for a range of input variables. The 
program provides for a maximum of 500 time cycles, instead of 300, to 
allow for the use of smaller time increments and a greater number of 
pile segments. Otherwise, the termination features are the same as for 
the BASIC PROGRAM. 

Program Features . The program is designed with the following 
features that are different from or in addition to those of the basic 
program. 

1. The ultimate resistance (Ru) is varied over a preselected 
range of values in any increment desired. The set is computed for 
each value of Ru. The initial value of Ru, the increment of increase 
of Ru and the number of times the increment is applied are provided as 
input data. 

2. The effective length of embedment of the pile can be varied 
and is determined by input data. This is the length of the pile over 
which the resistance to driving acts. 

3. Five patterns of resistance distribution between point bear- 
ing and side friction can be applied for computation. These patterns 
vary from full point bearing to full side friction with intermediate 
distributions of 75, 50, and 25 per cent point bearing with the balance 
distributed as side friction. 
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U. The method of distributing side friction can be rectangu- 
lar or triangular as desired. Input data determine the selection of 
the method. 

5. The capability of recycling the entire program up to 
three times with the new values of any variable is provided. This 
feature was used to provide new values of Q and J, although other 
variables could be used. This capability is provided with a c ontrol 
statement, and the new values of the variables to be changed must be 
put into the program itself, a somewhat awkward arrangement. 

6. Any number of sets of data can be used in a single run, 
as provided for in the input data. This is also a feature of the 
basic program. 

Input Data . The input data is the same as for the basic 
program except for the additional data required. 

CARD ORDER DESCRIPTION 

2a After the second card of the basic program, 

areas of the pile head, center and tip in square 
inches are punched on one card in the first 
three ten-column fields. Values are punched with 
a decimal point. 

After the last card In the first four columns punch 1 to indicate 
of the basic program triangular distribution of resistance on the side 

or -1 to indicate rectangular distribution. 

Next card In the first 10 columns punch the initial value 

of Ru in pounds, and in the second 10-column 
field the increment of Ru in pounds. A decimal 
point must be used. These are variables RU1 
and ADD. 

Next card In the first four columns punch the number rep- 

resenting the number of times the increment of 
Ru is to be applied; this establishes the.maximum 
value of Ru. In the second four columns (5-8) 
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punch the number of blocks in the mathematical 
model to which ground resistance is not to be 
applied; this establishes the effective embedded 
length of the pile. These variables are ICYCLE 
and NOLOAD and are punched without a sign or 
decimal point. 

Next card A number of 1 through 5 is punched in each of the 

first three four-column fields for the variables 
LA, LB, and I£ respectively. These numbers pro- 
vide the values for the indexing parameters of 
a DO statement that determines the distribution 
of resistance between point and side. 

Next and last card A number with a value between 1 through U is 

punched in the first four columns. This number 
controls the number of times the program is to 
be repeated with new variables. No sign or deci- 
mal is used. This variable is called IRES. If 
the number 1 is used the program will cycle the 
full three times, and if U is used it will stop 
after all computations using the initial data 
are made. 

Output Data . In addition to listing all input data, the follow- 
ing output is provided: 

1. Value of Ru in pounds for the particular computation. 

2. Set in blows per inch. 

3. Compression and tension stresses at the head, center and 
tip of the pile, in pounds per square inch. 

U. The number of time cycles run. 

Use of the Program . This program is designed to provide the 
data necessary to plot Ultimate Resistance vs. Set curves for various 
values of Q, J,and J*, and for one or more of the five different dis- 
tributions of resistance between point and side. The values of Q, J, 
and J* are initially provided in the input data. If additional values 
are desired for computations, they are provided through the use of the 
"computed-go-to" staterrent following statement number 600. The 
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particular values must be placed in the program following statements 
700, 701, and 702. This feature could easily be expanded if desired. 

If it is desired to use all five distributions of resistance between 
point and side, the values of LA, LB, and LC must be 1, 5, and 1 re- 
spectively. When the index variable of the DO statement (five state- 
ments after 777) takes different values, the distribution of resistance 
is as indicated below: 

Value Distribution 

1 All point bearing 

2 75 per cent point bearing, 25 per cent friction 

3 50 per cent point bearing, 50 per cent friction 

h 25 per cent point bearing, 75 per cent friction 

5 All side friction 

By assigning proper values to LA, LB, and LC, the desired distributions 
are obtained in accordance with the operation of a DO statement as ex- 
plained in reference (21). 
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REESE C2165 CIVIL ENGR G WAVE EQ 2 $ 250 VARY RU 

XEQ 

LABEL 

WAVE EQUATION VARY RU 101 

DIMENSION DS(IOO), DM(IOO), DL(lOO), VS(lOO), VM(lOO), VL(lOO) 
DIMENSION DPSM(lOO), DPSL(lOO) 

DIMENSION S(100) W(100), DP(500) R(100), SP(100),CS(100) 
DIMENSION CM(lOO), CL(100),FS(100), FM(IOO), FL(lOO),DPSS(lOO) 
READ INPUT TAPE 5,27,NCASE,ITURN 
FORMAT (110/110) 

READ INPUT TAPE 5,71,ABEA1,AREA2,AREA3 
FORMAT( 3F10.lt) 

READ INPUT TAPE 5, 70 , M,IM1,T,(W(I),I»1,M),(S(I),I»1,IM1) 
,(SP(l),I-l,M), RES1, RES2, V,Q,Z,ZP,CAP,RUB 
FORMAT (2Ilt/(F10.lt)) 

READ INPUT TAPE 5,60lt,MUD 
FORMAT(lh) 

READ INPUT TAPE 5,72,RU1,ADD,ICYCLE,N0L0AD 
FORMAT (2F10.lt/2IlO 
READ INPUT TAPE 5, 73, LA, LB, LC, IRES 
FORMAT (3IU/IU) 

WRITE OUTPUT TAPE 6,30,NCASE 
FORMAT (12H CASE NUMBER 15) 

WRITE OUTPUT TAPE 6,28, (I,W(l),S(I),SP(I),I»l,M) 

FORMAT (5H DATA//58H M W S 

SP/(IU,3F20.2)) 

WRITE OUTPUT TAPE 6,29,RES1,RES2,ZP,Z,Q,V 

FORMAT (36H COEF. OF REST. OF CAPBLOCK (RESl) = Fit. 2/ 

36H COEF. OF REST. OF PILECAP (RES2) = Fit. 2/ 

36H SIDE DAMPING FACTOR (ZP) = Fit. 2/ 

36H POINT DAMPING FACTOR (Z) = Fit. 2/ 

36H GROUND QUAKE (Q) = Fit. 2/ 

36H INITIAL RAM VELOCITY (V) - F6.2) 

WRITE OUTPUT TAPE 6,6Q5,MUD 
FCRM/IT (6H MUD = 16) 

WRITE OUTPUT TAPE 6,31,T,CAP,RUB 

FORMAT (36H TIME INTERVAL = F8.6/ 

36 H CAP « F6.2/ 

36 H RUB = F6.2//) 

WRITE OUTPUT TAPE 6,608,RU1,ADD,ICYCLE,N0L0AD 

FORMAT (5H RU = E20.It,10H ADD = E20.lt/9H ICYCLE = 16, 

13H NOLOAD = 16) 

WRITE OUTPUT TAPE 6,609,LA,LB,LC,IRES 

FORMAT (5H LA = Iit,10H LB = IU,1CH LC = Ilt,10H IRES 

IU) 

CHECK = 2 .* V 
CONTINUE 
RU = RU1 

DO 600 IAM = 1, ICYCLE 
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ITEM = M - NOLOAD 
AITEM = ITEM 
DO 602 J = LA, LB, LC 
FMAX1 = 0.0 
FMAX2 =0.0 
FMAX3 =0.0 
FTEN1 =0.0 
FTEN2 =0.0 
FTEN3 =0.0 
DO 500 I = 1,100 
DS(I) = 0.0 
DM(I) =0.0 
DL(I) = 0.0 
VS(I) = 0.0 
VM(I) = 0.0 
VL(I) = 0.0 
DPSM(I) =0.0 
DPSL(I) = 0.0 
R(I) = 0.0 
SP(I) = 0.0 
CS(I) = 0.0 
CM(I) = 0.0 
CL(I) = 0.0 
FS(l ) = 0.0 
FM(I ) = 0.0 
FL(I) = 0.0 

500 DPSS(l) =0.0 
DO 501 I = 1,500 

501 DP (I) = 0.0 

DS(l) = V*12.*T 

VS(l) = V +(-DS(l)*S(l))*T*32.2/W(l) 

VS(2 ) = (DS(1)*S(1))*T*32.2/W(2) 

DM(l) = DS(l) + ?S(i)*12.*-T 
DM(2) = VS(2)#12.*T 
CM(2 ) = DM(2) 

VM(l) = VS(l) -((DM(l)-DM(2))*S(l))*T*02.2/W(l) 

VM( 2 ) = VS ( 2 )+( (DM( 1 )-DM( 2 ) )*S ( 1 ) - (DM( 2 ) -DM( 3 ) )*S ( 2 ) )*T*32 . 2/W( 2 ) 
VM(3) = VS ( 3 ) + ( ( DM ( 2 ) -DM ( 3 ) ) *-S ( 2 ) -R ( 3 ) ) #T*3 2 . 2/W ( 3 ) 

CM(l) = DM(1) - DM(2) 

AJ = J 

PART = (5.-AJ)/U. 

SP(M) = PART*RU/Q 
SIDE = RU - PART*RU 
LOAD = NOLOAD + 1 

if(mud)65o,65i,652 

650 DO 601 N = L0AD,IM1 
601 SP(N) = SIDE/(AITEM*Q) 

SP(M) = SP(M) + SIDE/ ( AI TEM-*Q ) 

GO TO 651 
DO 603 N = 3,11*11 
NT = N-2 
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NT1 = 2*NT - 1 
IM2 = M-2 
ANTI = NT1 
AIM2 = IM2 

603 SP(N) = SIDE*ANTl/( (AIM2**£)*Q ) 

SP(M) = SP(M) + SIDE-*( 2 .-*AIM2-1 . )/( (AIM2**2 )*Q ) 

651 CONTINUE 

DO 101 N = 3,500 
MLESS1 = M-l 
DO 130 I = IB 

.130 DL(I) = DM(I) + VM(I)*12.*T 
DO 131 I = 1 MLESS1 

131 CL(I) « DL(I) - DL(I+1) 

13 IF(DL(M) ) 12, 12, lii 

12 DE = 0.0 
GO TO 2h 

111 IF(DL(M) - Q) 16,16,18 
16 DE = 0.0 
GO TO 2h 

18 VALUE = DL(M) - Q 

IF (VALUE-DE) 22,22,20 
20 DE = VALUE 
GO TO 2k 
22 DE = DE 

2k DP(N) = DE 

R(M) = (DL(M)-DP(N ) )*SP(M)*(l. +Z*VM(M) ) 

IF (DP(N)) 25,25,26 
: 26 IF (DP(N) - DP(N-l) ) 197,197,25 
25 CONTINUE 
3 VAL = S(1)*CL(1) 

IF(CL(l)-CM(l) ) 5,5,u 
k FL(1) = VAL 
GO TO 33 

5 CL(1) = CM(1) 

FL(l ) =VAL/ (RES1#**2 )-(l./(RESl#*2)-l. )*S(l)-!fCL(l) 

33 VAL2 = S(2)*CL(2) 

IF(CL(2)-CM(2)) 35,39,3k 
3k FL(2) = VAL2 
GO TO 37 

35 CL(2) = CM(2) 

FL(2) = VAL2 / (RES2**2 ) - (1 . /(RES2**2 )-l . )*S ( 2 )*CL( 2 ) 

37 IF (CAP) 38,36,36 

38 IF(FL(2)) 39,36,36 

39 FL(2) = 0.0 

36 CONTINUE 

7 IF(FL(1)) 8,133,133 

8 FL(l) = 0.0 

133 DO 132 I = 3.MLESS1 

132 FL(I) = CL(I)*S(l) 

IF (FL(2 ) )213, 211, 210 

213 IF(FTEN1-FL(2 ) )211, 211, 2lU 
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2lU FTEN1 = FL(2 ) 

GO TO 21 1 

210 IF(FMAXl-FL(2))2l5,211,211 
205 FMAX1 = FL(2 ) 

211 CONTINUE 

MC = M/2 + 1 
IF(FL(MC))2l6,2l7,2l8 

216 IF ( FTEN 2 -FL ( MC ) ) 217 , 21 7 , 219 

219 FTEN2 = FL(MC) 

GO TO 217 

218 IF (FMAX2-FL(MC ) )220, 217, 217 

220 FMAX2 = FL(MC) 

217 CONTINUE 

IF (FL(IMl) )221, 222, 22U 

221 IF(FTEN3-FL(na))222,222,223 
223 FTEN3 = FL(IMl) 

GO TO 222 

22h IF(FMAX3-FL(IMl))225, 222,222 
225 FMAX3 = FL(IMI) 

222 CONTINUE 

k2 IF(RUB)U9,U9,H* 
hk DO hB I = 3,MLESS1 
DPSL(I) = DPSM(I) 

CHANGE = DL(I) - Q 
SUM = DL(I) + Q 
IF(DPSL(I )-change)U5,U6,U6 
U5 DPSL(I) = CHANGE 
h 6 IF(DPSL(I)-SUM)U8,ii8,U7 

U7 DPSL(I) = SUM 

U8 R(I) = (DL(I) - DPSL(I))-*SP(I)-:Kl. + ZP*VM(l)) 

U9 DO 50 I = 1,M 

VL(I) = VM(I) + (FL(I-l) - FL(I) - R(I))*T*32.2/W(I) 

50 CONTINUE 

DO 75 K = 1,M 
X - 1. 

IF (VL(K)) 75,75,800 

75 CONTINUE 

GO TO 199 
800 CONTINUE 

IF(VL(2) - CHECK)99,99,190 
99 IF ( VL (M )-CHECK ) 102 , 102 , 19U 

102 DO 98 K = 1,M 
DS(K) = DK(K) 

DM(K) = DL(K) 

DL(K) =0.0 
VS(K) = VM(K) 

VM(K) = VL(K) 

VL(K) = 0.0 
CS(K) = CM(K) 

CM(X) = CL(K) 

CL(K) = 0.0 
FS(K) = FM(K) 
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VM(K) = FL(K) 

FL(K) = 0.0 
DPSS(K) = IPSM(K) 

DPSM(K) = DPSL(K) 

98 DPSL(K) = 0.0 

101 CONTINUE 
GO TO 196 

190 WRITE OUTPUT TAPE 6,191,N 

191 0~ FORMAT ( 6 L|H VELOCITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS 13) 

GO TO 196 

19it WRITE OUTPUT TAPE 6,195,N 

195 0 FORMAT ( 6 I 4 H VELOCITY OF PILE TIP EXCEEDED MCE THE RAM VELOCITY W 
1 HEN N WAS 113) 

GO TO 196 

197 WRITE OUTPUT TAPE 6,198,N 

198 FORMAT (26H DP BECAME CONSTANT AT N = 13) 

GO TO 196 

199 WRITE OUTPUT TAPE 6,200,N 

200 FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = lit) 

196 CONTINUE 

WRITE OUTPUT TAPE 6,607,RU 

607 FORMAT (27H ULTIMATE RESISTANCE (RU) = E15.U) 

WRITE OUTPUT TAPE 6,330,Q,W(l), V 
330 FORMAT (ItH Q = E10.2,7H W(l) = E15.U, ItH V = El5.lt) 

WRITE OUTPUT TAPE 6,606, (SP(L),L=1,M) 

60 6 FORMAT (17H NEW VALUES OF SP/(8El5.it)) 

WRITE OUTPUT TAPE 6,305,DP(N) 

305 FORMAT (7H SET = ElO.It) 

BLOW = l./DP(N) 

WRITE OUTPUT TAPE 6,805,BL0W 
805 FORMAT ( 66 H 

1 BLOWS/EN= E10.3) 

FTEN1 = FTENl/AREAl 
FTEN2 = FTEN2/AREA2 
FTEN3 = FTEN3/AREA3 

WRITE OUTPUT TAPE 6,307,FTEN1,FTEN2,FTEN3 
307 0 FORMAT (3CH MAXIMUM TENSION AT HEAD = El5.lt/ 

1 30H MAXIMUM TENSION AT CENTER = El5.lt/ 

2 3CH MAXIMUM TENSION AT TIP = El5.lt) 

FMAX1 = FMAX1/AREA1 

FMAX2 = FMAX2/AREA2 
FMAX3 = FMAX 3 // 1 REA 3 

VJRITE OUTPUT TAPE 6,306,FMAXl,FMAX2,FMi'tX3 

306 0 FORMAT (26H MAXIMUM FORCE AT HEAD = El5.lt/ 

1 26H MAXIMUM FORCE AT CENTER = El5.lt/ 

2 26H MAXIMUM FORCE AT TIP = El5.1t//) 

602 CONTINUE 

RU = RU + ADD 
600 CONTINUE 
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GO TO (700, 701, 702, 703), IRES 

700 CONTINUE 
Z = 0.20 
ZP = 0.067 
IRES = 2 
GO TO 777 

701 CONTINUE 
Q = 0.20 
IRES = 3 
GO TO 777 

702 CONTINUE 
Z = 0.15 
ZP = 0.05 
IRES = U 
GO TO 777 

703 CONTINUE 
IF(ITURN)205,205,207 

205 CALL EXIT 
END 

* DATA 
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APPENDIX C 

DETAt LS OF THE "RESEARCHER COMPUTER PROGRAM" 

WHICH USES THE WAVE EQUATION 

General : This program is written in Fortran language f or the IBM 

7090 computer, but it can be used equally well with the 70 9 or 709U 
computer. It has as its objective the systematic production of data 
from which set vs. ultimate ground resistance curves can be plotted 
for any case. Its chief advantage lies in the fact that for every set 
of data run, the program automatically begins with a value of ground 
resistance of 20 tons and finds the set associated with all the load 
carried by the pile tip, 75$ of load carried by the pile tip, 50$ of 
load carried by the pile tip, 25$ of load carried by the pile tip, and 
finally all the load carried by friction along the sides of the pile. 

It then assumes an increase of ground resistance of 20 tons and repeats 
the calculations. It automatically assigns new values of ground re- 
sistance in this manner up to and including 280 tons. Essentially, 
then, it progresses through the basic program five times at each value 
of ultimate ground resistance, or 70 times for each set of data. It 
also writes out the maximum compression and tension stresses in pounds 
per square inch for each of these 70 conditions at the head, mid-length, 
and tip of the pile. The side friction distribution may be chosen as 
rectangular or triangular in shape. In both of these cases the top 10 
feet of the pile is considered free of frictional resistance. 

Terminati on ; Termination of each of the 70 conditions for each 
set of data is programmed exactly like the basic program. The reason 
for each termination is written out, and the computer immediately moves 
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to the next set of conditions. No over-all termination between sets 
of data should occur. 

Stability : Besides the checks to detect instability in the 

calculations as described under the basic program, a further check 
is obtained when plotting set vs. ultimate ground resistance. This 
plot should also produce smooth curves. If not, a set of data for 
each suspicious point can be prepared and run with the basic program 
to obtain a detailed analysis so that evaluation can be made of the 
validity of the results. Reduced time intervals can be used to elim- 
inate the instability, or smaller pile lengths in the mathematical 
model together with reduced time intervals may be used. 

Preparation of Input Data : All data is punched on standard 

cards exactly as described in the basic program except that two addi- 
tional cards are required far this program. 

The first extra card is inserted between cards 2 and 3. In 
the first 10 spaces of this new card are punched, with a decimal point, 
the cross sectional area of the pile at its top, within spaces 11 
through 20 is punched in the same manner the pile cross-sectional area 
at mid-length of the pile, and similarly the area of the pile tip in 
spaces 21 through 30. 

The second extra card is included as the very last card of 
data. If side distribution is rectangular, punch -1 in spaces 3 and 
1*; if side distribution is to be triangular, punch +1 in spaces 3 and U. 

In all other respects, data is prepared exactly as described 
under instructions for the "basic computer program." 
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Running the program : The data deck of cards just described 

is placed immediately after the program deck (after the * DATA card) 
and is then ready to be processed by a machine operator. 

Results : The input data will be repeated in the output 

information. The time cycle that the calculations were terminated 
for the first condition and the reason therefor will first appear. 

The ultimate ground resistance; the soil spring coefficients; set; 
and the maximum tension and compression at the pile head, mid-length, 
and tip will be written out for each condition. The results are 
printed in decimal form with a multiplier notation. For example, 
0.1511 E01 means 0.1511 x 10 1 = 1.5ll. 

The Program : The complete program is presented so that any 

who wish to use it may do so freely. Each line of the progran is 
punched on a separate card and must be assembled in the exact order 
given below. Extreme care must be used in punching to prevent errors 
from creeping in, and for that reason it is recommended that test 
cases be run to check the program before relying on its results. 

'The notation used is as described under this section of 
instructions for the "basic computer program." 
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C WAVE EQUATION VARY RU 

DIMENSION DS(IOO), DM(lOO), DL(100),VS(100), VM(lOO), VL(lOO) 
DIMENSION DPSM(lOO), DPSL(lOO) 

DIMENSION S(lOO). W(lOO), DP(300), R(lOO), SP(100),CS(100) 
DIMENSION CM(lOO), CL(lOO),FS(lOO), FM(lOO), FL(l00),DPSS(100) 
207 READ INPUT TAPE 5,27,NCASE,ITURN 

27 FORMAT (110/110) 

READ INPUT TAPE 5, 71,AREA1,AREA2,AREA3 
71 FORMAT ( 3F10 . U ) 

0 READ INPUT TAPE 3, 70, M,IM1,T, (W(I),I-l,M),(S(I),I«l,IMl) 

1 ,(SP(I),I= 1 M), RES1, RES2, V,Q,Z,ZP,CAP,RUB 
70 FORMAT (2IU/(F10.U)) 

READ INPUT TAPE 5,60U,MUD 
60U FORMAT (Il|) 

WRITE OUTPUT TAPE 6,30,NCASE 

30 FORMAT (12H CASE NUMBER 15) 

WRITE OUTPUT TAPE 6,28, (I,W(l),S(l),SP(l),I-l,M) 

28 FORMAT (5H DATA//58H M W S 

1 SP/(IU,3F20.2)) 

WRITE OUTPUT TAPE 6,2?,RES1,RES2,ZP,Z,Q,V 

29 FORMAT (36H COEF. OF REST. OF CAPBLOCK (RESl) = FU.2/ 

1 36H COEF. OF REST. OF PILECAP (RES2) = FU.2/ 

2 36H SIDE DAMPING FACTOR (ZP) = FU.2/ 

3 36H POINT DAMPING FACTOR (Z) = FU.2/ 

U 36H GROUND QUAKE (Q) = FU.2/ 

3 36H INITIAL RAM VELOCITY (V) = F6.2) 

WRITE OUTPUT TAPE 6.605,MUD 
605 FORMAT ( 6 H MUD = 16) 

WRITE OUTPUT TAPE 6 , 31, T, CAP, RUB 

31 FCRMAT (36H TIME INTERVAL - F 8 . 6 / 

1 36 H CAP = F6.2/ 

2 36 H RUB = F6.2//) 

CHECK = 2.* V 

RU = uoooo. 

DO 600 I = 1,1U 
ITEM = M-3 
AITEM = ITEM 
DO 602 J = 1,5 
FMAX1 = 0.0 
FMAX2 =0.0 
FMAX3 = 0.0 
FTEN 1 =0.0 
FTEN2 =0.0 
FTEN3 = 0.0 
DO 500 I = 1,100 
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DS(I) = 0.0 
DM(l) =0.0 
DL(I) =0.0 
VS(I) = 0.0 
VM(I) = 0.0 
VL(I) = 0.0 
DPSM(I) =0.0 
DPSL(I) =0.0 
R(I) = 0.C 
SP(I) = 0.0 
CS(I) = 0.0 
CM(I) = 0.0 
CL(I) = 0.0 
FS(I) = 0.0 
FM(I) = 0.0 
FL(I) = 0.0 

500 DPSS(l) =0.0 

DO 501 I = 1,300 

501 DP(I) =0.0 
DS(1) = V*12.-*T 

VS(1) = V +( -DS (1 )#S (l ) )*T*32 . 2/W( 1 ) 

VS(2) = (DS(1)«S(1))«T*32.2/W(2) 

DM(1) = DS(1) + VS(l)*12.*T 
DM(2) = VS(2)-:i-12.-:: J T 
CM( 2 ) = DM(2) 

VM(1) = VS(1) -((DM(l)-DM(2))*S(l))#T*32.2/W(l) 

VM( 2 ) = VS ( 2 ) + ( (DM( 1 )-DM(2 ) )*S (l )-(DM( 2 ) -DM( 3 ) )*S ( 2 ) )^T-x-32 . 2/W( 2 ) 
VM(3) = VS(3)+((DM(2)-DM(3))-::-S(2)-R(3))*T*32.2/W(3) 

CM(1) = DM(l) - DM(2) 

AJ = J 

PART = (5.-AJ)/U. 

SP(M) = PART-x-RU/Q 
SIDE = RU - PART-x-RU 

if(mud)65o,65i,652 

650 do 6oi n = u,im 

601 SP(N) = SIDE/ ( AITEM-xQ ) 

SP(M) = SP(M) + SIDE/(AITEMkQ) 

GO TO 65l 

652 DO 603 N = 3,im 

NT = N-2 
NT1 = 2-x-NT - 1 
IM2 = M-2 
ANTI = NT1 
AIM2 = IM2 

603 SP(N) = SIDE-x-ANTl/( ( AIM2-X-X-2 )-x-Q ) 

SP(M) = SP(M) + SIDE#(2.*AIM2-1. )/( (AIM2*»t2 )«Q ) 

651 CONTINUE 

DO 101 N = 3, 300 
I4LESS1 = M-l 
DO 130 I = 1,M 

130 DL(I) = DM(I) + VM(I )-x-12.-x-T 



DO 131 I = 1.MLESS1 

131 CL(l) = DL(I) - DL(l+l) 

13 IF(DL(M>) 12,12,11; 

12 DE = 0.0 
GO TO 2h 

ll; IF(DL(M) - Q) 16,16,18 
16 DE = 0.0 
GO TO 2h 

18 VALUE = DL(M) - Q 

IF (VALUE-DE) 22,22,20 
20 DE = VALUE 
GO TO 2h 
22 DE = DE 

2k DP(N) = DE 

R(M) = (DL(M)-DP(N ) )-x-SP(M)-x-(l. +Z-xVM(M> ) 

IF (DP(N ) ) 25,25,26 
26 IF (DP(N) - DP(N-l)) 197,197,25 
25 CONTINUE 
3 VAL = S(l)*CL(l) 

IF(CL(l)-CM(l) ) 5,5,U 
U FL(l) = VAL 

GO TO 33 

5 CL(1) = CM( 1 ) 

FL( 1 ) =VAL/(RES1**2 )- ( 1 . /(RES 1**2 )-l . )-x-S (l)*CL(l) 

33 VAL2 = S(2)*CL(2) 

IF(CL(2)-CM(2)) 35,35,31* 

3l* FL(2) = VAL2 

GO TO 37 

35 CL(2) = CM(2) 

FL(2) = VAL2/(RES2^-2)-(l./(RSS2^2)-l.)-x-S(2)-»CL(2) 

37 IF(CAP) 38,36,36 

38 IF(FL(2) ) 39,36,36 

39 FL(2) = 0.0 

36 CONTINUE 

7 IF(FL(1)) 8,133,133 

8 FL(l) = 0.0 

133 DO 132 I = 3.MLESS1 

132 FL(I) = CL(I>x-S(l) 

IF(FL(2))213,211,210 

213 IF (FTEN1-FL( 2 ) )211, 211, 2lU 

22h FTEN1 = FL(2) 

GO TO 211 

210 IF(FMAX1-FL(2))215,211,211 

215 FMAX1 = FL(2) 

211 CONTINUE 

MC = M/2 + 1 
IF(FL(MC ) )2l6, 217, 218 

216 IF (FTEN2-FL(MC) >217,217,219 

219 FTEN2 = FL(MC) 

GO TO 217 

218 IF(FMAX2-FL(MC) >220,217,217 

220 FMAX2 = FL(MC) 
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217 CONTINUE 

IF(FL(IM1))221,222,22U 

221 IF(FTEN3-FL(IM1))222,222,223 
223 FTEN3 = FL(IMI) 

GO TO 222 

22U IF(FMAX3-FL(IM1))225,222,222 

225 FMAX3 = FL(IMl) 

222 CONTINUE 

U2 IF(RUB)U9,1|9,UU 
Ub DO U 8 I = 3,MLESS1 
DPSL(I) = DPSM(I) 

CHANGE = DL(I) - Q 
SUM = DL(l) + Q 
IF (DPSL( I ) -CHANGE )U5,U6,U6 
1*5 DPSL(I) = CHANGE 

U6 if(dpsl(i)-sum)U8,U8,U7 

1*7 DPSL(I) = SUM 

U8 r(i) = (dl(i) - dpsl(i))-*sp(iMi. + zpwm(i)) 

h9 DO 50 I = 1,M 

VL(I) = VM(I) + (FL(I-l) - FL(I ) - R(l))*T*32.2/W(D 
50 CONTINUE 

DO 75 K = 1,M 
X = 1. 

IF (VL(K) ) 75,75,800 

75 CONTINUE 
GO TO 199 
800 CONTINUE 

IF(VL(2) - CHECK )99, 99, 190 
99 IF(VL(M)-CHECK)102, 102,191* 

102 DO 98 K = 1,M 
DS(K) * DM(K) 

DM(K) = DL(K) 

DL(K) = 0.0 
VS(K) = VM(K) 

VM(K) = VL(K) 

VL(K) = 0.0 
CS(K) = CM(K) 

CM(K) = CL(K) 

CL(K) = 0.0 
FS(K) = FM(K) 

FM(K) = FL(K) 

FL(K) = 0.0 
DPSS(K) = DPSM(K) 

DPSM(K) = DPSL(K) 

98 DP3L(K) =0.0 
101 CONTINUE 
GO TO 196 

190 WRITE OUTPUT TAPE 6,191,N 

191 0 FORMAT ( 6 I 4 H VELXITY OF PILE CAP EXCEEDED TWICE THE RAM VELOCITY W 
1 HEN N WAS 13) 

GO TO 196 
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19il WRITE OUTPUT TAPE 6,195, N 

195 0 FORMAT (6iiH VEIDCITY OF PILE TIP EXCEEDED TV/ICE THE RAM VELOCITY W 
1 HEN N WAS 13) 

GO TO 196 

197 WRITE OUTPUT TAPE 6,198,N 

198 FORMAT (26H DP BECAME CONSTANT AT N =13) 

GO TO 196 

195 WRITE OUTPUT TAPE 6,200,N 

200 FORMAT (52H ALL VL WERE SIMULTANEOUSLY NEGATIVE OR ZERO AT N = IH) 

196 CONTINUE 

WRITE OUTPUT TAPE 6,607,RU 
607 FORMAT (27H ULTIMATE RESISTANCE (RU) = El5.il) 

WRITE OUTPUT TAPE 6,606, (SP(L).L=1,M) 

606 FORMAT (17H NEW VALUES OF SP/(8El5.ii)) 

WRITE OUTPUT TAPE 6,305, DP(N) 

305 FORMAT (7H SET = ElO.li) 

BLOW = l./DP(N) 

WRITE OUTPUT TAPE 6,805,BL0W 
805 FORMAT (66H 

1 blowsAn= E10.3) 

FTEN1 = FTENl/AREAl 
FTEN2 = FTEN2/AREA2 
FTEN3 = FTEN3/AREA3 

WRITE OUTPUT TAPE 6,307, FTEN1 , FTEN2 , FTEN3 
307 0 FORMAT (30H MAXIMUM TENSION AT HEAD = El5.il/ 

1 30H MAXIMUM TENSION AT CENTER = El5.il/ 

2 3 OH MAXIMUM TENSION AT TIP = E15.U) 

FMAX1 = FMAXl/AREAl 

FMAX2 = FMAX2/AREA2 
FMAX3 - FMAX3/AREA3 

WRITE OUTPUT TAPE 6,306,FMAX1,FMAX2,FMAX3 

306 0 FORMAT (26H MAXIMUM FORCE AT HEAD = E15.U/ 

1 26H MAXIMUM FORCE AT CENTER = El5.il/ 

2 26H MAXIMUM FORCE AT TIP - El5.ii//) 

602 CONTINUE 

RU = RU + UOOOO. 

600 CONTINUE 

IF(ITURN)205, 205,207 
205 CALL EXIT 
END 

* DATA 
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DETAILS OP THE !l ill LEY FORMULA COMPUTER PROGRAM" 

General : This program is written in Fortran Language for the 

IBM 7090 computer, but may be used equally well with the 709 or 709U 
computer. It has as its objective the systematic production of data of 
set vs. ultimate ground resistance using the Hiley type formula. The 
program will take given pile and hammer information and beginning with 
an ultimate ground resistance of 5 tons will range in increments of 
5 tons until an ultimate ground resistance of 2lj0 tons is reached. At 
each value of ground resistance the set is computed and written out. 

The program is arranged so that any type of hammer or pile can be used. 
The coefficients for use with the Hiley formula are as recommended by 
Chellis (ref. 5) and are built into the program except as indicated 
under input data preparation section below. 

Termination : No conditions for termination are built into 

this program. If extensive use were planned for this program, an "if" 
statement should be added so that it will cause the calculations to 
terminate vdien the computed set becomes negative. 

Input Data Preparation : All data is punched on standard data 

cards in the exact order indicated below. Upon completion of punching 
the data deck, it is placed immediately after the * DATA card of the 
source program. 

CARD ORDER DESCRIPTION 

1 If another set of data follows this set for 

another case, punch +1 in spaces 9 and 10; if 
not, punch -1. Use no decimals. 



If the hammer is double acting, punch +1 in 
spaces 9 and 10; if not," punch -1. Use no 
decimals . 

If the pile is steely punch +1 in spares 9 and 
10; if wood, -1 in spaces 9 and 10; if concrete, 
punch 0 in space 10. Use no decimals. 

If triangular side distribution is desired, punch 
+1 in spaces 9 and 10; if rectangular side dis- 
tribution is desired, puirh -1 in spaces 9 and 
10. 

In the first 10 spaces punch the hammer effi- 
ciently using a decimal point . 

In the first 10 spaces punch the weight of the 
ram in pounds using a decimal point. 

In the first 10 spaces punch the coefficient 
of restitution using a decimal point. 

In the first 10 spaces punch the weight of the 
pile including the shoe and driving cap for 
drop hammers and single acting steam hammers. 
Weight of pile including weight of anvil in case 
of double acting or differential acting steam 
hammers. All in pounds and with a decimal point. 

In the first 10 spaces punch the cross sectional 
area in square inches of the pile at its head 
or the area of the capblock if one is used. Use 
a decimal point. 

Same as card 9 but for the cross sectional pile 
area at mid-length of the pile. 

Same as card 9 but for the cross sectional pile 
area at the pile tip. 

Punch in the first 10 spaces the ground quake 
in inches, using a decimal point. 

Punch in the first 10 spaces the modulus of 
elasticity in pounds per square inch using a 
decimal point. 

Punch in the first 10 spaces the pile length in 
feet using a decimal point. 
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Punch in the first 10 spaces the height of free 
fall of the ram, in inches, for drop hammers 

or 

the normal (shortest) stroke of ram in inches 
for single acting steam hammers 

or 

for double acting, differential acting steam and 
diesel hammers use the rated energy per blow in 
foot pounds as published by the manufacturers. 

Running the Program : The data deck of cards just described is 

placed immediately after the source program deck and is ready for 
delivery to a machine operator for processing. 

Results : In the output will be repeated the original input 

information, values of the ultimate ground resistance, per cent of 
load carried by the point, and the set. The results are printed in 
decimal form with a multiplier so that 0.6192 EOO, for example, means 
0.6192 x 10°, or 0.6192. 

The Program : The complete program is presented so that any who 

wish to use it may do so. Each line of the program should be punched 
on a separate card and assembled in the exact order indicated. Anyone 
preparing cards from this print-out is cautioned to verify the program 
with test cases before relying on the results obtained. 
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*ID REESE C2165 CIVIL ENGR G HILEY 7 150 

* XEQ 

* NOBIN 

* LABEL 

C HILEY FORMULA PROGRAM 

206 READ INPUT TAPE 5, 20,ITURN,IEQUIP,IMAT,MUD 

20 FORMAT (110) 

READ INPUT TAPE 5,21,E,WR,RES,WP,APH,APAV,APT,Q,ELAS,PL,H 

21 FCRMAT (F10.ii) 

IF(IMAT) 1,3,5 

1 WRITE OUTPUT TAPE 6, 2 

2 FORMAT (12H WOODEN PILE) 

GO TO 7 

3 WRITE OUTPUT TAPE 6,1* 

ll FORMAT (ll|H CONCRETE PII£) 

GO TO 7 

5 WRITE OUTPUT TAPE 6,6 

6 FORMAT (11H STEEL PILE) 

7 IF(IEQUIP) 8,12,10 

8 WRITE OUTPUT TAPE 6,9 

9 FORMAT (3ilH SINGLE ACTING OR DROP HAMMER USED) 

GO TO 12 

10 WRITE OUTPUT TAPE 6,11 

11 FORMAT (26H DOUBLE ACTING HAMMER USED) 

12 IF(MUD) 13,205,15 

13 WRITE OUTPUT TAPE 6,ll* 

111 FORMAT (33H SIDE DISTRIBUTION IS RECTANGULAR) 

GO TO 17 

15 WRITE OUTPUT TAPE 6,l6 

16 FORMAT (32H SIDE DISTRIBUTION IS TRIANGULAR) 

17 WRITE OUTPUT TAPE 6,l8,E 

18 FORMAT (13H EFFICIENCY 3 ElO.l*) 

WRITE OUTPUT TAPE 6,19,WR 

19 FORMAT (16H WEIGHT OF RAM= E10.U) 

WRITE OUTPUT TAPE 6, 22, RES 

22 FORMAT (29H COEFFICIENT OF RESTITUTION 3 ElO.ii) 

WRITE OUTPUT TAPE 6,23,WP 

23 FORMAT (13H WT OF PILE 3 ElO.ii) 

WRITE OUTPUT TAPE 6,2i|,APH 

2i| FORMAT (26H AREA PILE HEAD IN SQ IN 3 ElO.ii) 

WRITE OUTPUT TAPE 6,25,APAV 

25 FORMAT (20H AVERAGE PILE AREA 3 ElO.ii) 

WRITE OUTPUT TAPE 6, 26, APT 

26 FORMAT (16H AREA PILE TIP 3 ElO.ii) 

WRITE OUTPUT TAPE 6,27.Q 

27 FORMAT (3H QUAKE 3 E10.U) 

WRITE OUTPUT TAPE 6,28,ELAS 

28 FORMAT (17H YOUNGS MODULUS 3 ElO.ii) 

WRITE OUTPUT TAPE 6, 29, PL 

29 FORMAT (20H PILE LENGTH IN FT 3 ElO.ii) 

IF (IEQUIP) 31,35,33 
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31 WRITE OUTPUT TAPE 6,32,H 

32 FORMAT (2101 HAMMER DROP IN INCHES= ElO.U) 

GO TO 35 

33 WRITE OUTPUT TAPE 6 , 3h, H 

3h FORMAT (UOH RATED ENERGY OF HAMMER IN FOOT POUNDS= ElO.U) 

35 RU = 10000. 

DO 20U J = 1,U8 
WRITE OUTPUT TAPE 6,30,RU 
30 FORMAT (20H THE VALUE OF RU IS ElO.U//) 

PI = RU/APH 
P2 = RU/APAV 
P3 = RU/APT 
DO 203 I = 1,3 
IF(I-l) 60,U0,U2 
UO PLMOD = PL 
WPC = WP/2. 

WRITE OUTPUT TAPE 6,Ul 
Ul FORMAT (29H ALL OF LOAD CARRIED BY POINT) 

GO TO 60 

U2 IF (1-2) 60,UU,52 

UU WRITE OUTPUT TAPE 6, U7 

U7 FORMAT (31H ONE HALF LOAD CARRIED BY POINT) 

IF(MUD) U6,60,U8 
U6 PLMOD =0.75 *PL 
WPC = WP 
GO TO 60 

U8 PLMOD = U. * PL / 6. 

WPC = WP 
GO TO 60 

52 IF(I-3) 60, 5U, 60 
5U WRITE OUTPUT TAPE 6, 55 

55 FORMAT (29H ALL LOAD CARRIED BY FRICTION) 

IF (MUD) 56,60,58 

56 PLMOD = PL/2 
WPC = WP 

GO TO 60 

58 PLMOD = 2. -x- PL / 3. 

WPC = WP 
60 C3 = Q 

IF(IMAT) 70,80,90 
70 Cl- 0.00010 * PI 

C2 = 0 . 00008*P2*PLM0D*l5 00000 . /ELAS 
GO ‘TO 100 

80 Cl- 0.00025 % PI 

C2 = 0 . 00000U-*P2*PLM0D*3000000 . /ELAS 
GO TO 100 

90 Cl = 0.00008 -x- PI 

C2 = 0 . 000000U-xP2-x-PLM0D-x-30000000. /ELAS 
100 DP=( ( (E-xWR-x-H ) J x-(WR+WPC'X-RES-x-x-2 ) )/((WR+WPC )-x-RU) )-(Cl+C2+C3)/2. 

IF(IEQUIP) 10U,203,102 
102 DP = DP + (Cl+C2+C3)/2. 



150 



DP = 12.-X-DP 

DP = DP -(Cl+C2+C3)/2. 

10U BLOW = l./DP 

WRITE OUTPUT TAPE 6, 201, DP 

201 F0RMAT(8H SET IS E20.U) 

WRITE OUTPUT TAPE 6, 202, BLOW 

202 FORMAT (2 OH BLOWS PER INCH ARE E20.U/) 

203 CONTINUE 

20U RU = RU + 10000. 

IF (ITURN) 205,205,206 
205 CALL EXIT 
END 

* DATA 



- 
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APPENDIX E 

This appendix includes some of the graphs of computer solutions 
which were used both for correlation with actual pile cades and for 
studying the effect of varying the size of hammer, length of pile, 
ground quake, point damping, and side damping. On each sheet several 
curves are drawn for various values of frictional resistance expressed 
as a percentage of the total ultimate bearing capacity of the pile and 
indicated as follows: 

Curve A — End bearing; no side friction 

Curve B — 75 per cent end bearing; 25 per cent side friction. 

Curve C — 50 per cent end bearing; 50 per cent side friction. 

Curve D — 25 per cent end bearing; 75 per cent side friction. 

Curve E — 100$ friction. 

Rectangular side distribution was used. The description of the driving 
equipment, pile, and soil properties for which the computations were 
made appear on each sheet, but for convenience of use in future correla- 
tion work an index is included which sumnarizes this information for 
each sheet. The figure numbers and page numbers are identical in this 



appendix. 
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INDEX OF FIGURES IN APPENDIX E 



1. Piles Driven by a 


Vulcan #1 Hammer 
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■foot 12BP53 Piles cont'd. ) 
Point Damping Side Damping 


Ground Quake 


Figure 
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2. Piles Driven with a Vulcan 0 Hammer 

a. 12BP53 Piles 80 feet long: 

Point Damping Side Damping Ground Quake 



0.20 0.067 0.10 

0.U0 0.133 0.10 

b. 12BP53 Piles U5 feet long: 
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